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I. INTRODUCTION 

In his pioneering treatment of characteristic energy 
losses of fast electrons passing through thin metal films, 
Ritchie predicted the existence of self-sustained collec- 
tive excitations at metal surfaces It had already 
been pointed out by Pines and Bohm E 0] that the long- 
range nature of the Coulomb interaction between valence 
electrons in metals yields collective plasma oscillations 
similar to the electron-density oscillations observed by 
Tonks and Langmuir in electrical discharges in gases [4L 
thereby explaining early experiments by Ruthemann [5[ 
and Lang [6[ on the bombardment of thin metallic films 
by fast electrons. Ritchie investigated the impact of the 
film boundaries on the production of collective excita- 
tions and found that the boundary effect is to cause the 
appearance of a new lowered loss due to the excitation 
of surface collective oscillations Two years later, in 
a scries of electron energy-loss experiments Powell and 
Swan [?J demonstrated the existence of these collective 
excitations, the quanta of which Stern and Ferrell called 
the surface plasmon 0]. 

Since then, there has been a significant advance in 
both theoretical and experimental investigations of sur- 
face plasmons, which for researches in the field of con- 
densed matter and surface physics have played a key 
role in the interpretation of a great variety of exper- 
iments and the understanding of various fundamental 
properties of solids. These include the nature of Van 
der Waals forces [E [IE (HI , the classical image potential 
acting between a point classical charge and a metal sur- 
face 01 E|H El, the energy transfer in gas-surface in 



teractions [16[ , surface energi es ll7l[TR Il9j , the damping 
of surface vibrational modes 2(i"2lj|. the energy loss of 
charged particles moving outside a metal surfa ce [22l . [23| , 
and the de-excitation of adsorbed molecules (24[. Sur- 
face plasmons have also been employed in a wide spec- 
trum of studies ranging from electrochemistry [25j . wet- 
ting [26[, and biosensing [IE [H, HE, to scanning tun- 
neling microscopy [3(3], the ejection of ions from sur- 
faces 31 1 , nanoparticle growth [IE HE] , surface-plasmon 
microscop y I34l . 35 1, and surface-plasmon resonance tech- 
nology [IE [33, HE SE SH, S3 • Renewed interest in 



surface plasmons has come from recent advances in the 
investigation of the electromagnetic properties of nanos- 
tructured materials [IE S3 > onc of the most attractive 
aspects of these collective excitations now being their 
use to concentrate light in subwavelength structures and 
to enhance transmission through periodic arrays of sub- 
wavelength holes in optically thick metallic films [IE S(| . 

The so-called field of plasmonics represents an excit- 
ing new area for the application of surface and interface 
plasmons, and area in which surface-plasmon based cir- 
cuits merge the fields of photonics and electronics at the 
nanoscale [47|. Indeed, surface-plasmon polaritons can 
serve as a basis for constructing nanoscale photonic cir- 
cuits that will be able to carry optical signals and electric 
currents [IE SE • Surface plasmons can also serve as a 



basis for the design, fabrication, and characterization of 
subwavelength waveguide components [50l . l5ll . I52I l53l . H3, 
HE HE HI \M, HE WL M, HE Hi, H . In the framework 
of plasmonics, modulators and switches have also been 
investigated [gE [66| , as well as the use of surface plas- 
mons as mediators in the transfer of energy from donor to 
acceptors molecules on opposite sides of metal films fo7j . 

According to the work of Pines and Bohm, the 
quantum energy collective plasma oscillations in a free 
electron gas with equilibrium density n is hu> p = 
h(4:Trne 2 /m e ) 1 / 2 , uj p being the so-called plasmon fre- 
quency [68[ ■ In the presence of a planar boundary, there 
is a new mode (the surface plasmon), the frequency of 
which equals in the nonretarded region (where the speed 
of light can be taken to be infinitely large) Ritchie's 
frequency lj s = uj p /^/2 at wave vectors q in the range 
lu s /c << q << qF (<?f being the magnitude of the Fermi 
wave vector) and exhibits some dispersion as the wave 
vector is increased. In the retarded region, where the 
phase velocity u> s /q of the surface plasmon is comparable 
to the velocity of light, surface plasmons couple with the 
free 1 electromagnetic field. These- surface-plasmon polari- 
tons propagate along the metal surface with frequencies 
ranging from zero (at q — 0) towards the asymptotic 
value oj s = oj p /v2, the dispersion relation uu(q) lying to 
the right of the light line and the propagating vector be- 
ing, therefore, larger than that of bare light waves of the 
same energy. Hence, surface-plasmon polaritons in an 
ideal semi-infinite medium are nonradiative in nature, 
i.e., cannot decay by emitting a photon and, conversely, 
light incident on an ideal surface cannot excite surface 
plasmons. 

In the case of thin films, the electric fields of both 
surfaces interact. As a result, there are (i) tangential 
oscillations characterized by a symmetric disposition of 
charge deficiency or excess at opposing points on the two 
surfaces and (ii) normal oscillations in which an excess 
of charge density at a point on one surface is accompa- 
nied by a deficiency at the point directly across the thin 
film. The phase velocity of the tangential surface plas- 
mon is always less than the speed of light, as occurs in 
the case of a semi-infinite electron system. However, the 
phase velocity of normal oscillations may surpass that 
of light, thereby becoming a radiative surface plasmon 
that should be responsible for the emission of light [IE] • 
This radiation was detected using electron beam bom- 
bardment of thin films of Ag, Mg, and Al with thick- 
nesses ranging between 500 and 1000 A [zE0|. More re- 
cently, light emission was observed in the ultraviolet from 
a mctal-oxidc-metal tunnel diode and was attributed to 
the excitation of the radiative surface plasmon [72| . 

Nonradiative surface plasmons in both thin and thick 
films can couple to electromagnetic radiation in the pres- 
ence of surface r oug hness or a grating, as suggested by 
Teng and Stern [73j. Alternatively, prism coupling can 
be used to enhance the momentum of incident light, 
as demonstrated by Otto (74[ and by Kretchmann and 
Raether 1751 . Since then, this so-called attenuated refiec- 
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tion (ATR) method and variations upon it have been 
used by several workers in a large variety of applica- 
tions 0, [rE M, Izl, M, ED, SI . 

During the last decades, there has also been a signifi- 
cant advance in our understanding of surface plasmons in 
the nonretarded regime. Ritchie |83l | and Kanazawa (84j 
were the first to attack the problem of determining the 
dispersion uj(q) of the nonretarded surface plasmon. Ben- 
nett [HI used a hydrodynamical model with a continu- 
ous decrease of the electron density at the metal surface, 
and found that a continuous electron-density variation 
yields two collective electronic excitations: Ritchie's sur- 
face plasmon at u> ~ u> s , with a negative energy disper- 
sion at low wave vectors, and an upper surface plasmon 
at higher energies. In the direction normal to the surface, 
the distribution of Ritchie's surface plasmon consists of 
a single peak, i.e., it has a monopole character; however, 
the charge distribution of the upper mode has a node, i.e., 
it has a dipole character and is usually called multipolc 
surface plasmon. 

Bennett's qualitative conclusions were generally con- 
firmed by microscopic descriptions of the electron gas. 
On the one hand, Feibelman showed that in the long- 
wavelength limit the classical result ui s = uo p /^/2 is cor- 
rect for a semi-infinite plane-bounded electron gas, irre- 
spective of the exact variation of the electron density in 
the neighborhood of the surface [H| . On the other hand, 
explicit expressions for the linear momentum dispersion 
of the conventional monopole surface plasmon that are 
sensitive to the actual form of the electron-density fluc- 
tuation at the surface were derived by Harris and Grif- 
fin [U using the equation of motion for the Wigner dis- 
tribution function in the random-phase approximation 
(RPA) and by Flores and Garci'a-Moliner [88[ solving 
Maxwell's equations in combination with an integration 
of the field components over the surface region. Quan- 
titative RPA calculations of the linear dispersion of the 
monopole surface plasmon were carried out by several 
authors by using the infinite-barrier model (IBM) of the 
surface [83| , a step potential [90l . [9l| , and the more realis- 
tic Lang-Kohn self-consistent surface potential [9^ ]. 
Feibelman's calculations showed that for the typical elec- 
tron densities in metals (2 < r s < 6) the initial slope of 
the momentum dispersion of monopole surface plasmons 
of jellium surfaces is negative (93|, as anticipated by Ben- 
nett (H). 

Negative values of the momentum dispersion had been 
observed by high-energy electron transmission on unchar- 
acterized Mg surfaces [94[ and later by inelastic low- 
energy electron diffraction on the (100) and (111) sur- 
faces of Al jaiHi]. Nevertheless, Klos and Raether (97| 
and Krane and Raether [98[ did not observe a negative 
dispersion for Mg and Al films. Conclusive experimental 
confirmation of the negative surface plasmon dispersion 
of a variety of simple metals (Li, Na, K, Cs, A l , and Mg ) 
did not come until several years later [93, llOOl . llOll . Il02| , 
in a series of experiments based on angle -resolved low- 
energy inelastic electron scattering |l03| . These ex- 



periments showed good agreement with self-consistent 
dynamic al-re sponse calculations carried out for a jellium 
surface jl04j in a time-dependent adiabatic extension 
of the density-fu nctio nal theory (DFT) of Hohenbcrg, 
Kohn, and Sham |l05| |. Furthermore, these experiments 
also showed that the multipole surface plasmon was ob- 
servable, its energy and dispersion being in quantitative 
agreement with the self-consistent .jelli um calculations 
that had been reported by Liebsch [l06j ] . 

Significant deviations from the dispers ion of surface 
plas mons at j elliu m surfaces occur on Ag |l07l . 1 108L Il09l 
llldj ] and Hg due to the presence of filled Ad and 

5d bands, respectively, which in the case of Ag yields 
an anomalous positive dispersion. In order to describe 
the observed features of Ag surface plasmons, various 
simplified mode ls fo r the screening of d electrons have 
been developed [nj, EH, ELI El EH • Most recently, 
calculations have been found to yield a qualitative un- 
derstanding of the existing electron energy-loss measure- 
ments by combining a self-consistent jellium model for 
valence 5s electrons with a so-called dipolium model in 
which the occupied Ad bands are represented in terms of 
polarizable spheres located at t he sit es of a semi-infinite 
face-cubic-centered (fee) lattice |l 17| . 

Ab initio bulk calculations of the dynamical response 
and plasmon dispersions of noble metals w ith occu- 
pied d bands have been carried out recently |l!8t Ill9l 
Il20l |. However, first-principles calculations of the surface- 
plasmon energy and linewidth dispersion of real solids 
have been carried out only in the case of the si mple - meta l 
prototype surfaces Mg(0001) and Al(lll) [m|, H^ . 
These calculations lead to an accurate description of the 
measured surface-plasmon energy dispersion that is supe- 
rior to that obtained in the jellium model, and they show 
that the band structure is of paramount importance for 
a correct description of the surface-plasmon linewidth. 

The multipole surface plasmon, which is originated 
in the selvage electronic structure at the surface, has 
been obse rved in a varie ty of simple metals at lu ~ 
0.8w p [9l [l00|, EqU, [Toll, in agreement with theoreti- 
cal predictions. Nevertheless, electron energy-loss spec- 
troscopy (EELS) measurements of Ag, Hg, and Li re- 
vealed no clear evidence of the multipole surface plas- 
mon. In the case of Ag, high-resolution energy-loss spec- 
troscopy low-energy electron diffraction (ELS-LEED) 
measurements indicated that a peak was obtained at 
3.72 eV by su btracting the data for two different impact 
energies |123| . which was interpreted to be the Ag mul- 
tipole plasmon. However, Liebsch argued that the fre- 
quency of the Ag multipole surface plasmon should be 
in the 6 — 8 eV range above rather than below the bulk 
plasma frequency, and suggested that the observed peak 
at 3.72 eV mi ght n ot be associated with a multipolc sur- 
face plasmon (124| . 

An alternative spectroscopy technique to investigate 
multipolc surface plasmons is provided by angle - and 
energy-resolved photoyield experiments (AERPY) [l25l | . 
In fact, AERPY is more suitable than electron energy- 
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loss spectroscopy to identify the multipole surface plas- 
mon, since the monopole surface plasmon of clean flat 
surfaces (which is the dominant feature in electron- 
energy loss spectra) is not excited by photons and thus 
the weaker multipole surface mode (which intersects the 
radiation line in the retardation regime) can be observed. 
A large increase in the surf ace p hotoyield was observed 
at u) = 0.8uj p from Al(100) and Al(lll) [H|. Re- 

cently, the surface electronic structure and optic al re - 
sponse of Ag has been studied using this technique jl27l ] . 
In these experiments, the Ag multipole surface plasmon 
is observed at 3.7 eV, while no signature of the multipole 
surface plasmon is observed above the plasma frequency 
(u) p = 3.8 eV) i n dis agreement with the existing theoret- 



ical prediction [124j. Hence, further theoretical work is 
needed on the surface electronic response of Ag th at g o 
beyond the s-d polarization model described in Rcf . [1241 ] . 

Another collective electronic excitation at metal sur- 
faces is the so-called acoustic surface plasmon that has 
been predicted to exist at solid surfaces where a par- 
tially occupied quasi-two-dimensional surface-state band 
coexi sts with the underlying three-dimensional contin- 
uum [l28l . ll2°| . This new low-energy collective excitation 
exhibits linear dispersion at low wave vectors, and might 
therefore affect elect ron-h ole (e-h) and phonon dynamics 
near the Fermi level [130| . It has been demonstrated that 
it is a combination of the nonlocality of the 3D dynamical 
screening and the spill out of the 3D electron density into 
the vacuum which allows the formation of 2D electron- 
density acoustic oscillations at metal surfaces, since these 
oscillations would otherwise b e co mpletely screened by 
the surrounding 3D substrate 131|. This novel surface- 
plasmon mode has been observed recently at the (0001) 
surface of Be, showing a linear energy dispersion that 
is in very good agreement with first-principles calcula- 
tions |l32j . 

Finally, we note that metal-dielectric interfaces of ar- 
bitrary geometries also support charge density oscilla- 
tions similar to the surface plasmons characteristic of 
planar interfaces. These are localized Mie plasmons oc- 
curring at frequenc ies w hich are characteristic of the 
interface geometry |l33| . The excitation of localized 
plasmons on small particles has attracted great inter- 
est over t he years i n scanning transmissi on electron mi- 
croscopy [Hi EM, fl36l [l37l . EM, EM] and near-field 



optical spectroscopy 140| . Recently, new advances in 



structuring and manipulating on th e nanometer scale 



have rekindled interest this field [l4l( . In nanostructured 
metals and carbon-based structures, such as fullerenes 
and carbon nanotubes, localized plasmons can be ex- 
cited by light and can therefo re be easily d etected as 
pronounced optical resonances |142| . Il43l . Il44| . Further- 
more, very localized dipolc and multipole modes in the 
vicinity of highly coupled structures are resp onsible for 
surface-enhanced Raman scattering jl45l . 1 14=61 ] and other 
striki ng pr operties like, e.g., the blackness of colloidal 
silver jl47| . 

Collective electronic excitations in thin adsorbed over- 



layers, semiconductor heterostructures, and parabolic 
quantum wells have also attracted attention over the last 
years. The adsorption of thin films is important, because 
of the drastic changes that they produce in the electronic 
properties of the substrate and also bec ause of related 
phenomena such as catalytic promotion 148]; however, 
the understanding of adso rbate-induced collective exci- 
tations is still incomplete [Hi, EM EM, EH EM EM 
Il55l 1 1 56l | . The excitation spectrum of collective modes in 
semiconducto r qu a ntum wel ls has been described by sev- 
eral authors \}M, EM, EM, EM, EH1. These systems, 
which have been grown in semiconductor heterostruc- 
tures with the aid of Molecular Beam Epitaxy [162] , form 
a nearly ideal free-electron gas and have been, therefore, 
a pl ayground on which to test existing many-body theo- 
ries [161E63]. 

Major reviews on the theory of collective elec- 
tronic e xcita tions at met al su rfaces have bee n giv en by 



Ritchie [165j], Feibelman [166J, and Liebsch [167J. Ex- 
perimental reviews are also av ailab le, which focus on 
high-energy EELS experiments [l68l |. surface plas mons 
on smooth and rough surfaces and on gratings 1 1691 1 , and 
angle-resolved low-energy EELS investigations [l70l . ll7ll ]. 
An extensive review on plasmons and magnetoplasmons 
in semiconductor he tero structures has been given re- 
cently by Kushwaha |l72j |. 

This review will focus on a unified theoretical descrip- 
tion of the many-body dynamical electronic response of 
solids, which underlines the existence of various collective 
electronic excitations at metal surfaces, such as the con- 
ventional surface plasmon, multipole plasmons, and the 
acoustic surface plasmon. We also review existing cal- 
culations, experimental measurements, and some of the 
most recent applications including particle-solid inter- 
actions, scanning transmission electron microscopy, and 
surfacc-plasmon based photonics, i.e., plasmonics. 



II. SURFACE-PLASMON POLARITON: 
CLASSICAL APPROACH 

A. Semi-infinite system 

1. The surface-plasmon condition 

We consider a classical model consisting of two 
semi-infinite nonmagnetic media with local (frequency- 
dependent) dielectric functions t\ and ti separated by a 
planar interface at z = (see Fig. 1). The full set of 
Maxwell's equations in the abse nce of external sources 
can be expressed as follows (1731 ] 



V x H, 



VxE, 



1 d „ 

c at 



1 d_ 

c dt 



H; 



(2.1) 



(2.2) 
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FIG. 1: Two semi-infinite media with dielectric functions ei 
and £2 separated by a planar interface at 2 = 0. 



V-(e i E i ) = 0, 



and 



V-H ?; 



0. 



(2.3) 



(2.4) 



where the index i describes the media: i = 1 at z < 0, 
and i = 2 at z > 0. 

Solutions of Eqs. (|2.ip - (|2.4[) can generally be classified 
into s-polarized and p-polarized electromagnetic modes, 
the electric field E and the magnetic field H being paral- 
lel to the interface, respectively. For an ideal surface, if 
waves are to be formed that propagate along the interface 
there must necessarily be a component of the electric field 
normal to the surface. Hence, s-polarized surface oscil- 
lations (whose electric field E is parallel to the interface) 
do not exist; instead, we seek conditions under which a 
traveling wave with the magnetic field H parallel to the 
interface (p-polarized wave) may propagate along the sur- 
face (z = 0), with the fields tailing off into the positive 
(z > 0) and negative (z < 0) directions. Choosing the 
x-axis along the propagating direction, we write 



Ei = (E ix ,0,E iz )e 



jt) 



and 



H; 



-Ki\z\ i{qix — uji) 



(2.5) 



(2.6) 



where qi represents the magnitude of a wave vector that 
is parallel to the surface. Introducing Eqs. (|2.5p and (|2.6p 
into Eqs. (|2~1> (|2~4|) . one finds 



(2.7) 



i ki Hi = H — ei E u 
c 



The boundary conditions imply that the component 
of the electric and magnetic fields parallel to the sur- 
face must be continuous. Using Eqs. (|2.7|) and (|2.8|) . one 
writes the following system of equations: 



— Hi + —H 2 







and 



Hi 



Ho. = 0, 



(2.10) 



(2.11) 



which has a solution only if the determinant is zero, i.e., 

(2.12) 



£l £2 
— + — = 0. 

Kl K 2 



This is the surface-plasmon condition. 

From the boundary conditions also follows the conti- 
nuity of the 2D wave vector q entering Eq. (|2.9p . i.e., 
qi = q 2 = q. Hence, the surface-plasmo n co ndition 
[Eq. (|2T2p ] can also be expressed as follows jl74| . 



q(uj) = - J — . 

c V £ i + £ 2 



(2.13) 



where co/c represents the magnitude of the light wave 
vector. For a metal-dielectric interface with the dielectric 
characterized by e 2 , the solution io(q) of Eq. (|2.13p has 
slope equal to c/ y/e^ at the point q = and is a mono- 
tonic increasing function of q, which is always smaller 
than cq/y/e^ and for large q is asymptotic to the value 
given by the solution of 



(-1 + e 2 = 0. 



(2.14) 



This is the nonretarded surface-plasmon condition 
[Eq. (|2.12p with Ki = k 2 = q], which is valid as long 
as the phase velocity oj/q is much smaller than the speed 
of light. 



2. Energy dispersion 

In the case of a Prud e semi-infinite metal in vacuum, 
one has £2 = 1 and Il75l ] 



£i=l- 



uj(ui + irf) ' 



(2.15) 



r\ being a positive infinitesimal. Hence, in this case 
Eq. (j2~T3l) yields 



i k 2 H 2 = e 2 E 2x , 

c 



and 



(2.8) 



(2.9) 



q(w) = - 
c 



2 2 



(2.16) 



We have represented in Fig. 2 by solid lines the dis- 
persion relation of Eq. I|2.16p . together with the light line 
lu = cq (dotted line). The upper solid line represents the 
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FIG. 2: The solid lines represent the solutions of Eq. (|2.16|) 
with uj p — 15 eV: the dispersion of light in the solid (upper 
line) and the surface-plasmon polariton (lower line). In the 
retarded region (q < uj s /c), the surface-plasmon polariton 
dispersion curve approaches the light line uj — cq (dotted 
line). At short wave lengths (q » uj a /c), the surface-plasmon 
polariton approaches asymptotically the nonretarded surface- 
plasmon frequency lj s = u p jy/2 (dashed line). 
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FIG. 3: Schematic representation of the electromagnetic 
field associated with a surface-plasmon polariton propagating 
along a metal-dielectric interface. The field strength E; [see 
Eq. p.5[) ] decreases exponentially with the distance \z\ from 
the surface, the decay constant Ki being given by Eq. (|2.18l l. 
+ and - represent the regions with lower and higher electron 
density, respectively. 

dispersion of light in the solid. The lower solid line is the 
surface-plasmon polariton 

L0 2 (q) = u 2 p /2 + c 2 q 2 - y/<4/4 + cV , (2-17) 

which in the retarded region (where q < lo s /c) cou- 
ples with the free electromagnetic field and in the non- 
retarded limit (q >> uj s /c) yields the classical nondis- 
persive surface-plasmon frequency lu s = d p / 

We note that the wave vector q entering the dispersion 
relation of Eq. (|2.17|) (lower solid line of Fig. 2) is a 2D 
wave vector in the plane of the surface. Hence, if light 
hits the surface in an arbitrary direction the external ra- 
diation dispersion line will always lie somewhere between 
the light line c q and the vertical line, in such a way that it 
will not intersect the surface-plasmon polariton line, i.e., 
light incident on an ideal surface cannot excite surface 
plasmons. Nevertheless, there are two mechanisms that 
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FIG. 4: Attenuation length U — versus q, as obtained 

from Eq. (|2.18|l at the surface-plasmon polariton condition 
[Eq. (|2.17[) ] for a Drude metal in vacuum. e\ has been taken 
to be of the form of Eq. (|2.15[) with lo v = 15 eV and £2 has 
been set up to unity. The dotted line represents the large-q 
limit of both l\ and I2, i.e., 1/q. 



allow external radiation to be coupled to surface-plasmon 
polaritons: surface roughness or gratings, which can pro- 
vide the requisite momentum via umklapp processes [73| , 
and attenuated total reflection (ATR) which provides the 
external radiation with an imaginary wave vector in the 
direction perpendicular to the surface [7J, [75| . 



3. Skin depth 

Finally, we look at the spatial extension of the electro- 
magnetic field associated with the surface-plasmon po- 
lariton (see Fig. [3]). Introducing the surface-plasmon con- 
dition of Eq. (f27T3|) into Eq. ([279]) (with q 1 = q 2 = q), one 
finds the following expression for the surface-plasmon de- 
cay constant Ki perpendicular to the interface: 

c y ei + e 2 

which allows to define the attenuation length Zj = at 
which the electromagnetic field falls to 1/e. Fig. [4] shows 
U as a function of the magnitude q of the surface-plasmon 
polariton wave vector for a Drude metal [ei of Eq. (|2.15l) ] 
in vacuum (e% = 0). In the vacuum side of the interface, 
the attenuation length is over the wavelength involved 
(h > 1/9)7 whereas the attenuation length into the metal 
is determined at long-wavelengths (q — ► 0) by the so- 
called skin-depth. At large q [where the nonretarded 
surface-plasmon condition of Eq. (|2.14p is fulfilled], the 
skin depth is Z j ~ 1/q thereby leading to a strong concen- 
tration of the electromagnetic surface-plasmon field near 
the interface. 



B. Thin films 

Thin films are also known to support surface collective 
oscillations. For this geometry, the electromagnetic fields 
of both surfaces interact in such a way that the retarded 
surface-plasmon condition of Eq. (|2.12[) splits into two 
new conditions (we only consider nonradiative surface 
plasmons), depending on whether electrons in the two 
surfaces oscillate in phase or not. In the case of a thin 
film of thickness a and dielectric fu nctio n ei in a medium 
of dielectric function e 2 , one finds [l69j | : 



ki tanh(/cia/2) k 2 



^=0 



and 



Ki coth(/cia/2) K2 



^ = 0. 



(2.19) 



(2.20) 



Instead, if the film is surrounded by dielectric layers of 
dielectric constant eo and equal thickness t on either side, 
one finds 



Ki^tanh(Kia/2) kq 



and 



where 



with 



and 



K\v coth(Kia/2) kq 



1 - Ac~ 2K0 * 



^ = 0, 



A 



1 + Ae" 


-2re t ' 


k 2 £o - 




«2e + 


n e 2 



;o = \l q 2 - % c 2 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



Electron spectrometry measurements of the dispersion 
of the surface-plasmon polariton in oxidiz ed A l films were 
reported by Pettit, Silcox, and Vincent |l76l |. spanning 
the energy range from the short-wavelength limit where 
uj ~ all the way to the long-wavelength limit 

where lo ~ cq. The agreement between the experimen- 
tal measurements and the prediction of Eqs. (|2.21|) - (|2.25|) 
(with a Drude dielectric function for the Al film and a 
dielectric constant eo = 4 for the surrounding oxide) is 
found to be very good, as shown in Fig. [5] 

In the nonretarded regime (q » u> s /c), where n\ = 
k 2 = q, Eqs. (|2~T9|) and (|2~20|) take the form 



ei + e 2 
ei - (2 



-qa 



(2.26) 




q(A ) 



FIG. 5: Dispersion w(q) of the surface-plasmon polariton of 
an Al film of thickness a = 120 A surrounded by dielectric lay- 
ers of equal thickness t = 40 A. The solid lines represent the 
result obtained from Eqs. (|2.21[1 - (|2.25|1 w i tn £ 2 = 1, eo = 4, 
and a frequency-dependent Drude dielectric func tion ei [see 
Eq. J2J5J)] with uj p = 15 eV and 77 = 0.75 eV ^f% . The solid 
circles represent the electron spectro metry measurements re- 
ported by Petit, Silcox, and Vincent [176( | . The dashed line 
represents the nonretarded surface-plasmon frequency lo p /^/E, 
which is the solution of Eq. (|2.14p with e2 = 4 and a Drude 
dielectric function ei . The dotted line represents the light line 
uj — cq. 



which for a Drude thin slab [ei of Eq. (|2. 15[) ] in vacuum 
(e 2 = 1) yields Q: 



1/2 



(2.27) 



This e quat ion has two limiting cases, as discussed by 
Ferrell [60| . At short wavelengths [qa » 1), the sur- 
face waves become decoupled and each surface sustains 
independent oscillations at the reduced frequency tu s = 
£j p /\/2 characteristic of a semi-infinite electron gas with 
a single plane boundary. At long wavelengths (qa << 1), 
there are normal oscillations at lu p and tangential 2D os- 
cillations at 



uj 2 d = {2irnaq) 1/2 , 



(2.28) 



[17 



which were later discussed by Stern 
in artificially structured semiconductors [17' 
rece ntly i n a metallic surface-state band on a silicon sur 
face Il80l. 




and observed 
and more 



III. 



NONRETARDED SURFACE PLASMON: 
SIMPLIFIED MODELS 



The classical picture leading to the retarded Eq. (|2 . 1 2[) 
and nonretarded Eq. ()2.14p ignores both the nonlocality 
of the electronic response of the system and the micro- 
scopic spatial distribution of the electron density near 
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the surface. This microscopic effects can generally be ig- 
nored at long wavelengths where q << qp; however, as 
the excitation wavelength approaches atomic dimensions 
nonlocal effects can be important. 

As nonlocal effects can generally be ignored in the re- 
tarded region where q < u> s /c (since uj s /c « qp), here 
we focus our attention on the nonrctarded regime where 
lu s /c < q. In this regime and in the absence of external 
sources, the w-components of the time-dependent electric 
and displacement fields associated with collective oscilla- 
tions at a metal surface satisfy the quasi-static Maxwell's 
equations 



V • E(r,w) = -47rcm(r,u;), 



or, equivalcntly, 



and 



(r,w) = 47T(5n(r,w), 



V-D(r,w) 



0, 



(3.1) 



(3.2) 



(3.3) 



Sn(r, w) being the fluctuating electron density associated 
with the surface plasmon, and </>(r, ui) being the to com- 
ponent of the time-dependent scalar potential. 



A. Planar surface plasmon 

1. Classical model 

In the classical limit, we consider two semi-infinite me- 
dia with local (frequency-dependent) dielectric functions 
ei and £ 2 separated by a planar interface at z = 0, as 
in Section [II Al (see Fig. 1). In this case, the fluctuating 
electron density 5n(r, uj) corresponds to a delta-function 
sheet at z = 0: 



Sn(r,uj) = <5n(r||,a;) S(z), 



(3.4) 



where r 1 1 defines the position vector in the surface plane, 
and the displacement field D(r, uj) takes the following 
form: 



D(r,«) = 



ei E(r, to), z < 0, 
£ 2 E(r,w), z > 0. 



(3.5) 



Introducing Eq. (|3.4|) into Eq. (|3.2|) . one finds that self- 
sustained solutions of Poisson's equation take the form 



(r.w) =0 o e q - r " 



-<i\~-\ 



(3.6) 



where q is a 2D wave vector in the plane of the surface, 
and q = |q|. A combination of Eqs. (|3.3p . (|3.5[) . and 
(|3.6|) with E(r, uj) = — W(f)(r,uj) yields the nonretarded 
surface-plasmon condition of Eq. (|2.14p . i.e., 



£l + £ 2 = 0. 



(3.7) 



2. Nonlocal corrections 

Now we consider a more realistic jellium model of 
the solid surface consisting of a fixed semi-infinite uni- 
form positive background at z < plus a neutralizing 
nonuniform cloud of interacting electrons. Within this 
model, there is translational invariance in the plane of 
the surface; hence, we can define 2D Fourier transforms 
E(z; q,uj) and D(z; q,u>), the most general linear relation 
between them being 



D(z; q,0j) = J dz'e{z, z'; q, uj) ■ E(z'; q, u), 



(3.8) 



where the tensor e(z, z';q,uj) represents the dielectric 
function of the medium. 

In order to avoid an explicit calculation of e(z, z'; q, w), 
one can assume that far from the surface and at low wave 
vectors (but still in the nonretarded regime, i.e., lu s /c < 
q < qp) Eq. (|3.8[) reduces to an expression of the form of 
Eq. (1331) : 



D(z;q,uj) 



eiE(z;q,w), z < zi, 



(3.9) 



I £ 2 E(z;<7,w), z > z 2 , 
where zi << and z 2 >> 0. Eqs. (|3~2"]) and ([33)) with 



E(r, uj) = —\7(f>(r,uj) then yield the following integration 
of the field components E z and D x in terms of the po- 
tential 4>{z) at zi and z 2 [where it reduces to the classical 
potential of Eq. (|33|) ]: 



dzE z (z;q,uj) = (f)(z 2 ;q,uj) - (f)(z 1 ;q,uj) (3.10) 



and 



i I dzD x (z;q,uj) = e 2 4>(z 2 ; q,uj) - e X (j){zi;q,ijj). 

'zi 

(3.11) 

Neglecting quadratic and higher-order terms in the 
wave vector, Eqs. (|3.10[) and (|3.11|) are found to be com- 
patible under the surface-plasmon condition [88j 



ei + £ 2 
£l - £2 



q [d±(uj) - d||(w)] 



(3.12) 



d±(uj) and d\\{ui) bei ng th e so-called e?-parameters intro- 
duced by Feibelman [166J: 

/d f d 

dz z —E z (z,u)) I dz —E z (z,u) 



and 

d\\(u>) 



dzz5n(z,uj)/ / dzSn(z,uj) (3.13) 



dz z —D x (z, uj) I 
dz 



dz —D x (z,uj), (3.14) 
dz 



where E z (z,w), D x (z,uj), and 8n(z,ui) represent the 
fields and the induced density evaluated in the q — > 
limit. 
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For a Drude semi-infinite metal in vacuum [e2 = 1 
and Eq. (|2.15[) for ei], the nonretarded surface-plasmon 
condition of Eq. (|3.12j) yields the nonretarded dispersion 
relation 

Gj = Lj s {l-qRe[d±(u s )-d\\(uj s )] /2 + ...}, (3.15) 

where lj s is Ritchie's frequency: lo s = u> p /y/2. For neu- 
tral jellium surfaces, d\\ (u>) coincides with the jellium edge 
and the linear coefficient of the surface-plasmon disper- 
sion u>(q), therefore, only depends on the position d±(ui s ) 
of the centroid of the induced electron density at u> 8 [see 
Eq. (|3 . 1 3[) ] with respect to the jellium edge. 



3. Hydrodynamic approximation 

In a hydrodynamic model, the collective motion of elec- 
trons in an arbitrary inhomogeneous system is expressed 
in terms of the electron density n(r, t) and the hydrody- 
namical velocity v(r, t), which assuming irrotational flow 
we express as the gradient of a velocity potential i/j(r,i) 
such that v(r,t) = —Vil>(r,t). First of all, one writes 
the basic hydrodynamic Bloch's equations (the continu- 
ity equation and th e Ber noulli's equation) in the absence 
of external sources 1 1 8 ill : 



dt 



n(r,t)=V-[n(r,t) VV>(r,t)] 



(3.16) 



and 



|v(r, t) = \ |W(r, t)| 2 + ^1 + 0(r , t l (3.17) 



and Poisson's equation: 

V 2 (j)(r,t) = 4irn(r,t), 



(3.18) 



where G[n] is the internal kinetic energy, which is typi- 
cally approximated by the Thomas-Fermi functional 



G[n] = - (Stt 2 ) 2 / 3 [n(r,t)] 



5/3 



(3.19) 



The hydrodynamic equations [Eqs. (|3.16|) - (|3.18|) ] are 
nonlinear equations, difficult to solve. Therefore, one 
typically uses perturbation theory to expand the electron 
density and the velocity potential as follows 



and 



n(r, t) = n (r) + m(r, t) + 



4>(r,t) = + Vi(r,t) 



(3.20) 



(3.21) 



so that Eqs. (|3.16|) - (|3.18|) yield the linearized hydrody- 
namic equations 



dt 



m(r,i) = V- Mr) V^i(r,t)] 



(3.22) 



^ 1 (r,t) = [/3(r)] 2 ^M+^ 1 (r,t), 
at n Q (r) 



(3.23) 



and 



V 2 ^i(r,i) = 47rm(r,i), (3.24) 
where no(r) is the unperturbed electron density and 

/3(r) = \/l/3 [37r 2 n (r)] represents the speed of prop- 
agation of hydrodynamic disturbances in the electron sys- 
tem dn. 

We now consider a semi-infinite metal in vacuum con- 
sisting of an abrupt step of the unperturbed electron den- 
sity at the interface, which we choose to be located at 
z = 0: 

{n, z < 0, 
(3.25) 
0, 2 >0. 

Hence, within this model no(r) and /3(r) are constant at 
z < and vanish at z > 0. 

Introducing Fourier transforms, Eqs. (|3.22p - (|3.25[) 
yield the basic differential equation for the plasma nor- 
mal modes at z < 0: 



(3.26) 



V 2 (uj 2 -Lu 2 p +l3 2 V 2 )M*^) (z<0), 

and Laplace's equation at z > 0: 

V 2 0i(r,w) = O (z>0), (3.27) 

where both rii(r, w) and ^i(r,w) vanish. Furthermore, 
translational invariance in the plane of the surface allows 
to introduce the 2D Fourier transform ipi(z; q, lu), which 
according to Eq. (|3.26|) must satisfy the following equa- 
tion at z < 0: 

(-q 2 + d 2 /dz 2 ) [uj 2 -uj 2 ~ (3(-q 2 + d 2 /dz 2 )} q, u), 

(3.28) 

where q represents a 2D wave vector in the plane of the 
surface. 

Now we need to specify the boundary conditions. Rul- 
ing out exponential increase at z — > oo and noting that 
the normal component of the hydrodynamical velocity 
should vanish at the interface, for each valu e of q on e 
finds solutions to Eq. ([3~28]) with frequencies flM EH] 



> u 2 p + P 2 q 2 



and 



1 



u; 2 + (3 2 q 2 + Pqyj2u% + pq 



(3.29) 



(3.30) 



Eqs. (|3.29p and (|3.30|) represent a continuum of bulk nor- 
mal modes and a surface normal mode, respectively. At 
long wavelengths, where /3q/u p << 1 (but still in the 
nonretarded regime where u> s /c < q), Eq. (|3 . 30[) yields 
the surface-plasmon dispersion relation 

uj = lu p /V2 + (3q/2, (3.31) 

which was first derived by Ritchie |83l| using Bloch's 
equations, an d lat er by Wagner |l85| and by Ritchie 
and Marusak 186] by assuming, within a Boltzmann 
transport-equation approach, specular reflection at the 
surface. 
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B. Localized surface plasmons: classical approach 

Metal-dielectric interfaces of arbitrary geometries also 
support charge density oscillations similar to the surface 
plasmons characteristic of planar interfaces. In the long- 
wavelength (or classical) limit, in which the interface sep- 
arates two media with local (frequency-dependent) di- 
electric functions ex and 62, one writes 

D i (r,w) = e i E i (r,cj), (3.32) 

where the index i refers to the media 1 and 2 separated 
by the interface. In the case of simple geometries, such as 
spherical and cylindrical interfaces, Eqs. (|3.1[) - (|3.3p can 
be solved explicitly with the aid of Eq. (|3.32|) to find 
explicit expressions for the nonrctarded surfacc-plasmon 
condition. 



1. Simple geometries 

a. Spherical interface. In the case of a sphere of di- 
electric function ex in a host medium of dielectric func- 
tion £2, the classical (long- wavelength) planar surfacc- 
plasmon c ondi tion of Eq. (|3.7p is easily found to be re- 
placed by [l33j 

lei + (l + l)e 2 = 0, J = 1,2,..., (3.33) 

which in the case of a Drudc metal sphere [ex of 
Eq. (|2.15|) ] in vacuum (e 2 = 1) yields the Mie plasmons 
at frequencies 

Wl = ^{wTi- (3 ' 34) 

b. Cylindrical interface. In the case of an infinitely 
long cylinder of dielectric function ex in a host medium 
of dielectric function e 2 , the classical (long- wavelength) 
surfacc-plasmon condition depends on the direction of 
the electric field. For electromagnetic waves with the 
electric field normal to the interface (p-polarization), 
the corresponding long-wavelength (and nonrctarded) 
surfacc-plas mon condition c oincides with that of a planar 
surface, i.e., [l87 l [l88l . fl89l ] 

ei+e 2 = 0, (3.35) 

which for Drudc cylinders [ex of Eq. (|2. 15|) ] in vacuum 
(e2 = 1) yields the planar surfacc-plasmon frequency 

For electromagnetic waves with the electric field par- 
allel to the axis of the cylinder (s-polarization), the pres- 
ence of the interface does not modify the electric field 
and one easily finds that only the bulk mode of the host 
medium is present, i.e., one finds the plasmon condition 

e 2 = 0. (3.36) 

In some situations, instead of having one single cylinder 
in a host medium, an array of parallel cylinders may be 



present with a filling fraction /. In this case and for 
electromagnetic waves polarized along the cylinders (s- 
polarization), t he p lasmon condition of Eq. (|3.36p must 
be replaced by [l9Cl| 

/ex + (l-/)e 2 = 0, (3.37) 

which for Drude cylinders [ex of Eq. (|2.15p ] in vac- 
uum (e2 = 1) yields the reduced plasmon frequency 

2. Boundary-charge method 

In the case of more complex interfaces, a so-called 
boundary-charge method (BCM) has been used by sev- 
eral authors to determine numerically the classical (long- 
wavelength) frequencies of localized surface plasmons. In 
this approach, one first considers the w-component of the 
time-dependent surface charge density arising from the 
difference between the normal components of the electric 
fields inside and outside the surface: 

^( r > w ) = ^ [E(r,^) • n| r=r _ + E(r,w) • n| r=r+ ] , 

(3.38) 

which noting that the normal component of the displace- 
ment vector [see Eq. (|3.32[) ] must be continuous yields 
the following expression: 

a s (r, u>) = ±- E(r, u) ■ n| r=r+ , (3.39) 

4ir ei 

where n represents a unit vector in the direction perpen- 
dicular to the interface. 

An explicit expression for the normal component of the 
electric field at a point of medium 2 that is infinitely close 
to the interface (r = r + ) can be obtained with the use of 
Gauss' theorem. One finds: 

E(r,w) -n| r=r+ = -n-V0(r,w) + 27r<7 s (r,w), (3.40) 

where <p(r, uj) represents the scalar potential. In the ab- 
sence of external sources, this potential is entirely due to 
the surface charge density itself: 

0(r )W )=/dV^>l. (3.41) 

Combining Eqs. (|3.39[) - (|3.41[) . one finds the following in- 
tegral equation: 

2 7r £l±!E <7s ( r)a; )_ /> r /_L^£_. n(Ts(r ' jtj) = 0i 

ex- £2 J I r — r' I J 

(3.42) 

which describes the self-sustained oscillations of the sys- 
tem. 

The boundary-charge method has been used by sev- 
eral authors to det e rmin e the normal-mode frequen- 
cies o f a cube |l9ll . Il92j and of bodies of arbitrary 
shape [l93l . fl94] . More recent applications of this method 
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include investigations of t he surface modes of channels 
cut on planar surfa ces [1951 ] , the surface modes of coupled 
parallel wires 196], and t he e l ectro n energy loss near in- 
homogeneous dielectrics jl97l 1 198| ] . A generalization of 
this procedure that i nclud es relativistic corrections has 
been reported as well jl99| . 



3. Composite systems: effective-medium approach 

Composite systems with a large number of inter- 
faces can often be replaced by an effective homogeneous 
medium that in the long-wavelength limit is charac- 
terized by a lo cal effective dielec tric function e e //(w). 
Bergman [200l ] and Milton [20Tj showed that in the 
case of a two-component system with local (frequency- 
dependent) dielectric functions ei and e 2 and volume 
fractions / and 1—/, respectively, the long-wavelength ef- 
fective dielectric function of the system can be expressed 
as a sum of simple poles that only depend on the mi- 
crogeometry of the composite material and not on the 
dielectric functions of the components: 



WE— 



where u is the spectral variable 

u=[l- ei/e 2 ]~ , 



(3.43) 



(3.44) 



m v are depolarization factors, and B v are the strengths 
of the corresponding normal modes, which all add up to 
unity: 



5> = i- 



Similarly, 



with 



5> 



(3.45) 



(3.46) 



(3.47) 



The optical absorption and the long-wavelength energy 
loss of moving charged particles are known to be dic- 
tated by the poles of the local effective dielectric function 
teffito) and inverse dielectric function e~h(uj), respec- 
tively. If there is one single interface, these poles are 
known to coincide. 

In particular, in the case of a two-component isotropic 
system composed of identical inclusions of dielectric func- 
tion e\ in a host medium of dielectric function e 2 , the 
effective dielectric function e e //(<x>) can be obtained from 
the following relation: 



where E is the macroscopic electric field averaged over 
the composite: 



E = /E m + (1 - /)E 



(3.49) 



Ej„ and E out representing the average el ectri c field inside 
and outside the inclusions, respectively [202j |. 

a. Simple geometries. If there is only one mode with 
strength different from zero, as occurs (in the long- 
wavelength limit) in the case of one single sphere or cylin- 
der in a host medium, Eqs. (|3.43|) and (|3.46D yield 



1-/- 



1 



and 



}Au>) 



1 + /- 



(3.50) 



(3.51) 



normal modes occurring, therefore, at the frequencies 
dictated by the following conditions: 



and 



m t\ + (1 — m) e 2 = 



n ei + (1 — n) e 2 = 0. 



(3.52) 



(3.53) 



For Drude particles [ei of Eq. (|2.15[) ] in vacuum (e 2 = 1), 
these frequencies are easily found to be w = ^fmuj v and 
u) = s/ntOp, respectively. 

Indeed, for a single 3D spherical or 2D circular [203| 
inclusion in a host medium, an elementary analysis shows 
that the electric field E;„ in the interior of the inclusion 
is 



E, 



E, 



(3.54) 



where m = 1/D, D representing the dimensionality of 
the inclusions, i.e., D = 3 for spheres and D = 2 for 
cylinders. Introduction of Eq. (j3~54"]l into Eq. ([3~I5]) 
leads to an effective dielectric function of the form of 
Eq. ([3~50l) with m = 1/D, which yields [see Eq (j3~52]) ] 
the surface-plasmon condition dictated by Eq. (|3.33|) with 
I = 1 in the case of spheres (D = 3) and the surface- 
plasmon condition of Eq. (|3 .35[) in the case of cylinders 
(D = 2). This result indicates that in the nonretarded 
long-wavelength limit (which holds for wave vectors q 
such that Lu s a/c < qa << 1, a being the radius of the in- 
clusions) both the absorption of light and the energy-loss 
spectrum of a single 3D spherical or 2D circular inclusion 
exhibit one single strong maximum at the dipole reso- 
nance where e\ + 2e 2 = and t\ + e 2 = 0, respectively, 
which for a Drude sphere and cylinder [ei of Eq. (|2.15D ] 
in vacuum (e 2 = 1) yield uj = oj p /^/3 and u> = uj p /^/2. 

In the case of electromagnetic waves polarized along 
one single cylinder or array of parallel cylinders (s- 
polarization), the effective dielectric function of the com- 
posite is simply the average of the dielectric functions of 
its constituents, i.e.: 



~eff 



e 2 )E = /(ei-e 2 )E;, 



(3.48) 



e e ff =/ei + (1 -/) e 2 , 



(3.55) 
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which can also be written in the form of Eqs. (|3.50[) 
and (|3.5ip . but now with m = and n = /, respec- 
tively. Hence, for this polarization the absorption of 
light exhibits no maxima (in the case of a dielectric host 
medium with constant dielectric function 62) and the 
long-wavelength energy-loss spectrum exhibits a strong 
maximum at frequencies dictated by the plasmon condi- 
tion of Eq. (|3.37p , which in the case of Drudc cylinders 
[ei of Eq. (|2 . 1 5|) ] in vacuum (c2 — 1) yields the reduced 
plasmon frequency uj = y/J u p . 

b. Maxwell- Garnett approximation. The interaction 
among spherical (or circular) inclusions in a host medium 
can be introduced approximately in the framework of th e 
well-known Maxwell-Garnett (MG) approximation |l33l |. 

The basic assumption of this approach is that the av- 
erage electric field Ej„ within a particle located in a sys- 
tem of identical particles is related to the average field 
E OM t in the medium outside as in the case of a single 
isolated (noninteracting) particle, thereby only dipole in- 
teractions being taken into account. Hence, in this ap- 
proach the electric field E in is taken to be of the form of 
Eq. (|3.54p but with the macroscopic electric field E re- 
placed by the electric field E out outside, which together 
with Eqs. (|3~48)l and (|3~49)) yields the effective dielec- 
tric function and effective inverse dielectric function of 
Eqs. (|3.50p and (|3.5ip with the depolarization factors 
in = n = 1/D (corresponding to the dilute limit, where 
/ — * 0) replaced by 



m = — (1 - /) 
D J ' 



and 



1 

D 



(3.56) 



(3.57) 



4- Periodic structures 

Over the years, theoretical studies of the normal modes 
of complex composite systems had been generally re- 
stricted to mean-field theories of the Maxwell-Garnett 
type, which approximately acc ount for the behavior of 
localized dipole plasmons [l33| . Nevertheless, a num- 
ber of methods have been developed recently for a 
full s olution of Maxwell ' s equatio ns in periodic struc- 
tures [2^,[20i,[2^,[203,[20l[lQi]. The transfer matrix 
method has been used to determine the nor mal- mode 
frequ enci es of a lattice of metallic cylinders [21dl | and 
rods |21l| , a so-called on-shell method has been employed 
by Yannopapas et al. to investigate the plasmon modes 
of a lat tice of metallic spheres in the low filling fraction 
regime [212| . and a finite difference time domain (FDTD) 
scheme has been adapt ed t o extract the effective response 
of metallic structures |209l | . 

Most recently, an embedding method [207t ] has been 
employed to solve Maxwell's equations, which has al- 
lowed to calculate the photonic band structure of three- 
and two-dimensional lattices of nanoscalc metal spheres 




■ O 



















FIG. 6: Complementary systems in which the regions of 
plasma and vacuum are interchanged. The top panel rep- 
resents the general situation. The bottom panel represents a 
half-space filled with metal and interfaced with vacuum. The 
surface-mode frequencies u> ai and u S2 of these systems fulfill 
the sum rule of Eq. (|3.58|) . 



and c ylind ers in the frequency range of the Mie plas- 
mons [l47j . For small filling fractions, there is a surface- 
plasmon polariton which in the non-retarded region 
yields the non-dispersive Mie plasmon with frequency 
lu p /\/D. As the filling fraction increases, a continuum of 
plasmon modes is found to exist betwe en ze ro frequency 
and the bulk metal plasmon frequency 147| , which yield 
strong absorption of incident light and whose energies 
can be tu ned according to the particle-particle separa- 
tion Eil. 



5. Sum rules 

Sum rules have played a key role in providing insight 
in the investigation of a variety of physical situations. A 
useful sum rule for the surface modes in complementary 
media with arbitrary geometry was introduced by Apell 
et al. I214j . which in t he sp ecial case of a metal/ vacuum 



interface implies that [2151 ] 



(3.58) 



where lo Si is the surface-mode frequency of a given sys- 
tem, and lo S2 represents the surface mode of a second 
complementary system in which the regions of plasma 
and vacuum are interchanged (see Fig. [5]). 

For example, a half-space filled with a metal of bulk 
plasma frequency lj p and interfaced with vacuum maps 
into itself (see bottom panel of Fig. 0, and therefore 
Eq. ((3381) yields 

u Sl = uj S2 = u) p /V2, (3.59) 

which is Ritchie's frequency of plasma oscillations at a 
metal/vacuum planar interface. 

Other examples are a Drude metal sphere in vac- 
uum, which sustains localized Mie plasmons at frequen- 
cies given by Eq. (|3.34p . and a spherical void in a Drudc 
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metal, which shows Mie plasmons at frequencies 



IV. DYNAMICAL STRUCTURE FACTOR 



l + l 
21 + 1' 



(3.60) 



The squared surface-mode frequencies of the sphere 
[Eq. (dfnU)] and the voi d [(Eg . (j3~50")) ] add up to u 2 p for 
all I. as required by Eq. (|3.58[) . 

The splitting of surface modes that occurs in thin films 
due to the coupling of the electromagnetic fields in the 
two surfaces [see Eq. (|2.26[) ] also occurs in the case of 
localized modes. Apell et al. |214| proved a second sum 
rule, which relates the surface modes corresponding to 
the in-phase and out-of-phasc linear combinations of the 
screening charge densities at the interfaces. In the case 
of metal/vacuum interfaces this sum rule takes the form 
of Eq. (|3.58[) . but now uj Si and uj S2 being in-phase and 
out-of-phase modes of the same system. 

For a Drude metal film with equal and abrupt planar 
surfaces, the actual values of the nonretarded lu Si and uj S2 
are those given by Eq. (|2.27p . which fulfill the sum rule 
dictated by Eq. (|3 . 58[) . For a spherical fullcrcne molecule 
described by assigning a Drude dielectric function to ev- 
ery point between the inner and outer surfaces of radii 
T\ and r2, one finds the following freq uenci es for the in- 
phase and out-of-phase surface modes [2 id ]: 



1± 



yi + 4Z(Z + D(ri/r 2 ) 



2i+l 



21 + 1 

(3.61) 

also fulfilling the sum rule of Eq. (|3.58p . 

Another sum rule has been reported recently |210Ll21l| . 
which relates the frequencies of the modes that can be 
excited by light [as dictated by the poles of the effec- 
tive dielectric function of Eq. (|3.43p ] and those modes 
that can be excited by moving charged particles [as dic- 
tated by the poles of the effective inverse dielectric func- 
tion of Eq. (|3.46p ]. Numerical calculations for various 
geometries have shown that the depolarization factors 
m v and n v e nter ing Eqs. (|3.43p and (|3.46|) satisfy the 
relation [111, HO] 



n v = 1 — (D — 1) m v 



(3.62) 



where D represents the dimensionality of the inclusions. 

Furthermore, combining Eqs. (|3.50|) and (|3.51[) (and 
assuming, therefore, that only dipole interactions are 
present) with the sum rule of Eq. (|3.62[) yields Eqs. (|3.56p 
and (|3.57[) . i.e., the MG approximation. Conversely, as 
long as multipolar modes contribute to the spectral rep- 
resentation of the effective response [see Eqs. (|3.43[) and 
()3.46|) ]. the strength of the dipolar modes decreases [see 
Eqs. (pT45)) and (|3~47) ] and a combination of Eqs. ([3743)) 
and (|3.46[) with Eq. (|3.62p leads to the conclusion that the 
dipolar resonances must necessarily deviate from their 
MG counterparts dictated by Eqs. (|3.56p and (|3.57|) . 
That a nonvanishing contribution from multipolar modes 
appears together with a deviation of the frequencies of 
the dipolar modes with respe ct to their MG counterparts 
was shown explicitly in Ref . [2ldj . 



The dynamical structure factor S(r,r';uj) represents a 
key quantity in the description of both single-particle and 
colle ctive electronic excitations in a many-electron sys- 
tem [2171 ] . The rate for the generation of electronic exci- 
tations by an external potential, the inelastic differential 
cross section for external particles to scatter in a given di- 
rection, the inelastic lifetime of excited hot electrons, the 
so-called stopping power of a many-electron system for 
moving charged particles, and the ground-state energy 
of an arbitrary many-electron system (which is involved 
in, e.g., the surface energy and the understanding of Van 
der Waals interactions) are all related to the dynamical 
structure factor of the system. 

The dynamical structure factor, which accounts for the 
particle-density fluctuations of the system, is defined as 
follows 

S(v,v';lo) S Pon( r i) S Pno(r2)S{uj~E n + E ). (4.1) 



Here, Sp n o(r) represent matrix elements, taken between 
the many-particle ground state |^o) of energy Eq and 
the many-particle excited state \^f n ) of energy E n , of the 
operator p(r) — no(r), where p(r) is the electron-density 
operator (2li| 



JV 



P(r) = ^<5(r-r,), 



(4.2) 



with S and describing the Dirac-delta operator and 
electron coordinates, respectively, and no(r) represents 
the ground-state electron density, i.e.: 



no(r) =< *o|p(r)|* > 



(4.3) 



The many-body ground and excited states of a many- 
electron system arc unknown and the dynamical struc- 
ture factor is, therefore, difficult to calculate. Never- 
theless, one can use the zero-t empe rature limit of the 
fluctuation-dissipation theorem 219], which relates the 
dynamical structure factor 5(r, r';w) to the dynami- 
cal density-response function x(r, r'; ui) of linear-response 
theory. One writes, 

S(r, r'; u) = -- Im X (r, r'; u) 0(u), (4.4) 

7T 

where f2 represents the normalization volume and 9{x) is 
the Heaviside step function. 



V. DENSITY-RESPONSE FUNCTION 

Take a system of N interacting electrons exposed to a 
frequency-dependent external potential <jf xt (r, to). Keep- 
ing terms of first order in the external perturbation and 
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neglecting retardation effects, time-dependent perturba- 
tion theory yield s the following expression for the induced 
electron density [218l |: 



5n(T,u)= / dr'x(r,r';w)0 ext (r», 



(5.1) 



where x(r, r'; u>) represents the so-called density-response 
function of the many-electron system: 



X(r,r';w) = ^ p* (r)p ra0 (r') 

n 

1 



i 



E + E n + h{uj + irj) 



E - E n + H(lu + IT]) 
(5.2) 



77 being a positive infinitesimal. 

The imaginary part of the true density-response func- 
tion of Eq. (|5.2[) . which accounts for the creation of both 
collective and single-particle excitations in the many- 
electron system, is known to satisfy the so-called /-sum 
rule: 

/oo 
fWImx(r,r';w) = -ttV-V [n (r)5(r, r')] , (5.3) 
-00 

with 7io( r ) being the unperturbed ground-state electron 
density of Eq. (|4~3j) . 



where x°( r ! r 'i w ) denotes the density-response function 
of nonintcracting Hartrce electrons: 



X °(r,r» 



1,3 



^(r^r^rQ^r') 

u — Sj + Ei + irj 



(5.7) 



Here, /, are Fermi-Dirac occupation factors, which at 
zero temperature take the form fi = 6(sf — £f being 
the Fermi energy, and the single-particle states and ener- 
gies ^i(r) and £j are the eigenfunctions and eigenvalues 
of a Hartrce Hamiltonian, i.e.: 



1 , 
--V 2 



■v H [n ](r) 



ipi(r) = Eiipiir), 



where 



(5.8) 



(5.9) 



with vo(r) denoting a static external potential and no(r) 
being the unperturbed Hartree electron density: 



UH-[n ](r) = v (r) + / dr' v(r, r')n (r'), 



N 



>l 



w = Ei^wi 2 - 



(5.10) 



A. Random-phase approximation (RPA) 

In the so-called random-phase or, equivalently, time- 
dependent Hartrce approximation, the electron density 
5n{r,u>) induced in an interacting electron system by a 
small external potential cj> ext (r, 10) is obtained as the elec- 
tron density induced in a noninteracting Hartree system 
(of electrons moving in a self-consistent Hartree poten- 
tial) by both the external potential (f> ext (r,u>) and the 
induced potential 

6<f> H (r,u)= I dr' v(t,t') Sn(r',u}), (5.4) 



with w(r, r') representing the bare Coulomb interaction. 
Hence, in this approximation one writes 



5n{v,uj) = J dr'x°(r,r'; W ) 

b ext (r',Lu)+ f dr"v(r',r")5n(r",u>) 



(5.5) 



which together with Eq. (|5.ip yields the following Dyson- 
type equation for the interacting density-response func- 
tion: 



X(r,r';w) = x°(r,r';w)+ / dr 1 J dr 2 x° (r, ri ; w) 

x u(ri,r 2 )x(r 2 ,r';u>), (5.6) 



B. Time-dependent density-functional theory 

In the framewor k of time-dependent density-functional 
theory (TDDFT) [220j |. the exact density-response func- 
tion of an interacting many-electron system is found to 
obey the following Dyson-type equation: 

X(r,r';w) = x°(r,r';o;) + J dr 1 J dr 2 x°(r,ri;w) 

x {«(n,r 2 ) + /scM(ri,r 2 ;w)}x(i"2,r';a;). (5.11) 

Here, the noninteracting density-response function 
X°(r, r'; ui) is of the form of Eq. (|5.7p . but with the single- 
particle states and energies ipi( r ) an d £i being now the 
eigenfunctions and eigenvalues of the Kohn-Sham Hamil- 
tonian of DFT, i.e: 



1 9 
"2 V 



VKs[no](r) 



ipi(r) = eiipi(r), 



where 



VKs[no](r) = v H [no\(r) + v xc [n ](r), 



with 



v xc [n ] (r) 



SE r 



5n{r) 



(5.12) 



(5.13) 



(5.14) 



^c[ n ] represents the unknown XC energy functional and 
no(r) denotes the exact unperturbed electron density of 
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Eq. (|4.3|) , which DFT shows to coincide with that of 
Eq. (|5.10p but with the Hartree eigenfunctions ipi(r) of 
Eq. (|5.8j) being replaced by their Kohn-Sham counter- 
parts of Eq. (|5.12[) . The xc kernel f xc [no](r, r'; u>) denotes 
the Fourier transform of 



fxc[no](r,t;r';t') = 



Sv xc [n](r,t) 



(5.15) 



with v xc [n](r,t) being the exact time-dependent xc po- 
tential of TDDFT. 

If short-range XC effects are ignored altogether by set- 
ting the unknown XC potential i^c [no] (r) and XC ker- 
nel fxc[no](r,r';uj) equal to zero, the TDDFT density- 
response function of Eq. (|5 . reduces to the RPA 
Eq. (1531) . 



1. The XC kernel 

Along the years, several approximations have been 
used to evaluate the unknown XC kernel of Eq. (|5.15[) . 

a. Random-phase approximation (RPA). Nowa- 
days, one usually refers to the RPA as the result of 
simply setting the XC kernel f xc [no](r,r';uj) equal to 
zero: 



/£ p > ](r,r';u;) = 0, 



(5.16) 



but still using in Eqs. (|5.7p and (|5. 10|) the full single- 
particle states and energies tpi{r) and Si of DFT [i.e., the 
solutions of Eq. (|5. 12[) ] with v xc [?io](r) set different from 
zero. This is sometimes called DFT-bascd RPA. 

b. Adiabatic local-density approximation (ALDA). 
In this approximation, also called tim e-dependent local- 
density approximation (TDLDA) [22 1| , one assumes that 
both the unperturbed no(r) and the induced Sn(r,u>) 
electron densities vary slowly in space and time and, 
therefore, one replaces the dynamical XC kernel by the 
long-wavelength (Q — > 0) limit of the static XC kernel of 
a homogeneous electron gas at the local density: 



j^> ](r,r» 



ns.. 



=(«)] 



dn 2 



5{r 



n=n (r) 



(5.17) 

where e xc {n) is the XC energy per particle of a homoge- 
neous electron gas of density n. 

c. PGG and BPG. In the spirit of the optimized 
effective-potential meth od [222], Pctcrsilka, Gossmann, 
and Gross (PGG) [223j ] derived the following frequency- 
independent exchange-only approximation for inhomogc- 
neous systems: 



A|2 



/f GG [n ](r,r';c) = - 



IEi/i^W;(r')l 



|r-r'| 7io(r)n (r') 

(5.18) 

where ipi(r) denote the solutions of the Kohn-Sham 
Eq. (I5~T21 . 



More rec ently, Burke, Pctcrsilka, and Gross 
(BPG) [H| devised a hybrid formula for the XC 
kernel, which combines expressions for symmetric and 
antisymmetric spin orientations from the exchange-only 
PGG scheme and the ALDA. For an unpolarized 
many-electron system, one writes |224l |: 



M(r,r' 



If. 



TT, PGG 



1l,LDA~\ 



(5.19) 



where /JJ and represent the XC kernel for electrons 
with parallel and antiparallcl spin, respectively. 

d. Average approximation The investigation of 
short-range XC effects in solids has been focused in a 
great extent onto the simplest possible many-electron 
system, which is the homogeneous electron gas. Hence, 
recent attempts to account for XC effects in inhomo- 
gene ous system s have adopted the following approxima- 
tion [225l.l22i: 



/^N(r,r'; W ) = / a *r(n 



(5.20) 



where h represents a function of the electron densities at 
points r and r', typically the arithmetical average 



h= \ ["o(r) +no(r')] , 



(5.21) 



and f x ° m (n; l r — r 'h ti; ) denotes the XC kernel of a ho- 
mogeneous electron gas of density n, whose 3D Fourier 
transform f^° m (h;Q,uj) is directly connected to the so- 
called local-field factor G(n;Q,uj): 



4tt 



f2r(n;Q 7 u) = -—G(n; 



Q 2 



(5.22) 



In the ALDA, one writes 



G 



ALDA 



(n;Q,cj) 



G(h; Q -> 0,lu = 0) 

_Cf_ d 2 [ne xc (n)} 
An 



dn 2 



,(5.23) 



which in combination with Eqs. (|5.20p - (|5.22p yields the 
ALDA XC kernel of Eq. (|5.17|) . However, more accu- 
rate nonlocal dynamical expressions for the local-field 
factor G(h;Q,oj) are available nowdays, which together 
with Eqs. (|5.20[) - (|5.22p should yield an accurate (beyond 
the ALDA) representation of the XC kernel of inhomo- 
geneous systems. 

During the last decades, much effort has gone 
into the determination of the static local-field factor 

QstaUc^ Q) Q^. Q ; w = 0) SH, SU, [HO, HI 

12321 |233| . 1 . the most recent work s includi ng diffu- 
sion Monte Carlo (DMC) calculations (2351. l236| and the 
parametrization of th e DMC data of Ref . |236| given by 
Corradini et al. |237|: 



G s 



\h; Q) = CQ 2 + BQ 2 /(g + Q 2 ) + a Q 4 e 



(5.24) 
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where Q = Q/qF, and the parameters B, C, g, a, and 
are the dimensionless functions of n listed in Ref. [237l |. 

Calculations of the frequency dependence of the local- 
field factor G(h;Q 1 uj) have been carried out mainly in 
the limit of long wavelengths (Q -> 0) [HI, SH, [mJ 
[24ll . l24l [24l. b ut work has also been done for finite 
wave vectors 1244 HH US [247 



VI. INVERSE DIELECTRIC FUNCTION 

In the presence of a many-electron system, the total 
potential cj>(r, oj) of a unit test charge at point r that is 
exposed to the external potential 4> ext (r,u>) can be ex- 
pressed in the following form: 

</>(r,w) = ext (r,oj)+5<j) H (r,oj), (6.1) 

where 5(f> H (r,oj) represents the induced potential of 
Eq. (pT4)) . Using Eqs. (pTTj) and (jO)) , the total potential 
</>(r, oj) of Eq. (|6.1|) is easily found to take the following 
form: 



)= / dr'e-^r.r';^ ^(r.u,), 



(6.2) 



where 



(r, v'-uj) = 5{r -r')+ / dr" v(r - r") X (r", r'; uj). 



(6.3) 

This is the so-called inverse longitudinal dielectric func- 
tion of the many-electron system, whose poles dictate 
the occurrence of collective electronic excitations and 
which can be evaluated in the RPA or in the framework 
of TDDFT from the knowledge of the density-response 
function of Eqs. (|5.6p and (|5.11|) . respectively. 

Some quantities, such as the optical absorption and 
the electron energy loss of charged particles moving in 
arbitrary inhomogencous media, can be described by the 
so-called effective inverse dielectric function e~y^(Q,cj), 
which is defined as a 3D Fourier transform of the inverse 
dielectric function e _1 (r, r';oj): 

e-} f (Q,w) = ^ Jdv J dr'c- i Q'( r - r ')e- 1 (r,r';c), 

(6.4) 

and which at long wave lengths (Q — » 0) should take the 
form of Eq. (I3"^D [H|. 

In particular, in the case of a homogeneous system and 
in the classical long- wavelength limit, where the total 
potential </>(r, oj) of a unit test charge at point r only 
depends on the external potential ip ext (r, oj) at that point, 
the inverse dielectric function takes the following form: 

e- 1 (r,r';o;) = e- 1 H5(r-r'), (6.5) 

which in combination with Eq. (|6.2[) yields the classical 
formula 



(r,oj) = ^ xt {r,oj)/e{oj), 



(6.6) 



e(u>) representing the so-called local dielectric function of 
the medium. 



VII. SCREENED INTERACTION 

Another key quantity in the description of electronic 
excitations in a many-electron system, which also dic- 
tates the occurrence of collective electronic excitations, 
is the frequency-dependent complex screened interac- 
tion W(r, r'; oj). This quantity yields the total potential 
4>(r, oj) of a unit test charge at point r in the presence of 
an external test charge of density n ext (r',uj) at point r': 

0(r,w) = J dr'W(r,r';oj)n ext (r',oj). (7.1) 

The potential (jf xt (r, ui) due to the external test charge 
density n ext (r,uj) is simply 

c/} ext {r,Lu) = J dr' v(r,r')n ext (r',u;). (7.2) 

Hence, a comparison of Eqs. (|6.2|) and (|7.1|) yields 

W(t,t';uj) = j dr" e'^r, r"; ui) v(r", r'), (7.3) 

and using Eq. (|6.3p : 

W(r, r';oj) = u(r,r')+ J dri j dr 2 

x u(r,ri)x(ri,r 2 ,w)v(r 2 ,r'). (7.4) 

From Eqs. (|6 .4[) and (|7.3p one easily finds the following 
representation of the effective inverse dielectric function: 

e~f f (Q,uj) = dr J dr' W(r, r'; lj), 

(7.5) 

where vq = 4tt/Q 2 denotes the 3D Fourier transform of 
the bare Coulomb interaction v(r, r'). 



A. Classical model 

In a classical model consisting of two homogeneous me- 
dia characterized by local (frequency-dependent) dielec- 
tric functions ei and e 2 and separated by an interface of 
arbitrary geometry, the total potential at each medium 
is simply given by Eq. (|6.6|) and is, therefore, a solution 
of Poisson's equation 

_o , , , 47T 



e»(w) 



n ex (r, uj), 



(7.6) 



ei(u>) being ei or e 2 depending on whether the point r is 
located in medium 1 or in medium 2, respectively. Hence, 
the screened interaction W(r,r' ;oj) entering Eq. (|7.ip is 
a solution of the following equation: 



V 2 VF(r,r';w) 



47T 



S(r-r'). 



(7.7) 



For simple geometries, such as the planar, spherical, 
and cylindrical interfaces, Eq. (|7.7p can be solved ex- 
plicitly by imposing the ordinary boundary conditions of 
continuity of the potential and the normal component of 
the displacement vector at the interface. 
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1. Planar surface 

In the case of two semi-infinite media with local 
(frequency-dependent) dielectric functions ei (at z < 0) 
and e 2 (at z > 0) separated by a planar interface at z = 
(see Fig. [1]) , there is translational invariancc in two direc- 
tions, which we take to be normal to the z axis. Hence, 



one can define the Fourier transform W(z, z'; q,ui), q be- 
ing the magnitude of a 2D wave vector in the plane of 
the interface, and imposing the ordinary boundary condi- 
tions of continuity of the potential and the normal com- 
ponent of the displacement vector at the interface, one 
finds: 



2tt 

W{z,z';q,uj) = — < 

q 



e -«l*-*'l + ge -i(\ z \+\ z '\) /a, z<0,z'<0, 
2 ff e-9l 2 - 2 'l/(ei - e 2 ), z<<0,z>>0, 

/e 2 , z > 0, z' > 0, 
I 



e -q\z-z'\ _ g e -q(\z\ + \z'\) 



(7.8) 



where z < (z > ) is the smallest (largest) of z and z', and 
g is the classical surface-response function: 



ei(u>) + e 2 (u>)'' 



or, equivalcntly, 



(7.9) 



(7.10) 



where u is the spectral variable of Eq. (|3.44p and n= 1/2. 

An inspection of Eqs. (|7.8p and (|7.9[) shows that the 
screened interaction W(z, z';q,cu) has poles at the clas- 
sical bulk- and surface-plasmon conditions dictated by 
€i = and by Eq. (|3.7p . respectively. 



2. Spheres 



In the case of a sphere of radius a and local 
(frequency-dependent) dielectric function ei embedded 
in a host medium of local (frequency-dependent) dielec- 
tric function e 2 , we first expand the screened interaction 
W(t,t';lj) in spherical harmonics: 



W(v, r'; uj) = Wi{r, r'; d) Y^Cl) H >m (fi'), 

(7.11) 

and we then derive the coefficients of this ex pansi on by 
imposing the boundary conditions. One finds [249j : 



{r <y 



(r r') 1 

n 2l + l 



- + {l + l)gi 



[r>) l + 

( r <)l 



/ei, r,r'<a, 



(r>)'+ 



,2/+l I 



i - 1 91 



i r r i\l+i 
I 



r < < a, r> > a, 



(7.12) 



/e 2 , r,r' > a, 



where r< (r>) is the smallest (largest) of r and r', and 

ex(w) - e 2 (w) 



5zM 



Zei(w) + (/ + l)e 2 (w)' 



(7.13) 



or, equivalcntly, 



u — ni 



(7.14) 



with it being the spectral variable of Eq. <|3.44[) and 

I 



21 + 1 



(7.15) 



As in the case of the planar surface, the screened in- 
teraction of Eqs. (|7.11[) - (|7.13p has poles at the classical 
bulk- and surface-plasmon conditions, which in the case 
of a single sphere in a host medium are dictated by = 
and by Eq. (|3.33[) . respectively. 
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Introducing Eqs. I|7.11[) - (|7.13[) into Eq. (|7.5p . one finds 
the follo wing expression for the effective inverse dielectric 
function Il88ll: 



£e 7 1 / (Q,c) = e 2 - 1 + /( e r 1 -e 2 - 1 ) 



l + -J2(2l + l)g l j l (x) e,(a 



(=0 



where 



e z (x) 



' jl-x (x) ei — (/ + 1) e 2 



£l - £2 



(7.16) 



(7.17) 



and x = Qa. Here, / represents the volume fraction filled 
by the sphere and 7'; (x ) are spherical Bessel functions of 
the first kind |25C| . This equation represents the dilute 
(/ — > 0) limit of the effective inverse dielectric function 
derived by Barrera and Fuchs for a system co mpose d by 
identical interacting spheres in a host medium [25l|. 

In the limit as Qa « 1, an expansion of Eq. (|7.16p 
yields 



-eff 



; (Q,w) 



1-3/ 



ei 



2e 2 



(7.18) 



which is precisely the long-wavelength effective in- 
verse dielectric function obtained in Section IIII B 31 from 
Eqs. (pTl8")) and ([3~M]) with D = 3 and which admits the 
spectral representation of Eq. (|3.5ip with n = 1/3. This 
result demonstrates the expected result that in the limit 
as Qa << 1 a broad beam of charged particles interact- 
ing with a single sphere of dielectric function ei in a host 
medium of dielectric function e 2 can only create collec- 
tive excitations at the dipole resonance where t\ + 2e 2 = 



[Eq. (|3.33|) with 2 = 1], which for a Drudc sphere in vac- 
uum yields ui = lu p /\/3. 



3. Cylinders 

In the case of an infinitely long cylinder of radius a and 
local (frequency-dependent) dielectric function a embed- 
ded in a host medium of local (frequency-dependent) 
dielectric function e 2 , we expand the screened interac- 
tion W(r, r'; uj) in term s of the modified Bessel functions 
I m (x) and K m (x) [25(J |. as follows 

2 f°° 

W(r,r';uj) = - dq z cos [q z (z - z')] 

* Jo 

OO 

x fi"mW m (p,p';u) cos [m((j>- $)] , (7.19) 

m=0 

where z and p represent the projections of the position 
vector along the axis of the cylinder and in a plane per- 
pendicular to the cylinder, respectively, q z denotes the 
magnitude of a wave vector along the axis of the cylin- 
der, and m m are Neumann numbers 

r 1, m = 0, 
Mm = { (7.20) 
I 2, m>l. 

The coefficients W m (p, p'\ uj) are then derived by impos- 
ing the boundary conditions, i.e., by requiring that the 
total scalar potential and the normal component of the 
displacement vectors be continuous at the interface. One 
finds: 



* c ) 9m Im {lzP > )/I' m {x)\ /ei, p,p'<a, 



W m (p,p';uj) = I [xI' m {x)K m (x)] 1 g m I m {q z p < )K m (q z p > )/{ei- e 2 ), 



p < < a,p > > a, 



[I m (q z p < ) - I m (x) g m K m (q z p < ) / K m (x)] K m (q z p > )/e 2 , p,p' > a, 

I 



(7.21) 



where x = q z a and 



I' m (x)K m (x) [eriuj) - c { uj)] 



I^(x) K m {x) ei(u) - I m (x) K' m {x) e 2 (w) ' 

(7.22) 

or, equivalcntly, 



n, 



(7.23) 



u - n m 

with u being the spectral variable of Eq. (|3.44p and 

n m =xI' m {x)K m {x). (7.24) 



In the limit as x — > (p-polarization) [252j |. the de- 
polarization factors of Eq. (|7.24|) are easily found to 

be hq = (corresponding to the plasmon condition 

e 2 = 0) [253| and n m = 1/2 (corresponding to the planar 
surfacc-plasmon condition e\ + e 2 = 0) for all m =/= 0; 
in the limit as a; — > 00, Eq. (|7.24[) yields n m = 1/2 
for all m. For the behavior of the depolarization fac- 
tors n m of Eq. (|7.24[) as a function of x sec Fig. [7] This 
figure shows that the energies of all modes are rather 
close to the planar surface-plasmon energy (correspond- 
ing to n m = 1/2), except for m = 0. The m = mode, 
which corresponds to a homogeneous charge distribution 
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FIG. 7: Depolarization factors n m — x I' m {x) K m (x), as a 
function of x = q z a, for m = 0,1,2,3,4,5 (thin solid lines), 
and m — 10 (thick solid line). As m — > oo, the depolarization 
factor n m equals the planar surface-plasmon value n m = 1/2 
for all values of x = q z a. 



around the cylindrical surface, shifts downwards from the 
planar surface-plasmon energy (no = 1/2) as the adi- 
mcnsional quantity x — q z a decreases, as occurs with the 
symmetric low-energy mode in thin films [see Eqs. (|2.26[) 

and f2~27l) ]. 

Introducing Eqs. f7TTg|) - ([7^2|) into Eq. (fT3|) . one finds 
the follo wing expression for the effective inverse dielectric 
function I188ll: 



^ 1 f (Q,^) = ^ 1 + /(er 1 -^ 1 ) 

2 



1 



Hm.J m (y)g m <d m {x,y) 



m=0 



(7.25) 



where 



4(x) (*, y) ei + K'Jx) f { I ] (x, y) e 2 



r m (x)K m (x) [ei - e 2 ] 



(7.26) 



fm(x,y) = x J m (y) K m -i(x) +yJ m -i(y)K m (x), 

(7.27) 

fmi x ,y) =xJ m {y)Im-i(x)-yJ m -i{y)I m (x), (7.28) 

x = q z a, and y = qa, q z and g representing the compo- 
nents of the total wave vector Q along the axis of the 
cylinder and in a plane perpendicular to the cylinder, re- 
spectively. The volume fraction filled by the cylinder is 



denoted by /, an d Jm (x) are cylindrical Bcsscl functions 
of the first kind [250l |. A spectral representation of the 
effective inverse diele ctric function of Eqs. ()7.25j) - (|7.28j) 
was reported in Ref. jl89| . 

In the limit as Qa << 1, an expansion of Eqs. (|7.25|) - 
(TT281) yields 



(Q,w) 



1-/ 



ei - £2 



— + 2- 

£2 



y 



ei + £2 



or + 

(7.29) 

which admits the spectral representation of Eq. (|3.46p 
with two nonvanishing spectra strengths: Co = x 2 / (x 2 + 
y 2 ) and C\ = y 2 j{x 2 + y 2 ), the corresponding depolar- 
ization factors being no = and n\ = 1/2, respectively. 

Equation (|7.29[) demonstrates that in the limit as 
Qa « 1 and for a wave vector normal to the cylinder 
[x = 0), moving charged particles can only create collec- 
tive excitations at the dipolc resonance where ni = 1/2, 
i.e., ei +£2 = 0, which for a Drudc cylinder in vacuum 
yields Ritchie's frequency uo s = w p /y2. Conversely, still 
in the limit as Qa « 1 but for a wave vector along the 
axis of the cylinder (y — 0), moving charged particles 
can only excite the bulk mode of the host medium dic- 
tated by the condition u = (corresponding to n = 0), 
i.e., e2 = 0, in agreement with the discussion of Sec- 
tion mni] 



B. Nonlocal models: planar surface 

Nonlocal effects that are absent in the classical model 
described above can be incorporated in a variety of semi- 
classical and quantal approaches, which we here only de- 
scribe for a planar surface. 



1. Hydrodynamic model 

a. Semiclassical hydrodynamic approach. Within a 
semiclassical hydrodynamic approach, the screened in- 
teraction W(r, r';cj) [as defined in Eq. (|7.1[) ] can be ob- 
tained from the linearized hydrodynamic Eqs. (|3.22[) - 
(|3.24p . For a semi- infinite metal in vacuum consisting of 
an abrupt step of the unperturbed electron density no(z) 
[see Eq. (|3.25|) ]. we can assume translational invariancc 
in the plane of the surface, and noting that the normal 
component of the hydrodynamical velocity should van- 
ish at the interface Eqs. (|3.22j) - (|3.24j) yield the following 
expression for the 2D Fourier transform W(z, z'\ q, u>): 
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W(z,z';q,w) 



2tt 

q 



e s (z - z') + e s {z + z')-2g ^J'^ , z<0,z'<0, 



I e 



1 - £ 

— q\z — z 







fle -g(«+«') ) 
I 



z< < 0,z> > 0, 
z > 0,z' > 0, 



(7.30) 



where z K (z > ) is the smallest (largest) of z and z', 



^s{z;q,oj) 



Auj(uj + irj) e~ q ^ - qu; 2 e~ A|z| 
A [lu(uj + irj) — w 2 ] 



2/? 2 A(A + g) - u}% ' 



A "5 



Ju 2 +P 2 q 2 — + irj), 



(7.31) 



(7.32) 



(7.33) 



P = V /T73(37r 2 n ) 1/3 [as in Eqs. (|3~29l and (j3~3U)) ]. and 



e°(a,w) = e s (z = 0;q,u). 



(7.34) 



An inspection of Eqs. (|7.30j) - (j7.34|) shows that 
the hydrodynamic surface-response function g(q,uj) 
and, therefore, the hydrodynamic screened interaction 
W(z, z';q,u>) become singular at the hydrodynamic 
surface-plasmon condition dictated by Eq. (|3.30p . We 
also note that the second moment of the imaginary part 
of the hydrodynamic surface-response function g(q,u>) is 
found to be 



dui cj Img(q 7 ui) = 2tt 2 n, 



(7.35) 



where n represents the electron density: n = lo 2 /At:. 

Finally, we note that in the long-wavelength (q — > 
0) limit the hydrodynamic screened interaction of 
Eqs. (|7.30|) - (|7.34|) reduces to the classical screened inter- 
action of Eqs. (|7.8[) and (|7.9p with the dielectric functions 
ei and £2 being replaced by the Drude dielectric function 
[Eq. (|2.15[) ] and unity, respectively. The same result is 
also obtained by simply assuming that the electron gas 
is nondispersive, i.e., by taking the hydrodynamic speed 
P equal to zero. 

b. Quantum hydrodynamic approach. Within a 
quantized hydrodynamic model of a many-electron sys- 
tem, one first linearizes the hydrodynamic Hamiltonian 
with respect to the induced electron density, and then 
quantizes this Hamiltonian on the basis of the normal 
modes of oscillation (bulk and surface plasmons) corre- 
sponding to Eqs. $72% and (pT30|) . One finds: 



H = H r 



H B 



(7.36) 



where Hq represents the Thomas-Fermi grou nd state of 
the static unperturbed electron system |l8l| , and Hq 
and Hq are free bulk and surface plasmon Hamiltonians, 
respectively: 

H ° = h E t 1 / 2 + "q] «q (*)«Q (*) (7-37) 



and 



(7.38) 



Here, £1 and A represent the normalization volume and 
the normalization area of the surface, respectively, aQ(t) 
and bq(t) are Bose-Einstein operators that annihilate 
bulk and surface plasmons with wave vectors Q = (q, q z ) 
and q, respectively, and loq and ui^ represent the disper- 
sion of bulk and surface plasmons: 



P 2 Q 2 



and 



+ p 2 q 2 + Pq 



P 2 q 1 



(7.39) 



(7.40) 



Hence, within this approach one can distinguish 
the separate contributions to the imaginary part of 
the hydrodynamic surface-response function g(q,u>) of 
Eq. (|7.32p coming from the excit ation of either bulk or 
surface plasmons. One finds |254| : 



where 



lmg(q,uj) = lmg B (q,uj) +liag s (q,uj), 
1 f°° 

Im.9 B (<?,^) = 2 1 / dq z 5{uj - ujq) 

H/uq) £ 

X qi + q 2 (q 2 +^ 2 /P 2 )+^ p /(^) 



(7.41) 



and 



Img s (q,u>) 



7, 



2 q + 2 lq u)- q 
with Q = yjq 2 + q 2 and 

-pq + ^2^7^) . 



1 

2d 



(7.42) 



(7.43) 



(7.44) 
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For the second moments of Img B (q, u>) and lmg s (q,ui), 
one finds: 



d 7T q n 

dujujlmq (q.Lj) = - — w. 

o 4g + 2 lq p 



and 



y dwwlmg (q,u) = - 



7T 2 7g 2 



9 + 2 7g P ' 



(7.45) 



(7.46) 



which add up to the second moment of Eq. (|7.35|) . 

In the limit as q — > the bulk contribution to the 
so-called energy-loss function Img(q,uj) [see Eqs. (|7.41[) - 
(|7.44[) ] vanishes, and the imaginary part of both 
Eqs. (|7.32p and ()7.43[) yields the classical result: 



Img(q,uj) -> - u s 6(u - u s ), 



(7.47) 



which can also be obtained from Eq. (|7.9[) with ei re- 
placed by the Drude dielectric function of Eq. (|2.15p and 
£2 set equal to unity. Equation (|7.47[) shows that in the 
classical (long-wavelength) limit the energy loss is dom- 
inated by the excitation of surface plasmons of energy 



/\/2, as predicted by Ritchie. 



2. Specular-reflection model (SRM) 

An alternative scheme to incorporate nonlocal effects, 
which has the virtue of expressing the screened inter- 
action W(z, z';q,uj) in terms of the dielectric function 
e(Q,iv) of a homogeneous electron gas representing the 
bulk material, is the so-called spec ular- reflection model 
reported indep ende ntly by Wagner [l85l | and by Ritchie 
and Marusak 186]. In this model, the medium is de- 
scribed by an electron gas in which all electrons are con- 
sidered to be specularly reflected at the surface, thereby 
the electron density vanishing outside. 

For a semi-infinite metal in vacuum, the unperturbed 
electron density uq(z) is taken to be of the form of 
Eq. (|3.25|) , and the SRM yields a screened interaction of 
the form of Eq. (|7.30[) but with the quantities c s (z; q,tu) 
and g{q 1 uj) being replaced by the more general expres- 
sions: 



ts{z;q,uj) 



™%j«**e-\Q,u) (7.48) 



and 



l-e°(g,a;) 
l + c2(?.w)' 



(7.49) 



with e®(q,ui) defined as in Eq. (|7.34[) . and Q = \J q 2 + q 2 . 

The inverse dielectric function e~ 1 (Q,uj) entering 
Eq. (|7.48[) represents the 3D Fourier transform of the 
inverse dielectric function e _1 (r, r',oj) of a homogeneous 
electron gas. From Eq. (|6.3p . one finds: 

e- 1 (Q,u) = l + v QX (Q,u), (7.50) 



where x(Qi w ) represents the 3D Fourier transform of the 
density-response function x(r, r'; u>). 

In the framework of TDDFT, one uses Eq. ([5~TTjl to 
find 



x(Q^) = x°(Q^) + x (Q,uj) 
x {vq + Uc(n;Q,u)} x(Q,u), 



(7.51) 



with x°(Qi UJ ) an d fxc{n]Q,oj) being the 3D Fourier 
transforms of the nonintcracting density-response func- 
tion and the XC kernel of Eqs. (|5.7[) and (|5.15[) . respec- 
tively. For a homogeneous electron gas, the eigenfunc- 
tions tpi( r ) entering Eq. (|5.7[) arc all plane waves; thus, 
the integrations can be carried out analytic ally to yield 
the well-known Lindhard function x°(Qj uj ) [255| . If one 
sets the XC kernel f X c{n]Q,oj) equal to zero, introduc- 
tion of Eq. (|73T|) into Eq. (177511)) yields the RPA dielectric 
function 



(Q,u) = l-v Q x°(<9,w), 



(7.52) 



which is easy to evaluate. 

The RPA dielectric function e RPA (Q, ui) of a homo- 
geneous electron gas can be further approximated in the 
framework of the hydrodynamic scheme described in Sec- 
tion UTTAl One finds, 



^y dro {Q,uj) = i + 



ii 2 



u>(lu + i-q) ' 



(7.53) 



which in the classical (long-wavelength) limit yields the 
local Drude dielectric function of Eq. (|2.15p . Introduction 
of Eq. ([7753)1 into Eq. ((7748)) yields the hydrodynamic 
screened interaction of Eqs. (|7.30|) - ()7.34|) . 

We know from Eq. (|7.1|) that collective excitations 
are dictated by singularities in the screened interaction 
or, equivalently, maxima in the imaginary part of this 
quantity. For z and z' coordinates well inside the solid 
(z, z' — > -co), one finds 



W in (z,z';q,uj) 



Z7T 



(7.54) 

which in the case of the Drude dielectric function e(Q, lu) 
of Eq. (|2.15|) and for positive frequencies (lu > 0) yields 

lmW ln {z,z , ;q,uj) -> -— uj p 5{uj - u p ) c~ q \ z ~ z L (7.55) 

For z and z' coordinates both outside the solid (z, z' > 
0), one finds 

W out (z,z';q,uj) = — \e- q \ z - z '\ - g(q,oj) C" 1 ^^") 
q L 

(7.56) 

0) limit and for positive fre- 



which in the classical (q 
quencies (w > 0) yields 

lmW out (z,z';q,uj) -> 



■— lu s 5(uj-uj s )c- q(z+z '\ 

(7.57) 




CO (a.u.) co (a.u.) 



FIG. 8: The solid line represents the energy-loss function, 
Im[— W(z, z'; q, ui)], versus ui, as obtained at z — z' , q — 
0Aq F , and r s = 2.07 from Eq. (|7.54|l by using the full RPA 
dielectric function e RPA (Q, w). The thick dashed and dot- 
ted lines represent separate contributions from the excitation 
of bulk collective modes and e-h pairs occurring at energies 
u>Q =q < ui < u)q = q c and u> < qqF + q 2 /2, respectively. The 
vertical dotted line represents the energy u> p = 15.8 eV at 
which collective oscillations would occur in a Drude metal. 



Figures [5] and [5] show the energy-loss function 
Im[— W(z, z'; q, ui)] that we have obtained at z = z', 
q = 0Aq F , and r s = 2.07 from Eqs. ([7341) and lff36|) . 
respectively, by using the full RPA dielectric function 
e RPA (Q, u>). For z coordinates well inside the solid 
(Fig. [8} , instead of the single classical collective excita- 
tion at ui p (dotted vertical line) predicted by Eq. (|7.55p 
the RPA energy-loss spectrum (solid line) is composed 
of (i) a continuum of bulk collective excitation s (da shed 
line) occurring at energies u)q =q < uj < u>q = q c [256| |. and 
(ii) the excitation of electron-hole (e-h) pairs represented 
by a thick dotted line. 

For z coordinates that are outside the surface, it had 
been generally believed that only surface plasmons and 
c-h pairs ca n be excite d. H owever, it was shown explic- 
itly in Refs. [254| and |257j that the continuum of bulk- 
plasmon excitations dominating the energy-loss spectrum 
inside the solid [see Fig. [8] is still present for z c oordi 



nates outside, as shown in Fig. [S]for z = z' = Af |258 |. 
This continuum, which covers the excitation spectrum at 
energies u> > Wg =5 and is well separated from the lower- 
energy spectrum arising from the excitation of surface 
plasmons and e-h pairs, is accurately described by using 
the quantal hydrodynamic surface energy-loss function 
of Eq. (|7.42[) . which has been represented in Fig. [5] by a 
thick dotted line. 

Nonetheless, the main contribution to the energy-loss 
spectrum outside the solid comes from the excitation of 
surface plasmons, which are damped by the presence of 
c-h pairs. These c-h pair excitations are not present in 
the classical and hydrodynamic schemes described above, 
which predict the existence of long-lived surface plas- 



FIG. 9: The solid line represents the energy-loss function, 
Im[— W(z, z'\ q, ui)], versus u), as obtained at z = z' — Af, q = 
0.4q F , and r s = 2.07 from Eq. (fT56|l by using the full RPA di- 
electric function e RPA (Q,cu). The thick dotted line, which is 
nearly indistinguishable from the solid line covering the same 
part of the spectrum, represents the hydrodynamic prediction 
from Eq. (|7.42|l . The thin dotted vertical lines represent the 
energies ui s = 11.2 eV and = 15.4 eV of Eq. (|7^u] at which 
long-lived surface plasmons would occur in a semi-infinite 
metal described by a Drude and a hydrodynamic model, re- 
spectively. Introduction of the full RPA dielectric function 
e RPA {Q,u)) into Eq. (|7,49[) yields a maximum of the surface- 
loss function Im<?(g, ui) (and, therefore, lm[—W(z,z';q,uj)]) 
at the surface-plasmon energy u>q = 16.0 eV, which is slightly 
larger than its hydrodynamic counterpart. 



mons at the energies represented in Fig. [9] by thin dotted 
vertical lines: ui s = ui p /\/2 [sec Eq. (|7.57| ] and of 
Eq. fTTlCJj) [see Eqs. (J733J -([131])], respectively. 



3. Self-consistent scheme 

For an accurate quantal description of the electronic 
excitations that can occur in a semi-infinite metal, 
we need to consider the true density-response function 
x(r, r';ui) entering Eq. (|7.4j) . which is known to fulfill the 
/-sum rule of Eq. (|5.3[) . 

a. Jellium surface. In the case of a free-electron gas 
bounded by a semi-infinite positive background of density 



n + {z) 



n, z < 0, 
0, z > 0, 



(7.58) 



translationally invariance in the plane of the surface al- 
lows to define the 2D Fourier transform W(z, z'\ q, ui), 
which according to Eq. (|7.4[) can be obtained as follows 



W(z, z'; q, ui) = v(z, z'; q) + J dz\ J dz 2 v(z, Z\\q) 
xx{zi,Z2;q,uj)v(z 2 ,z';q), (7.59) 
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where v(z, z'; q) is the 2D Fourier transform of the bare 
Coulomb interaction v(r, r'): 



v(z,z';q) = — e"" 12 

q 



(7.60) 



and x{ z i z 'i <?i w ) denotes the 2D Fourier transform of the 
interacting density-response function %(r, r';u>). In the 
framework of TDDFT, one uses Eq. (f!HT|) to find: 

X( z , z ';<l,v) = x°{z,z';q 7 uj) + J dz 1 J dz 2 X°( z , Z H q, 

x {v(z 1 ,z 2 ;q) + f xc [no]( z i, z 2\ Q, ")} x( z 2, z '\ q, <*>), ( 7 - 61 ) 

where x°( z i z '; 9: w ) an d /xc[^o] i z , z'\ q, w) denote the 2D 
Fourier transforms of the noninteracting density-response 
function x°( r 5 r ';w) and the XC kernel f xc [n ](r,r';uj), 
respectively. Using Eq. (|5.7|) . and noting that the single- 
particle orbitals tpi(r) now take the form 



x n (z)5(z- z'). (7.70) 

Within this scheme, the simplest possible approxi- 
mation is to neglect XC effects altogether and set the 
XC potential v xc [no](z) and kernel f xc [no](z,z';q,Lj) 
equal to zero. In this case, the onc-dimcnsional single- 
particle wave functions V J i( z ) an d energies E\ are the 
self-consistent eigenfunctions and eigenvalues of a one- 
dimensional Hartree Hamiltonian. The calculation of 
the density-response function is further simplified if the 
Hartree potential vh [tlq] (z) of Eq. (|7.67|) is replaced by 



,ik-r, t 



^i( z ), 



(7.62) 



one finds: 



x°( z , z ';q,u) = 7 E M z )^( z )^( z 'Wi( z ') 



E 



A 

Jk,j — /k+qj 



-Ek.i - ^k+q,j + W + 27/ ' 



where 



e 4 



(7.63) 



(7.64) 



the single-particle orbitals i/ji( z ) an( i energies £j now 
being the solutions of the one-dimensional Kohn-Sham 
equation 



-~+VK8[no]{z) 



tpi(z) = e t ipi{z), (7.65) 



with 



VKs[no]{z) = v H [n Q ](z) + v xc [n ](z), (7.66) 
v H (z) = -2ir / dz' \z - z'\ [no(z') - n + (z')] , (7.67) 



Vxc[no]{z) = 



6E X 



dn(z) 



(7.68) 



and 



no 



( z ) = - V(£f - e t ) $( z ) °(cf - si). (7.69) 

7T * " 



From Eq. (|5.3j) . the imaginary part of the density- 
response function z'; q, u>) is easily found to fulfill 
the following sum rule: 



duj u> lmx(z, z\ q, u>) 



dzdz' 



vibm( z ) 



z < z 0, 



z > ZQ, 



(7.71) 



where the value z = (3/16)A_f is chosen so as to ensure 
charge neutra lity. This is the so-called inifinite-barrier 
model (IBM) 259], in which the single-particle orbitals 
ipi(z) are simply sines. If one further neglects interference 
between incident and scattered electrons at the s urface, 
this model yields the classical IBM (CIBM) [260| which 
can be shown to be equivalent to the SRM described in 
Section EHBU 

Alternatively, and with the aim of incorporating band- 
structure effects (such as the presence of energy gaps and 
surface states) approximately, the self-consistent jellium- 
like Kohn-Sham potential of Eq. (|7. 66[) can be replaced 
by a physically motivated model potential vmp{ z )- Ex- 
amples are the pa rame terized model potential reported 
by Chulkov el al. |26l| , which was successful in the de- 
scription of the lifetimes of image and Shockley states in 
a variety of metal surfaces |262tl263l 1263 , [265TI266I. [267l 
|268| , and the g-dependent model potential that has been 
reported recently to investigate the momentum-res olved 
lifetimes of Shockley states at the Cu(lll) surface [269j |. 

At this point, we note that for z and z' coordinates 
that are far from the surface into the vacuum, where the 
electron density vanishes, Eq. (|7.59|) takes the form of 
Eq. ([7T56]) (which within the SRM is true for all z, z' > 0), 
i.e.: 

W(z, z'; q, u) = v(z, z'; q) - — e~^ z+z > g(q, cj), (7.72) 

q 

but with the surface-response fu nctio n g(q,u>) now being 
given by the general expression |270j : 

g(q,u>) = -— J d Zl J dz 2 e^+^) x (z 1 ,z 2 ;g,o.), 

(7.73) 

which according to Eq. (|5.ip can be expressed as follows 



g(q,oj) = / dzc qz Sn(z;q 7 uj) 7 



(7.74) 



with Sn(z; q, oj) being the electron density induced by an 
external potential of the form 



<t> ext (z;q,uj) 



2tt 



(7.75) 
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In the framework of TDDFT, the induced electron den- 
sity is obtained as in the RPA [see Eq. (|5.5[) ] , but with the 
XC kernel f xc [n }(r,r';q,uj) added to the bare Coulomb 
interaction u(r,r'). Hence, after Fourier transforming 
one writes 

6n(z;q,u) = J dz' x°(z, z'; q, uj) jV^O'; q, uS) + J dz" 
x [v(z',z";q) + f xc [n ](z' , z";q,uj)]Sn(z";q,uj)} . (7.76) 



Using Eqs. <|7.70[) and (|7.73p . the surface loss function 
lmg(q, uj) is easily found to fulfill the following sum rule: 

J dujculmg(q,uj) = 2tt 2 q J dz e 2qz n a (z), (7.77) 

which for a step- like electron density hq(z) of the form 
of Eq. (|3.25p reduces to Eq. (|7.35[) . as expected. 

b. Periodic surface. For a periodic surface, single- 
particle wave functions are of the form 



-0k,n;i(l") = -0k,n(r||) ^(z) 

where i/jk t n(r\\) are Bloch states: 



(7.78) 



Vk,n(ri|) = ^=e k - r " Wk , n (r||), (7.79) 
V A 

with r II and k being 2D vectors in the plane of the sur- 
face. Hence, one may introduce the following Fourier 
expansion of the screened interaction: 



SBZ 



q s,g' 
x Wg, g /(z,z';q,w), 



(7.80) 



where q is a 2D wave vector in the surface Brillouin zone 
(SBZ), and g and g' denote 2D reciprocal-lattice vec- 
tors. According to Eq. (|7.4p , the 2D Fourier coefficients 
Wg ig ' (z, z'; q, uj) are given by the following expression: 

W s , g >(z, z';q, uj) = v g (z, z';q)8 g!g : + J dzi J dz 2 
xu g (z, zi; q) x g , g '(^i, z 2 ; q, to) v s >(z 2 ,z'; q), (7.81) 

where v s (z,z';q) denote the 2D Fourier coefficients of 
the bare Coulomb interaction v(r, r'): 



v g (z,z';q) = 



2vr 



q + g 



-|q+g| \z—z 



(7.82) 



and x gjg /(z, z'; q, w) are the Fourier coefficients of the 
interacting density-response function x( r , r '; w )- In the 
framework of TDDFT, one uses Eq. (|5.1ip to find: 

Xs,g'(z,z';<i,u) =x gtS >(z,z'\q,u) + J dz x J dz 2 
x X s , s '(z,Zi;q,uj) x [v gl (zx, z 2 ; q) S gug2 



+/ g J r,g2["-o](^i,^2;q,w)] Xg2,g'(^2,z';q,w), (7.83) 

where x°g, s > ( z i z '\ 1, u) and f™ g ,{n ](z, z'\ q, uj) denote 
the Fourier coefficients of the nonintcracting density- 
response function x ( r ? r '; w ) and the XC kernel 
fxc[no}(^, r'; uj), respectively. Using Eq. (|5.7|) . one finds: 



x" ,(z,z';q,u) 



L g,g 

SBZ 



-y 



^{z)r 3 {z)^{z'wd z ') 



EE 



/k , n ; i ,/k + q , n ' ; j 



k n,n 



/ 6 K 

-i(q+g)-ry 



^k+q,, 



h(tu + irj) 



x(^ k ,„|e- I ^- r ii|Vk +q , Il ' 

X^k+q.n'l^^'^IV'k,™) 



(7.84) 



the single-particle orbitals i/Ju. >n -i(r) = ^k,n(i"||) ipi(z) and 
energies £k,n;i being the eigenfunctions and eigenvalues of 
a 3D Kohn-Sham Hamiltonian with an effective potential 
that is periodic in the plane of the surface. 

As in the case of the jellium surface, we can focus on 
the special situation where both z and z' coordinates are 
located far from the surface into the vacuum. Eq. (|7.4j) 
shows that under such conditions the Fourier coefficients 
Wg !g ' (z, z'; q, uj) take the following form: 



W StS '(z,z';ci,uj) = v s (z,z';q)6 SyS > - 



2irq 



-g| |q + g I 

x.9 g , g ,(q, W )e-l^ e -|q+gV ; (7 . 85) 
where 

g s ,g>(<l,u) = -y J dzi J dz 2 c q(zi+Z2) Xg, g >(zuz 2 ;(i,uj). 

(7.86) 

In particular, 

5g=o, g =o(q, w) = J dz e qz 5ng= (z; q, uj), (7.87) 

with Sn s (z; q, uj) being the Fourier coefficients of the elec- 
tron density induced by an external potential of the form 



bl xt (z;q,uj) = -—e" z S s , . 
1 



(7. 



Finally, one finds from Eq. (|5.3[) that the Fourier coef- 
ficients Xg.g'( z i z ' \ w ) fulfill the following sum rule: 



/oo 

duj uj Imx g , s '(z,z';q,uj) = —n 
-oo 

x n°g_g,(z;q L )5(z- z'), 



dzdz' 



(7.89) 



where the coefficients n g (z;q) denote the 2D Fourier 
components of the ground-state electron density ?io( r )- 
Furthermore, combining Eqs. (|7.86p and (|7.89p . one 
writes 

J du> ujlmgg : g'(q,uj) = 2ir 2 q J dze 2q * n g _ g <(z; q), 

(7.90) 
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which is a generalization of the sum rule of Eq. (|7.35|) to 
the more general case of a real solid in which the crystal 
structure parallel to the surface is taken into account. 



VIII. SURFACE-RESPONSE FUNCTION 




The central quantity that is involved in a description 
of surface collective excitations is the surface-response 
function introduced in the preceding section: g(q,uj) for 
a jcllium surface [see Eq. Q7.73|) ] and ffg g '(q, u>) for a 
periodic surface [see Eq. (|7.86[) ]. 



A. Generation-rate of electronic excitations 

In the framework of TDDFT, an interacting many- 
electron system exposed to a frequency-dependent exter- 
nal potential (f> ext (r, u>) is replaced by a fictitious system 
of noninteracting electrons exposed to an effective self- 
consistent potential </> sc (r, uj) of the form 



(8.1) 



where 



S<f>{t,u) = / dr' [v{r,r') + f xc [n ](r,r';uj)} Sn{r';cj). 



(8.2) 

Hence, the rate at which a frequency-dependent exter- 
nal potential <fi ext (r, u>) generates electronic excitations in 
the many-electron system can be obtained, within lowest- 
order perturbation theory, as follows 

W ( W ) = 2^/ ! (l-/ 1 )|(^(r)|r(r,u;)^(r))| 2 S(e t -e f ) 7 

i,j 

(8.3) 

with ipi(r) and £j being the eigenfunctions and eigen- 
values of a 3D Kohn-Sham Hamiltonian. In terms of 
the i nduc ed electron density Sn(r, uS) of Eq. (|5.ip . one 
finds [163 



w(uj) = -21m / dr<f)* ext (r,uj)5n(r,uj), (8.4) 

which in the case of a periodic surface takes the following 
form: 



SBZ 



(8.5) 



g q 



where w g (q, uj) denotes the rate at which the external 
potential generates electronic excitations of frequency w 
and parallel wave vector q + g: 



w g (q,w) = -— Im / dz</) g xt (z,q)5n s (z,q), (8.6) 

</> g x * (z, q) and Sn g (z,q) being the Fourier coefficients of 
the external potential and the induced electron density, 
respectively. 




FIG. 10: A schematic drawing of the scattering geometry in 
angle-resolved inelastic electron scattering experiments. On 
exciting a surface mode of frequency ui(q), the energy of de- 
tected electrons becomes e/ = gj — ui(q), with the momentum 
q being determined by Eq. (|8.9|l . 



In particular, if the external potential is of the form of 
Eq. (|7.88p . the rate w; g (q, w) can be expressed in terms 
of the surface-response function <? g , g '(q, uj) of Eq. (|7.86[) . 
as follows 

4-7T 

w g (q, uj) = — Im ffgjg (q, u) <5 gi0 , (8.7) 
which for a jcllium surface reduces to 

47T 

w(g, w) = — - Img(q, uj), (8.8) 
qA 

with g(q,u>) being the surface-response function of 
Eq. (|7.73p . We note that although only the coefficient 
X g= o. g ' = o(2, z'; q, w) enters into the evaluation of the 
more realistic Eq. (|8.7|) . the full x g g ' (z, z'\ q, uj) matrix 
is implicitly included through Eq. (|7.83[) . 



B. Inelastic electron scattering 

The most commonly used experimental arrangement 
for the detection of surface collective excitations by the 
fields of moving charged particles is b ased on angle- 
resolved inelastic electron scattering 27l|. Fig. [TOl shows 
a schematic drawing of the scattering geometry. A 
monochromatic beam of electrons of energy Si, incident 
on a flat surface at an angle 9i, is back scattered and 
detected by an angle-resolved energy analyzer positioned 
at an angle 8f and energy Sf. Inelastic events can occur, 
either before or after the elastic event, on exciting a sur- 
face mode of frequency uj(q) = e% — £/ . The energy and 
lifetime of this mode are determined by the correspond- 
ing energy-loss peak in the spectra, and the momentum 
q parallel to the surface is obtained from the measured 
angles 6i and 8f, as follows: 



= V2[ 



e f smtf/j 



(8.9) 



The inelastic scattering cross section corresponding to 
a process in which an electronic excitation of energy w 



2G 



and parallel wave vector q is created at a semi-infinite 
solid surface can be found to be proportional to the rate 
given by Eq. (|8.7p [or Eq. (|8.8jl if the medium is rep- 
resented by a jellium surface] and is, therefore, propor- 
tional to the imaginary part of the surface-response func- 
tion, i.e., the so-called surface-loss function. Hence, apart 
from kinematic factor s (which can indeed vary with en- 
ergy and momentum [272l . |273| | ) , the inelastic scattering 
cross section is dictated by the rate u> g (q, w) [or w(q, uj) if 
the medium is represented by a jellium surface] at which 
an external potential of the form of Eq. (|7.88[) generates 
electronic excitations of frequency uj and parallel wave 
vector q. 

In the following sections, we focus on the behavior of 
the energy-loss function Imgg = o,g=o(q, uj) [or lmg(q, uj)} 
and its maxima, which account for the presence and 
momentum-dispersion of surface collective excitations. 



(3 representing the speed of propagat ion o f hydrodynamic 
disturbances in the electron system [l82j . 

In the SRM (with the bulk dielectric function being 
described within the RPA), surface plasmons, which oc- 
cur at a momentum-dependent energy slightly different 
from its hydrodynamic counterpart, are damped b y th e 
presence of e-h pair excitations, as shown in Fig. l9l [274| . 

The one-step hydrodynamic Eqs. (|8.12|) and (|8. 13|) and 
a numerical evaluation of the imaginary part of the SRM 
surface-response function of Eq. (|7.49[) (see Fig. \§§ both 
yield a positive surface-plasmon energy dispersion at all 
wave vectors. Nonetheless, Bennett used a hydrodynamic 
model with a continuum decrease of the electron den- 
sity at the metal surface, and found that a continuous 
electron-density variation yields a monopole surface plas- 
mon with a negative dispersion at low wave vectors [85[ . 



C. Surface plasmons: jellium surface 

1. Simple models 

In the simplest possible model of a jellium surface in 
vacuum, in which a semi-infinite medium with local di- 
electric function e(uj) at z < is terminated at z = 0, 
the surface-response function g(q,uj) is obtained from 
Eq. (JLU) with e 2 = 1, or, equivalently, from Eq. ([7749]) 
with q = 0, i.e.. 



9(q, u) 



e(oj) + 1' 



(8.10) 



which for a Drude dielectric function [see Eq. (|2.15p ] leads 
to the surface-loss function 



lmg(q,uj) = -uj s S(uj 



UJ S ) 



(8.11) 



peaked at the surface-plasmon energy uj s = uj v j\[2. 

The classical energy-loss function of Eq. (|8.10p repre- 
sents indeed the true long-wavelength (q — > 0) limit of 
the actual self-consistent surface-loss function of a jel- 
lium surface. Nevertheless, the classical picture leading 
to Eq. (|8.10p ignores both the nonlocality of the electronic 
response of the system and the microscopic spatial distri- 
bution of the electron density near the surface. Nonlocal 
effects can be incorporated within the hydrodynamic and 
specular-reflection models described in Sections IVIIB II 
and lVHB21 

Within a one-step hydrodynamic approach, the 
surface-loss function is also dominated by a delta function 
[see Eq. (|7.43p ] but peaked at the momentum-dependent 
surface-plasmon energy of Eq. (|7.40p [see also Eq. (|3.30p ]: 



uj 2 = - [u> 2 p + (3 2 q 2 + Pq^2uj 2 + f3 2 q 2 
which at long wavelengths yields 

uj = w p /V2 + /3q/2, 



.12) 



(8.13) 



2. Self- consistent calculations: long wavelengths 

Within a self-consistent long-wavelength description 
of the jcllium-surface electronic response, Feibelman 
showed that up to first order in an expansion in powers of 
the magnitude q of the wave vector, the surfa ce-response 
function of Eq. (|7.73p can be written as [l66j 



[eM-1] [l + grf±M] 
e(w) + 1 - [e(w) - 1] qd±(Lo)'' 



(8.14) 



where e(uj) represents the long- wavelen gth limit of the 
dielectric function of the bulk material [275j ] and dj_ (uj) 
denotes the centroid of the induced electron d ensity 
[see Eq. (|3.13p ] with respect to the jellium edge [2761 ]. 
Eq. (|8.14p shows that at long wavelengths the poles of 
the surface-response function g(q,uj) are determined by 
the non-retarded surface-plasmon condition of Eq. (|3.12p 
with 62 = 1, which for a semi- infinite free-electron metal 
in vacuum yields the surface-plasmon dispersion relation 
of Eq. (j3~T5)l . i.e.: 



with 



uj = uj s (1 + aq) , 



-Re[d ± (w a )]/2. 



(8.15) 



(8.16) 



Equations (|8.15p - (|8.16p show that the long- wavelength 
surface-plasmon energy dispersion is dictated by the po- 
sition of the centroid of the induced electron density with 
respect to the jellium edge. This can be understood by 
noting that the potential associated with the surface- 
plasmon charge attenuates on either side of Re [dj_(u; s )] 
with the attenuation constant q. If the fluctuating charge 
lies inside the jellium edge (Re[dj_(w s )] < 0), as q in- 
creases (thus the surface-plasmon potential attenuating 
faster) more of the field overlaps the metal giving rise to 
more interchange of energy between the electric field and 
the metal and resulting in a positive energy dispersion. 
However, if the fluctuating charge lies outside the jellium 
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FIG. 11: Long- wavelength surface-plasmon dispersion coef- 
ficient a entering Eq. (|8.15[) . versus the electron-density pa- 
rameter r s , as obtained from Eq. (|8.17[1 (dotted line) and 
from SRM (thick dashed line), IBM (thin dashed line), and 
self-consistent RPA and ALDA (thin and thick solid lines, re- 
spectively) calculations of the centroid of the induced elec- 
tron density at u = u> s . The solid circles represent the 
angle-resolved low-energy inelastic electron scattering mea- 
surements repo rted in Refs. for Na and K, i n R ef. [lOl| 
for Cs, in Ref. [Tol for Li and Mg, and in Ref. [H3| for Al. 
The IBM and self-consistent (RPA and A LDA) calculations 
have been taken from Refs. [89| and [l06f | . respectively. 



edge (Re[dj_(w s )] > 0), as q increases less of the metal 
is subject to the plasmon's electric field (i.e., there is a 
decreasing overlap of the fluctuating potential and the 
unperturbed electron density) which results in less inter- 
change of energy between the electric field and the metal 
and a negative dispersion coefficient [277l . l278j | . 

Quantitative RPA calculations of the surface-plasmon 
linear-dispersion coefficient a entering Eq. (|8.15p were 
carried out by several authors by using the specular- 
reflection and infinite-barrier models of the surface (83l . 
HH HI], a step potential [9(J HH, an d the more re- 
alist i c La ng-Kohn self-consistent surface potential (93l . 
llOrl 12791 ]. Both Feibelman's RPA self-consistent cal- 
culations [93[ an d th e ALDA calculations carried out 
later by Liebsch [l06| | and by Kempa and Schaich [279j | 
demonstrated that in the range of typical bulk densities 
(r s = 2 — 6) the centroid d± of the induced electron 
density at oj s lies outside the jellium edge, which leads 
to a negative long-wavelength dispersion of the surface 
plasmon. These calculations, which corroborated Ben- 
nett's prediction (8f|, also demonstrated that the long- 
wavelength surface-plasmon dispersion is markedly sen- 
sitive to the shape of the barrier and to the presence of 
short-range XC effects. 

Existing calculations of the long- wavelength dispersion 
coefficient a entering Eq. (|3.15[) are shown in Fig. II II and 
summarized in Table HI At one extreme, the dotted line 
of Fig. [11] gives the single-step hydrodynamic coefficient 



TABLE I: Long-wavelength surface-plasmon dispersion co- 
efficient a in A for various simple metal surfaces, as ob- 
tained from Eq. (pUT)) (HD) and from SRM, IBM, and self- 
consistent RPA and ALDA calculations of the centroid of the 
induced electron density at uj = uj s [see Eq. (|8.16[) ] , and from 
angle- resolved low-energy inelastic electron scattering mea- 
surements repo rted in Refs. [9^] for Na and K, i n R ef. [I0lf | 
for Cs, in Ref. [Tol for Li and Mg, and in Ref. [H| for Al. 
As in Fig. mi the IBM and self-consistent (RPA and ALDA) 
calculations have been taken from Refs. [89l [ and |l06l [. respec- 
tively. Also shown in this table are the measured values of the 
surface-plasmon energy ui a at q = 0, which are all slightly be- 
low the jellium prediction: w p /\/2 = y%3/2rfe 2 /ao due to 
band-structure effects. 





r a 




HD 


SRM 


IBM 


RPA 


ALDA 


Exp. 


Al 


2.07 


10.86 


0.46 


0.50 


0.36 


-0.21 


-0.32 


-0.32 


Mg 


2.66 


7.38 


0.52 


0.57 




-0.30 


-0.50 


-0.41 


Li 


3.25 


4.28 


0.58 


0.63 


0.33 


-0.40 


-0.70 


-0.24 


Na 


3.93 


3.99 


0.64 


0.70 




-0.45 


-0.85 


-0.39 


K 


4.86 


2.74 


0.71 


0.78 




-0.35 


-1.10 


-0.39 


Cs 


5.62 


1.99 


0.76 


0.84 


0.33 


-0.26 


-1.14 


-0.44 



a HD obtained from Eq. (|3.3ip . i.e. 



HD 



JL 

2lo, 



.17) 



with (3 = y / 375(37r 2 no) 1 / 3 ; within this model, the long- 
wavelength surface-plasmon dispersion is always positive. 
As in the single-step hydrodynamic approach, the SRM 
equilibrium density profile is of the form of Eq. (|3.25j) . 
and in the IBM the electron density still varies too 
rapidly; as a result, the corresponding SRM and IBM a 
coefficients (represented by thick and thin dashed lines, 
respectively) are both positive. At the other extreme are 
the more realistic self-consistent RPA and ALDA calcu- 
lations, represented by thin and thick solid lines, respec- 
tively, which yield a linear-dispersion coefficient a that 
is negative in the whole metallic density range. 

Conclusive experimental confirmation that the original 
Bennett's prediction [H[ was correct came with a series 
of measurements based on ang l e-resolved low- e nerg y in- 
elastic electron scattering [if [iM [lOll [To! [28(|. It 
has been demonstrated that the surface-plasmon energy 
of simple metals disperses downward in energy at small 
momentum q parallel to the surface, the dispersion coef- 
ficients (represented in Fig. [TT] by solid circles) being in 
reasonable agreement with self-consistent jellium calcu- 
lations, as shown in Fig. [TT1 and Table [T] 



3. Self- consistent calculations: arbitrary wavelengths 

At arbitrary wavelengths, surface-plasmon energies 
can be derived from the maxima of the surface loss 
function Ivag{q,uj) and compared to the peak positions 
observed in experimental electron energy-loss spectra. 
Fig. [T2l shows the self-consistent calculations of lmg(q, u>) 
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FIG. 12: Self-consistent RPA calculations (thin solid lines) 
of the surface loss function \m.g(q,uj), as obtained from 
Eq. (|7.73p . versus the energy uj for a semi-infinite free-electron 
gas with the electron density equal to that of valence electrons 
in Al (r 3 — 2.07) and for various magnitudes of the 2D wave 
vector q. The thick solid line represents the SRM surface loss 
function for r s = 2.07 and q = 0.5 A -1 . 




that we have obtained in the RPA for a semi-infinite free- 
electron gas (jellium surface) with the electron density 
equal to that of valence electrons in Al (r s = 2.07). 
For 2D wave vectors of magnitude in the range q = 
0—0.5 A , all spectra are clearly dominated by a surface- 
plasmon excitation, which is seen to first shift to lower 
frequencies as q increases [as dictated by Eqs. (|8.15[> - 
(I8.16p ] and then, from about q = 0.15 A on, towards 
higher frequencies. 

For the numerical evaluation of the spectra shown in 
Fig. [T2] we have first computed the interacting density- 
response function x( z > z '\ 9> w ) of a sufficiently thic k jel - 
lium slab by following the method described in Ref . |28l| , 
and we have then derived the surface-loss function from 
Eq. (|7.73|) . Alternatively, the self-consistent ener gy-loss 
spectra reported in Ref. [lOll ] (see also Ref. [167| |) were 
obtained by first computing the nonintcracting density- 
response function x°( z > z 'l 9i w ) °f a semi-infinite electron 
system in terms of Green's functions, then performing a 
matrix inversion of Eq. (|7.76[) with the external poten- 
tial (j> ext (z]q,ui) of Eq. (|7.75|) . and finally deriving the 
surface- loss function from Eq. (|7.74p . As expected, both 
approaches yield the same results. 

Figure [T3] shows the calculated and measured disper- 
sion of surface plasmons in Al, Mg, Li, Na, K, and Cs. 
The jellium calculations presented here do not include 
effects due to band-structure effects; thus, all frequen- 
cies have been normalized to the measured value tu s of 
the g = surface-plasmon energy (see Table UJ) . This 
figure shows that the single-step HD and SRM surface- 
plasmon dispersions (thin and thick dashed lines, respec- 
tively) are always upward and nearly linear. However, 
self-consistent RPA and ALDA calculations (thin and 



FIG. 13: Surface-plasmon energy dispersion for the simple 
metals Al, Mg, Li, Na, K, and Cs, as obtained from Eq. (|8.12|) 
(thin dashed lines) and from SRM (thick dashed lines) and 
self-consistent RPA (thin solid lines) and ALDA (thick solid 
lines) calculations of the surface loss function lmg(q, uj), and 
from the peak positions observed in experimental electron 
energy- loss spectra (solid circles). The dotted lines represent 
the initial slope of the HD surface-plasmon energy dispersion, 
as obtained from Eq. (|3.31l In the case of Mg, the thick solid 
line with open circles repre sents the jellium TDDFT calcula- 
tions reported in Ref. [121| and obtained with the use of the 
nonlocal (momentum-dependent) static XC local-field factor 
of Eq. (|5.24[) . All frequencies have been normalized to the 
measured value uj s of the q = surface-plasmon energy. 



thick solid lines, respectively), which are based on a self- 
consistent treatment of the surface density profile, show 
that the surface-plasmon dispersion is initially downward 
(as also shown in Fig. [TT] and Table QJ, then flattens out, 
and rises thereafter, in agreement with experiment (solid 
circles). A comparison between self-consistent RPA and 
ALDA calculations show that although there is no quali- 
tative difference between them XC effects tend to reduce 
the surface-plasmon energy, thereby improving the agree- 
ment with the measured plasmon frequencies. This low- 
ering of the surface-plasmon energies shows that dynamic 
XC effects combine to lower the energy of the electron 
system, which is due to the weakening of the Coulomb 
e-e interaction by these effects. 

Recently, XC effects on the surface-plasmon dispersion 
of Mg and Al were introduced still in the framework of 
TDDFT (but beyond the ALDA) by using the nonlo- 
cal (momentu m-depend ent) static XC local-field factor 
of Eq. (|5.24p |l2ll . Il22| . At low wave vectors, this cal- 
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culation (thick solid line with open circles) nearly coin- 
cides with the ALDA calculation (thick solid line), as 
expected. At larger wave vectors, however, the nonlocal 
calculation begins to deviate from the ALDA, bringing 
the surface-plasmon dispersion to nearly perfect agree- 
ment with the data for all values of the 2D wave vector. 
Ab initio calculations of the surface-plasm on di spers ion o f 
real Al and Mg were also reported in Refs. |l2l| and |l22| . 
For plasmon frequencies that are normalized to the mea- 
sured value lo s at q = (as in Fig. 1 13|) . ab initio and 
jellium calculations arc found to be nearly indistinguish- 
able; however, only the ab initio calculations account for 
an overall lowering of the surface-plasmon dispersion that 
is due to core polarization. 

In the case of the alkali metals Li, Na, K, and Cs, there 
is a mismatch in the linear (low q) region of the surface- 
plasmon dispersion curve between jellium ALDA calcu- 
lations and the experiment. The theoretical challenges 
for the future are therefore to understand the impact of 
band-structure and many-body effects on the energy of 
surface plasmons in these materials, which will require to 
pursue complete first-principles calculations of the sur- 
face e lectr onic r espon se at the level of those reported in 
Refs. [mj and [l2H for Mg and Al, respectively. 

Thin films of jellium have been considered recently, in 
order to investigate the influence of the slab thickness 
on the excitati on sp ectra and the surface-plasmon en- 
ergy dispersion [282 j. Oscillatory structures were found, 
corresponding to electronic interband transitions, and it 
was concluded that in the case of a slab thickness larger 
than ~ 100 a.u. surface plasmons behave like the surface 
plasmon of a semi-infinite system. 



4- Surface-plasmon linewidth 

At jellium surfaces, the actual self-consistent sur- 
face response function Img(g, u>) reduces in the long- 
wavelength (q — > 0) limit to Eq. (|8.11[) . so that long- 
wavelength surface-plasmons are expected to be infinitely 
long-lived excitations. At finite wave vectors, however, 
surface plasmons are damped (even at a jel lium surface) 
by the presence of e-h pair excitations (283[. 

Figure [TJ] shows the results that we have obtained from 
SRM (dashed lines) and self-consistent (solid lines) RPA 
calculations of the full width at half maximum (FWHM), 
Aoj, of the surface loss function Img(g, u>). Also shown 
in this figure are the corres pond i ng li n ewid ths that have 
been reported in Refs. [99|, llOll . 1 1021 . l280j from experi- 
mental electron energy-loss spectra at different scatter- 
ing angles. This figure clearly shows that at real surfaces 
the surface-plasmon peak is considerably wider than pre- 
dicted by self-consistent RPA jellium calculations, espe- 
cially at low wave vectors. This additional broadening 
should be expected to be mainly caused by the presence 
of short-range many-body XC effects and interband tran- 
sitions, but also by scattering from defects and phonons. 

Many-body XC effects on the surface-plasmon 




FIG. 14: Surface-plasmon full width at half maximum 
(FWHM), Alo, for the simple metals Al, Mg, Li, Na, K, 
and Cs, as obtained from SRM (thick dashed lines) and self- 
consistent RPA (thin solid lines) calculations of the surface 
loss function Img(q,u>), and from the experimental electron 
energy-loss spectra at different scattering angle s (so lid cir- 
cles) re porte d in Refs. [§tj] for Na and K, in R ef. [l0ll | for Cs, 
in Ref. [To2 for Li and Mg, and in Ref. HH for Al. In the 
case of Mg, the thick solid line with open circl es re presents 
the jellium TDDFT calculations reported in Ref. [l 2 ll ] and ob- 
tained with the use of the nonlocal (momentum-dependent) 
static XC local-field factor of Eq. (|5.24|) . 



linew idth of Mg and Al were incorporated in Refs. [121| 
and fl22T | in the framework of TDDFT with the use 



of the nonlocal (momentum-dependent) static XC local- 
field factor of Eq. (pTMj) . These TDDFT calculations 
have been plotted in Fig. HUby solid lines with circles; a 
comparison of these results with the corresponding RPA 
calculations (solid lines without circles) shows that short- 
range XC correlation effects tend to increase the finite-q 
surface-plasmon linewidth, bringing the jellium calcula- 
tions into nice agreement with experiment at the largest 
values of q. Nevertheless, jellium calculations cannot pos- 
sibly account for the measured surface-plasmon linewidth 
at small <?, which deviates from zero even at q = 0. 

In the long-wavelength (q — ► 0) limit, surface plasmons 
are known to be dictated by bulk properties through 
a long- wavelength bulk dielectric function e{uj), as in 
Eq. (|8.10p . Hence, the experimental surface-plasmon 
widths Alj at q = should be approximately described 
by using in Eq. (|8.10|) the measured bulk dielectric func- 
tion e(ui). Table HT1 exhibits the relative widths Au>/ui s 
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TABLE II: Relative widths Au/uj s of surface plasmons, as 
derived from the imaginary part of the surface-response func- 
tion of Eq. (|8,10[) with meas ured values of the bulk dielectric 



function e(w) (theory) [1671 ] and from th e su r face-loss mea- 
surements at q = reported in Refs. [H, [lOl], [tol (ex- 
periment). 





Al 


Mg 


Li Na 


K 


Cs Ag 


Hg 


Theory 


0.035 


0.16 


0.33 


0.035 


0.027 


0.18 


Experiment 


0.22 


0.19 


0.35 


0.07 


0.027 


0.16 



derived in this way [167| , together with available surface- 
loss measurements at q = 0. Since silver (Ag) and 
mercury (Hg) have partially occupied d-bands, a jellium 
model, like the one leading to Eq. (|8. 10[) . is not, in prin- 
ciple, appropriate to describe these surfaces. However, 
Table [II] shows that the surface-plasmon width of these 
solid surfaces is very well described by introducing the 
measured bulk dielectric function (which includes band- 
structure effects due to the presence of (i-electrons) into 
Eq. (|8.10|) . Nevertheless, the surface-plasmon widths of 
simple metals like K and Al, with no d-electrons, are con- 
siderably larger than predicted in this simple way |284l |. 
This shows that an understanding of surface-plasmon 
broadening mechanisms requires a careful analysis of the 
actual band structure of the solid. 

Approximate treatments of the impact of the band 
structure on the surface-plasmon energy dispersion have 
been developed by several authors, but a first-principles 
description of the surface-plasmon energy dispersion 
and linewidth has been reported only in th e ca se of 
the sim ple-m etal prototype surfaces Mg(0001) |l2l| and 
Al(lll) [l22l | : these calculations will be discussed in Sec- 
tion EnTDl 



5. Multipole surface plasmons 

In his attempt to incorporate the smoothly decreasing 
electron density profile at the surface, Bennett 85[ solved 
the equations of a simple hydrodynamic model with a 
density profile which decreases linearly through the sur- 
face region and found that in addition to Ritchie's surface 
plasmon at u> ~ u> s , with a negative energy dispersion at 
low q wave vectors, there is an upper surface plasmon at 
higher energies. This is the so-called multipole surface 
plasmon, which shows a positive wave vector dispersion 
even at small q. 

The possible existence and properties of multipole sur- 
face plasmons was later investigated in the framework 
of hydrodynamical models for various cho i ces o f th e 
electron-density profile at the surface [285l . 12861 . 12871 ]. 
According to these calculations, higher multipole exci- 
tations could indeed exist, for a sufficiently diffuse sur- 
face, in addition to the usual surface plasmon at w ~ u> s . 
However, approximate quantum-mechanical RPA calcu- 
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FIG. 15: Angle-resolved high-resolution electron energy- loss 
spectra of the Al(lll) surface at various scattering angles 9 S . 
The primary beam energ y is 50 eV and the incident angle is 
8i — 45° (from Ref. [280j , used with permission) . 



lations gave no evidence for the existence of multipole 
surface plasmons 0, 90] , thereby leading to the specula- 
tion that multipole surface plasmons might be an artifact 
of the hydrodynamic approximation [287l |288| . 

The first experimental sign for the existence of mul- 
tipole s urfac e modes was established by Schwartz and 
Schaich |289j in their theoretical analysis of the photoe- 
missio n yie ld spectra that had been reported by Levinson 



et al. 125]. Later on, Dobson and Harris |290| used a 
DFT scheme to describe realistically the electron-density 
response at a jellium surface, to conclude that multi- 
pole surface plasmons should be expected to exist even 
for a high-density metal such as Al (which presents a 
considerably abrupt electron density profile at the sur- 
face). Two years later, direct experimental evidence 
of the existence of multipole surface plasmons was pre- 
sented in inelastic reflection electron scattering exper- 
iments on smo oth films of the low-density metals Na, 
K, and Cs |l00| . the intensity of these multipole surface 
plasmons being in agreeme nt w ith the DFT calculations 
reported later by Nazarov [29 1| . The Al multipole sur- 
face plasmon has been detected only recently by means of 
angle-resolved high-r esolu tion electron-energy-loss spec- 
troscopy (HREELS) (280| . 

Figure IT51 shows the loss sp ectra of the Al(lll) surface, 
as obtained by Chiarello al. [28Cl| with an incident elec- 
tron energy of 50 eV and an incident angle of 45° with 
respect to the surface normal. The loss spectrum ob- 
tained in the specular geometry (6 S = 45°) is character- 
ized mainly by a single peak at the conventional surface- 
plasmon energy u> s = 10.55 eV. For off-specular scatter- 
ing angles (9 S ^ 45°), the conventional surface plasmon 
exhibits a clear energy dispersion and two other features 
arise in the loss spectra: the multipole surface plasmon 
and the bulk plasmon. The loss spectrum obtained at 
9 S = 53°, which corresponds to q = [see Eq. (|8.9p ] is 
represented again in Fig. 1161 but now together with the 
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FIG. 16: Angle-resolved high-resolution electron energy-loss 
spectra of the Al(lll) surface at 8 S — 53°, for 9i = 45° and 
Si = 50 eV (as in Fig. I15[) but now together with the de- 
convolution into contributions corresponding to the excita- 
tion of the conventional surface plasmon at u) s = 10.55 eV, 
multipole surface plasmon at 13.20 eV, and bulk plasmon at 
ujp = 15.34 eV (from Ref. [280( |. used with permission). 



TABLE III: Anglc-rcsolvcd low-energy inelastic electron scat- 
tering measurements of the energy u) m and width Aui of mul- 
tipole surface plasmons at q = 0, as report ed in Refs. [lOOl ] 
for Na and K, in Ref. [Tol for Cs, in Ref. [To2 for Li and 
Mg, and in Ref. |280l | for Al. Also shown is the ratio uJij 
calculated in the RPA and ALDA and reported in Ref. |101| . 





r s 


Wm (eV) 




RPA 


ALDA 


Auj (eV) Alo/uj 


Al 


2.07 


13.20 


0.86 


0.821 


0.782 


2.1 0.14 


Mg 


2.66 






0.825 


0.784 




Li 


3.25 






0.833 


0.789 




Na 


3.93 


4.67 


0.81 


0.849 


0.798 


1.23 0.21 


K 


4.86 


3.20 


0.84 


0.883 


0.814 


0.68 0.18 


Cs 


5.62 


2.40 


0.83 


0.914 


0.837 


0.64 0.22 



background subst ract ion and Gaussian-fitting procedure 
reported in Ref. |280| . The peak at uj p = 15.34 corre- 
sponds to the excitation of the Al bulk plasmon, and the 
multipole surface plasmon is located at 13.20 eV. 

The calculated and measured energies and linewidths 
of long-wavelength (q — > 0) multipole surface plasmons 
in simple metals are given in Table IIIII On the whole, 
the ratio u m /w p agrees with ALDA calculations. Good 
agreement between ALDA calculations and experiment is 
also obtained for the entire dispersion of multipole sur- 
face plasmons, which is found to be approximately linear 
and positive. This positive dispersion is originated in 
the fact that the centroid of the induced electron den- 
sity, which at ui ~ ui s is located outside the jcllium edge, 
is shifted into the metal at the multipole resonance fre- 
quency uj m . We also note that multipole surface plas- 
mons have only been observed at wave vectors well be- 
low the cutoff value for Landau damping, which has been 
argued to be due to an interplay between Coulomb and 



FIG. 17: Real part of the electron density induced at ui = ui a 
(solid line) and ui = ui m (dashed line) , as reported in Ref. [l67f | 
for a model unperturbed electron density profile with linear 
decay over a finite region. 



kinetic energies |292j. 

In a semi-infinite metal consisting of an abrupt step of 
the unperturbed electron density at the surface (as in the 
hydrodynamic and specular-reflection models described 
above), surface plasmons would be localized at the sur- 
face and would propagate like plane waves, as shown in 
Fig. [31 with positive and negative surface charge regions 
alternating periodically. In the real situation in which the 
electron density decays smoothly at the surface, surface 
plasmons also propagate along the surface as illustrated 
in Fig. [31 but the finite width of the electron-density pro- 
file yields fluctuating densities with finite widths, as il- 
lustrated in Fig. [17] for a model surface electron density 
profile with linear decay over a finite region. 

Figure [L7l qualitatively represents the real situation in 
which apart from small Friedel oscillations the distribu- 
tion of the conventional surface plasmon at u s consists in 
the direction normal to the surface of a single peak (i.e., 
it has a monopole character), while the charge distribu- 
tion of the upper mode at u) m has decreasing oscillating 
amplitude towards the interior of the metal, i.e., it has 
a multipole character. Along the direction normal to the 
surface, the electronic density associated with this mul- 
tipole surface plasmon integrates to zero. 

In the retarded region, where q < lu s /c, the surface- 
plasmon dispersion curve deviates from the non-retarded 
limit (where q » uj s /c) and approaches the light line 
uj = cq, as shown by the lower solid line of Fig. [21 thus 
going to zero at q = 0; hence, in a light experiment the 
external radiation dispersion line will never intersect the 
surfacc-plasmon line, i.e., in an ideal flat surface the con- 
ventional monopole surface plasmon cannot be excited 
in a photoyield experiment. However, the multipole sur- 
face plasmon dispersion curve crosses the light line at 
q ~ uj m /c and goes to u> m at q — 0. Consequently, anglc- 
and energy-resolved photoyield experiments are suitable 
to identify the multipole surface plasmon. A large in- 
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crease in the surface photoyield was observed fr om K 
and Rb at u> m = 3.15 and 2 .84 e V, respectively |293| . 
from Al(100) at uj m = 12.5 eV [HI, and from Al(lll) at 
(jj rn = 13 eV |l26l |. in nice agreement with the multipole- 
plasmon energy ob serv ed with the HREELS technique 
by Chiarcllo et al. [280]. HREELS measurements yield, 
however, a FWHM of 2.1 e V fo r the q = multipole 
surface plasmon in Al(lll) [280j . which is considerably 
smaller than the F WHM of 3 eV measured by photoyield 
experiments jl26l |. 



D. Surface plasmons: real surfaces 

1. Stabilized jellium model 

A simple way of including approximately the lattice 
potential that is absent in the jellium model is the so- 
called s tabi li zed jellium or structureless pseudopotcntial 
model |294l . |295| |. which yields energy stability against 
changes in the background density. 

In this model, a solid surface is assumed to be trans- 
lationally invariant in the plane of the surface, as in the 
jellium model. Hence, single-particle wave functions can 
be separated as in Eq. (|7.62|) into a plane wave along the 
surface and a component ipi(z) describing motion normal 
to the surface. In the framework of DFT and TDDFT, 
this component is obtained by solving self-consistently a 
Kohn-Sham equation of the form of Eq. (|7.65|) but with 
the effective Kohn-Sham potential of Eq. (|7.66|) being re- 
placed by 

vi<s[na](z) = v H [n ](z) + v xc [n a ](z)+ < Sv >ws, 

(8.18) 

< Sv >ws representing the difference between a local 
pseudopotcntial and the jellium potential: 



< Sv > 



ws- 



3r| 
2rl 



3Z 2 / 3 
10r s ' 



(8.19) 



where Z is the chemical valence of the solid and r c is a 
core radius that is chosen to stabilize the metal for given 
values of the parameters r s and Z. 

The s tabil ized jellium model was used by Ishida and 
Liebsch |296l | to carry out RPA and ALDA calculations 
of the dispersion of the energy and linewidth of sur- 
face plasmons in Mg and Li. It is well known that the 
stabilized jellium model gives considerably better work 
functions and surface energies than the standard jellium 
model; however, the impact of a structureless pseudopo- 
tential on the properties of surface plasmons is found 
to be very small. In particular, this simple pseudopo- 
tential cannot explain the presence of core polarization 
lowering the surface-plasmon frequency for all wave vec- 
tors and does not account for the fact that the measured 
q dependence of the Li surface-plasmon energy is very 
much flatter than calculated within the jellium model 
(see Fig. [T5|) . This discrepancy (not present in the case 



TABLE IV: Measured values of the surface-plasmon energy 
Us at q — and the long-wavelength surface-plasmon disper- 
sion coefficient a of Eq. (|8.16|l for the noble metal Ag and the 
transition metals Hg and Pd. Also shown in this table is the 
coefficient a obtained from RPA and ALDA calculations of 
the centroid of the electron density induced at u) = ui 3 [see 
Eq. (|8.16p ] in a homogeneous electron gas with the electron 
density equal to that of valence sp electrons in Ag and Hg. 
Note that the measured values of the surface-plasmon energy 
uj s at q = are considerably below the jellium prediction: 
u)p/V2 = y3/2rfe 2 /ao, mainly due to the presence of d elec- 
trons. uj s and a are given in eV and A, respectively. 

r s uj s Direction 
Ag(100) 3~M 
Ag(lll) 3.02 
Ag(110) 3.02 



Hg 
Pd 



2.65 



3.7 
3.7 
3.7 

6.9 
7.4 



< 100 > 

< 110 > 



a exp 


a RPA 


_ ALDA 

a 


0.377 


-0.370 


-0.609 


0.162 


-0.370 


-0.609 


0.305 


-0.370 


-0.609 


0.114 


-0.370 


-0.609 


-0.167 


-0.32 


-0.50 


-1.02 







of Al and Mg) is attributed to the presence of an in- 
tcrband transition in Li at 3.2 eV, which is only slightly 
below the measured surface-plasmon energy at q = 
(lo s = 4.28 eV) and allows, therefore, for interference be- 
tween bulk single-p articl e and surface collective modes, 
as discussed in Ref. |296| 



2. Occupied d-bands: simple models 

Significant deviations from the dispersion of surface 
plas mons at je l lium surfa ces occur on the noble metal 
Ag fl07l . fl08l Il09l. [Hot , and the transition metals 
Hg [llll ] and Pd [297j] (see Table HVj). These deviations 
are mainly due to the presence in these metals of filled 4c? 
and bd bands, which in the case of Ag yields an anoma- 
lous positive and strongly crystal-face dependent disper- 
sion. 

a. Ag. In order to describe the observed features 



of Ag surface plasmons, Liebsch [1141 ] considered a self- 
consistent jellium model for valence 5s electrons and ac- 
counted for the presence of occupied 4<i bands via a po- 
larizable background of d electrons characterized by a lo- 
cal dielectric function ed(u j) th at can be taken from bulk 
optical data (see Fig. [TH]) [2981 ], 

For a homogeneous electron gas, one simply replaces 
in this model the bare Coulomb interaction v(r,r') by 



u'(r,r';u>) = v(r, r') e d 1 (ui), 



.20) 



which due to translational invariance yields the follow- 
ing expression for the Fourier transform W'(q,u>) of the 
screened interaction W'(r, r'; cu) of the form of Eq. (|7.3[) : 



W'(q,u) 



v(q) 



e(q,to) + e d {ijj) - 1 



.21) 



with e(q, w) being the dielectric function of a homoge- 
neous system of sp valence electrons (5s electrons in the 
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ra(eV) 



FIG. 18: The thin solid line represents the real part of the 
d band contrib ution Cd(w) to the measured optical dielectric 
function of Ag [298f | . The thick solid line represents e(qw) + 
£d(v) — 1 at q = 0, i.e., £d(a>) — uip/u 2 . The vertical dotted 
lines represent the energies at which the plasmon conditions 
of Eqs. (pH2)) and (pT29j) are fulfilled: u)' p = 3.78 eV and uj' a = 
3.62 eV, respectively. At these frequencies: £d(u' p ) = 5.65 and 
ed(wp) = 5.15. 



case of Ag). Hence, the screened interaction W'(q, cu) has 
poles at the bulk-plasmon condition 



e(q,u) + e d (u>) -1 = 0, 



(8.22) 



which in the case of Ag (r s = 2.02) and in the absence of 
d electrons (e d = 1) yields the long- wavelength (q — > 0) 
bulk plasmon energy uj p = 8.98 eV. Instead, if d elec- 
trons are characterized by the frequency-dependent di- 
electric function e d (tu) represented by a thin solid line in 
Fig. [18l Eq. (|8.22[) yields the observed long-wavelength 
bulk plasmon at 



3.78 eV. 



.23) 



In the case of a solid surface, one still uses a modified 
(d-screened) Coulomb interaction v'(r,r';u>) of the form 
of Eq. (|8.20[) . but with e d (ui) being replaced by 



e d (z,uj) = 



z < z d 



(8.24) 



1. 



z > z d , 



which represents a polarizablc background of d electrons 
that extends up to a certain plane at z = z d . Us- 
ing E q (I8.24[) . the 2D Fourier transform of v'(r, r';w) 
yields [2991] : 



2tt 



q\\ td(z',uj) 
x <r d (uj)e 



[e~ 9 » I 2 " 2 '! +sgn(z d - z') 
(8.25) 



i\\ \ z - z d\ c -q\\ \zd-z I 



where 



@ d 



[e d H-l]/[e d H + l]. 



.26) 



Introduction of Eqs. (|8.25[) and (|8 . 26[) into and equation 
of the form of Eq. (|7.59|) yields the following expression 
for the modified RPA surface-response function: 



dze9Z M^) + 

e d (zuj) 



where 8n(z;q,u>) represents the RPA induced density of 
sp valence electrons [which is given by Eq. (|7.76p with 
fxc[no]{z,z';q,uj) = 0], and 



a(q,u)=<r d (u) [e«*« 

+ / dze- qlzd - zl sgn(z d - z) 



Sn(z; q, lu) 



£d{z,uj) 



(8.28) 



In the long- wavelength (q — > 0) limit, Eq. (|8.2T|) takes 
the form of Eq. (|8.10|) but with e(co) replaced by e(ui) + 
e d (u) — 1, which leads to the surface-plasmon condition: 



e(cu)+e d (Lu) = 0. 



(8.29) 



For a Drude dielectric function e(u>) with uj p = 8.98 eV 
(r s = 3.02) and in the absence of d electrons (e d = 1), 
the surface-plasmon condition of Eq. (|8.29|) yields the 
surface-plasmon energy lo s = 6.35 eV. However, in the 
presence of d electrons characterized by the frequency- 
dependent dielectric function e d (cu), Eq. (|8.29p leads to a 
modified ((i-screened) surface plasmon at 



3.62 eV, 



(8.30) 



which for 



8.98 eV and the dielectric function e d (ui) 



represented by a thin solid line in Fig. [18] yields tu' s = 
3.62 cV, only slightly belo w the energy 3.7 cV of the mea- 
sured Ag surface plasmon [l09j ] . 

At finite wavelengths, the surface-plasmon disper- 
sion was derived by Liebsch from the peak positions of 
the imaginary part of the surface-response function of 
Eq. (|8.27[) . as obtained with a self-consistent RPA cal- 
culation of the induced density 8n{z\ q, u>) of 5s 1 valence 
electrons. The surface-plasmon dispersion was found to 
be positive for z d < and to best reproduce the observed 
linear dispersion of Ag surface plasmons for z d = —0.8 A. 
Hence, one finds that (i) the s-d screened interaction is 
responsible for the lowering of the surface-plasmon en- 
ergy at q = from the free-electron value of 6.35 eV 
to 3.62 eV, and (ii) the observed blueshift of the surface- 
plasmon frequency at increasing q can be interpreted as a 
reduction of the s-d screened interaction in the " selvedge" 
region that is due to a decreasing penetration depth of 
the induced electric field. 

With the aim of describing the strongly crystal-face 
dependence of the s urface-plasmon energy dispersion in 
Ag, Feibelman [ll3j ] thought of the Ag surface plasmon 
as a collective mode that is split off the bottom of the 
4o?-to-5s electron-hole excitation band, and considered 
the relation between the surface-plasmon dispersion and 
the Ad-to-5s excitations induced by the surface-plasmon's 
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FIG. 19: Surface-plasmon energy dispersion for the low-index 
faces (100), (111), and (110) of Ag. In the case of Ag(110), 
the exhibited surface-plasmon energy dispersions correspond 
to wave vectors along the [0011 and [110] directions, (a) Ex- 
perimental data from Ref. [l09j |. (b) The r esult s obtained 
within the jellium-dipolium model (from Ref. |ll7l ], used with 
permission). 



field. He calculated the s-d matrix elements using a 
one-dimensional surface perturbation of the jellium Lang- 
Kohn potential, and he argued that for Ag(100) the dis- 
persion coefficient is increased relative to Ag(lll) be- 
cause in the case of Ag(100) the Ad electrons lie closer to 
the centroid of the oscillating free-electron charge, and 
the probability of 4d-to-5s excitation is, therefore, en- 
hanced. 

Most recently, a self-consistent jellium model for va- 
lence 5s electrons in Ag was combined with a so-called 
dipolium model in which the occupied 4c? bands are rep- 
resented in terms of polarizable sphe res located at the 
sites of a semiinfinite fee lattice [l!7| to calculate elec- 
tron energy-loss spectra for all three low-index faces of 
Ag. The surface-plasmon energy dispersions obtained 
from these spectra are exhibited in F ig. [TTJI together with 
the experimental measurements [109]. This figure shows 
that the trend obtained for the different crystal orien- 
tations is in qualitative agreement with the data. On 
the one hand, the surface-plasmon dispersion of Ag(100) 
lies above that of the (111) surface; on the other hand, 
the slope for Ag(110) is larger for wave vectors along the 
[001] direction than when the wave vector is taken to have 
the [110] direction, illustrating the effect of the intcrpla- 
nar geometry on the effective local fields. The observed 
overall variation of the measured positive slope with crys- 
tal orientation is, however, considerably larger than pre- 
dicted by the jellium-dipolium model. This must be a 
signature of genuine band structure effects not captured 
in the present model, which calls for a first-principles ab 
initio description of the electronic response of real Ag 
surfaces. 

The impact of the band structure on the surface- 
plasmon energ y dis persion in Ag was addr essed by 
Moresco et al. [300| and by Savio et al. [30 1| . These 
authors performed test experiments with the K/Ag(110) 
and O/Ag(001) systems and concluded that (i) the 
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FIG. 20: Hg surface-plasmon energy dispersion, as obtained 
from standard jellium calculations of the RPA and ALDA sur- 
face loss function g(q, uj) of Eq. (|7.74|l (upper dashed and solid 
curves, respectively), from stabilized jellium calculations of 
the RPA and ALDA surface loss function g'(q, lo) of Eq. (|8.27|) 
(upper dashed and solid curves, respectively), and from the 
peak positions obs erved in experimental electron energy-loss 
spectra (from Ref. |lll| . used with permission). 



surface-plasmon dispersion can be modified at will by 
manipulating the surface electronic structure near u s and 
(ii) surface interband transitions between Shockley states 
should be responsible for the anomalous linear behavior 
of the surface-plasmon dispersion in Ag(100). Further- 
more, HREELS experiments on sputtered and nanostruc- 
tured Ag(100) have shown that the anomaly exhibited by 
surface plasmons on Ag(100) can inde ed be eliminated by 
modifying the surface structure [302l ]. 

b. Hg. The s-d polarization model devised and used 
by Liebsch to describe th e po sitive energy dispersion of 
surface plasmons in Ag jl!4| was als o em ployed for a 
description of surface plasmons in Hg 

At q = 0, the measured surface-plasmon energy of this 
transition metal (~ 6.9 cV) lies about 1 eV below the 
value expected for a bounded electron gas with the den- 
sity equal to the average density of 6s 2 valence electrons 
in Hg (r s — 2.65). This can be explained along the lines 
described above for Ag [see Eq. (|8.30[) ] , with the use of a 
polarizable background with ed(y' s ) = 1.6. 

At finite wave vectors, the s-d screening in Hg is found 
to considerably distort the surface-plasmon dispersion, 
as in the case of Ag, but now the linear coefficient of the 
low-g dispersion being still negative though much smaller 
than in the absence of d electrons. This is illustrated 
in Fig. [201 where the calculated (RPA and ALDA) and 
measured energy dispersions of the Hg surface plasmon 
are compared to the corresponding dispersions obtained 
within a standard jellium model (with r s — 2.65) in the 
absence of a polarizable medium. This figure clearly 
shows that in addition to the overall lowering of the 
surface-plasma frequencies by about 12% relative to the 
jellium calculations (with no d electrons) the s-d screen- 
ing leads to a significant flattening of the energy disper- 
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sion, thereby bringing the jellium calculations into nice 
agreement with experiment. 

Simple calculations have also allowed to describe cor- 
rectly the measured broadening of the Hg surface plas- 
mon. As shown in Table [Til the Hg surface-plasmon 
width at q = is very well described by simply introduc- 
ing the measured bulk dielectric function into Eq. (18.101) . 
For a descriptio n of the surface-plasmon broadening at 
finite q, Kim at introduced into Eq. (|8.10[) a Drude 
dielectric function of the form of Eq. (|2.15[) . but with a 
q-dependent finite rj defined as 

t?(?)= 7/(0) e-«°, (8.31) 

with a = 3A and 77(0) = Auj s /uj s being the measured 
relative width at q = 0. This procedure gives a surface- 
plasmon linewidth of about 1 cV for all values of q, in 
reasonable agreement with experiment. 

c. Pd. Collective excitations on transition metals 
with both sp and d-bands crossing the Fermi level are 
typically strongly damped by the presence of interband 
transitions. Pd, however, is known to support collective 
excitations that are relatively well defined. 

Surface plasmons in Pd(110) were investigated 
with angle-re solve d electron-energy-loss spectroscopy by 
Rocca et al. 297], These authors observed a prominent 
loss feature with a strongly negative linear initial disper- 
sion, which they attributed to a surface-plasmon excita- 
tion. At q = 0, they found lo s = 7.37 eV and Auj s ~ 2 eV, 
in agreement with optical data [303]. At low wave vec- 
tors from q = up to q = 0.2 A , they found a linear 
surface-plasmon dispersion of the form of Eq. (|8.15|) with 
a = 1 A. This linear dispersion is considerably stronger 
than in the case of simple metals (see Tabic |TJ) , and calls 
for a first-principles description of this material where 
both occupied and unoccupied sp and d states be treated 
on an equal footing. 

d. Multipole surface plasmons. Among the noble 
and transition metals with occupied d bands, multipole 
surface plasmons have only been observed recently in the 
case of Ag. In an improvem ent over previous HREELS 
experiments, Moresco et al. jl23l | performed ELS-LEED 
experiments with both high-momentum and high-energy 
resolution, and by substracting the data for two differ- 
ent impact energies they found a peak at 3.72 eV, which 
was interpreted to be the Ag multipole plasmon. How- 
ever, Liebsch argued that the frequency of the Ag multi- 
pole surface plasmon should be in the 6-8 eV range above 
rather than below the bulk plasma frequency, and sug- 
gested that the observed peak at 3.72 cV m ight not be 
associated with a multipole surface plasmon [l24|. 

Recently, the surface electronic structure and optical 
response of Ag have been studied on the basis of angle- 
and energy-resolved photoyield experiments [127| . In 
these experiments, the Ag multipole surface plasmon is 
observed at 3.7 eV, but no signature of the multipole 
surface plasmon is observed above the plasma frequency 
(ojp = 3.8 eV) i n dis agreement with the existing theoret- 
ical prediction [l24| . 



3. First-principles calculations 

First-principles calculations of the surface-plasmon en- 
ergy and linewidth dispersion of real solids have been 
carried out only in the case of the simple- meta l proto- 
type surfaces Mg(0001) QUI and Al(lll) [HI]. These 
calculations were performed by employing a supercell 
geometry with slabs containing 16(27) atomic layers of 
Mg(Al) separated by vacuum intervals. The matrix 
Xg g >(z, z'; q, u>) was calculated from Eq. (|7.84|) . with 
the sum over n and n' running over bands up to en- 
ergies of 30 eV above the Fermi level, the sum over k 
including 7812 points, and the single-particle orbitals 
V'k.n;^ 1 ") = i } \t,n(?\\)'4 ! i{.z) being expanded in a plane- 
wave basis set with a kinetic-energy cutoff of 12 Ry. The 
surface-response function 3 g =o,g'=o( c l7 was then cal- 
culated from Eq. (|7.86|) including ~ 300 3D reciprocal- 
lattice vectors in the evaluation of the RPA or TDDFT 
Fourier coefficients x g;g '(z,z';q, u) of Eq. (|7.83|) . Finally, 
the dispersion of the energy and linewidth of surface plas- 
mons was calculated from the maxima of the imaginary 
part of ff g =o,g'=o(qi w ) for various values of the magni- 
tude and the direction of the 2D wave vector q. 

T he a b initio calculations reported in Refs. |l2l| 
and jl22| for the surface-plasmon energy dispersion of 
Mg(0001) and Al(lll) with the 2D wave vector along 
various symmetry directions show that (i) there is almost 
perfect isotropy of the surface-plasmon energy dispersion, 
(ii) there is excellent agreement with experiment (thereby 
accurately accounting for core polarization not presented 
in jellium models), as long as the nonlocal (momentum- 
dependent) static XC local-field factor of Eq. (|5.24p is 
employed in the evaluation of the interacting density- 
response function, and (iii) if the corresponding jellium 
calculations are normalized to the measured value lj s at 
q = (as shown in Fig. Q2] by the thick solid line with 
open circles) ab initio and jellium calculations are found 
to be nearly indistinguishable. 

Ab initio calculations also show that the band struc- 
ture is of paramount importance for a correct descrip- 
tion of the surface-plasmon linewidth. First-principles 
TDDFT calculations of the Mg(0001) and Al(lll) 
surface-plasmon linewidth dispersions along various sym- 
metry directions are shown in Figs. and [2H respec- 
tively. Also shown in this figures are the experimental 
measurements reported in Refs. |l02| and [96| (stars) 
and the corresponding jellium calculations (dashed lines). 
For small 2D wave vectors, the agreement between the- 
ory and experiment is not as good as in the case of 
the surface-plasmon energy dispersion (which can be at- 
tributed to finite-size effects of the supercell geometry). 
Nevertheless, ab initio calculations for Mg(0001) (see 
Fig. l2~lj) yield a negative slope for the linewidth dis- 
persion at small q, in agreement with experiment, and 
properly account for the experimental linewidth disper- 
sion at intermediate and large wave vectors. Fig. |2"T1 also 
shows that the Mg(0001) surface-plasmon linewidth dis- 
persion depends considerably on the direction of the 2D 
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FIG. 21: Linewidth of surface plasmons in Mg(0001), as a 
function of the magnitude of the 2D wave vector q. T he filled 
(open) diamonds represent ab initio calculations [1211 ] , as ob- 
tained with q along the TM (TK) direction by employing 
the static XC local-field factor of Eq. (|5.24|l in the evaluation 
of the Fourier coefficients of the interacting density-response 
function. The solid and long dashed line represent the best 
fit of the ab initio calculations. The short dashed and dotted 
lines correspond to jellium and ID model-potential calcula- 
tions, also obtained by employing the static XC local-field 
factor of Eq. (|5.24l) . Stars represent the experimental data 
reported in Ref. [l02t|. 
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FIG. 22: Linewidth of surface plasmons in Al(lll), as a func- 
tion of the magnitude of the 2D wave vectors q. The filled 
(open) diamonds represent ab initio A LDA calculations, as 
obtained with q along the TM (FK) [L22T ]. The solid and 
long dashed line represent the best fit of the ab initio calcu- 
lations. The dotted line corresponds to ID model-potential 
calculations, also obtained in the ALDA. Stars represent the 
experimen tal d ata reported in Ref. [9(|. The measured values 
from Ref. |280| ] are represented by dots with error bars. 
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FIG. 23: Schematic representation of the surface band struc- 
ture on Cu(lll) near the T point. The shaded region repre- 
sents the projection of the bulk bands. 

wave vector, and that the use of the paramet erize d one- 
dimensional model potential reported in Ref. |26l| (dot- 
ted line) docs not improve the jellium calculations. As 
for Al(lll), we note that band-structure effects bring the 
jellium calculations closer to experiment; however, ab ini- 
tio calculations are still in considerable disagre emen t with 
the measurements reported in Refs. [9|| and [280j |. espe- 
cially at the lowest wave vectors. At the moment it is 
not clear the origin of so large discrepancy between the 
theory and experiment. From the other hand, consid- 
erable disagreement is also found between the measured 
surface-plasmon line width at q ~ reported in Refs. [96| 
(Aw s ~ 1.9 eV) and (HO] (Aw s - 2.3 eV) and the value 
derived (see Table [TTJ) from the imaginary part of the 
surface-response function of Eq. (|8.10|) with measured 
values of the bulk dielectric function (Aw s = 0.38 eV). 



E. Acoustic surface plasmons 

A variety of metal surfaces, such as Be(0001) and the 
(111) surfaces of the noble metals Cu, Ag, and Au, arc 
known to support a partially occupied band of Shock- 
ley surface states within a wide energy gap around the 
Fermi level (see Fig. l23f [304l l305| . Since these states are 
strongly localized near the surface and disperse with mo- 
mentum parallel to the surface, they can be considered 
to form a quasi 2D surface-state band with a 2D Fermi 
energy e'jP equal to the surface-state binding energy at 
the f point (see Table Ki- 
ln the absence of the 3D substrate, a Shocklcy surface 
state would support a 2D collective oscillation, the energy 
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3D substrate 



TABLE V: 2D Fermi energy (e 2 F D ) of surface states at the f 
point of Be(0001) and the (111) surfaces of the noble metals 
Cu, Ag, and Au. VjP and m 2D represent the corresponding 
2D Fermi velocity and effective mass, respectively, vjp is 
expressed in units of the Bohr velocity no = e 2 /%. 





e¥(eV) 


2D 1 
V F /Vo 


m 2D 


Be(0001) 


2.75 


0.41 


1.18 


Cu(lll) 


0.44 


0.28 


0.42 


Ag(lll) 


0.065 


0.11 


0.44 


Au(lll) 


0.48 


0.35 


0.28 



of the corresponding plasmon being given by [17 



.32) 



with n 2U being the 2D density of occupied surface states, 
i.e, 



i 2D = ejP/ir. 



.33) 



Eq. (|8.32p shows that at very long wavelengths plas- 
mons in a 2D electron gas have low energies; however, 
they do not affect e-ph interaction and phonon dynamics 
near the Fermi level, due to their square-root dependence 
on the wave vector. Much more effective than ordinary 
2D plasmons in mediating, e.g., superconductivity would 
be the so-called acoustic plasmons with sound-like long- 
wavelength dispersion. 

Recently, it has been shown that in the presence of 
the 3D substrate the dynamical screening at the surface 
provides a mechanism for the existence of a new acous- 
tic collective mode, the so-called acoustic surface plas- 
mon, whose ener gy exhibi t s a l inea r dependence on the 
2D wave number [l28l . 1 129L Il30l . Il3l| . This novel surface- 
plasmon mode has been observed at the (0001) surface 
of Be, showing a linear energy dispersion that is in ve ry 
good agreement with first-principles calculations |l32| . 



1. A simple model 

First of all, we consider a simplified model in which 
surface-state electrons comprise a 2D electron gas at z = 
z d (see Fig. |2"4"|) , while all other states of the semi- infinite 
metal comprise a 3D substrate at z < represented by 
the Drude dielectric function of Eq. (|2.15[) . Within this 
model, one finds that both e-h and collective excitations 
occurring within the 2D gas can be described with the 
use of an effective 2 D di electric function, which in the 
RPA takes the form [ml 



'lff(l^) = 1 ~ W(z d ,z d ;q,uj)x2D(<l, UJ ), 



.34) 



W(z, z'\ q, u>) being the screened interaction of Eq. (17.591) . 
and X2£>(g,^) being the noninteracting den sity-r esponse 
function of a homogeneous 2D electron gas jl78| . 

In the absence of the 3D substrate, W(z,z';q, ui) re- 
duces to the bare Coulomb interaction v(z,z';q), and 




z<0 



z=0 



z>0 



FIG. 24: Surface-state electrons comprise a 2D sheet of in- 
teracting free electrons at z — zj,. All other states of the 
semi-infinite metal comprise a plane-bounded 3D electron gas 
at z < 0. The metal surface is located at z — 0. 



uS) coincides, therefore, with the RPA dielec- 
tric function of a 2D electron gas, which in the long- 
wavelength (g — > 0) limit has one single zero correspond- 
ing to collective excitations at ui = LU2D- 

In the presence of a 3D substrate, the long-wavelength 
limit of e^j(q,u>) is found to have two zeros. One zero 
corresponds to a high-frequency (u » VFq) oscillation 
of energy 



.35) 



in which 2D and 3D electrons oscillate in phase with one 
another. The other zero corresponds to a low- frequency 
acoustic oscillation in which both 2D and 3D electrons 
oscillate out of phase; for this zero to be present, the long- 
wavelength (q — > 0) limit of the low-frequency (u> — > 0) 
screened interaction W(z d , Zd', q, u>), 



I(z d ) = lim W(z d ,z d ;q,u; -> 0), 

9^0 



(8.36) 



must be different from zero. The energy of this low- 
frequency mode is then found to be of the form [1 3 lj 



to = a Vp D q, 



where 



ir[ir + 2I{z d )Y 



.37) 



(8.38) 



a. Local 3D response. If one characterizes the 3D 
substrate by a Drude dielectric function of the form of 
Eq. (|2.15p . the 3D screened interaction W(z d} z d ; q, avpq) 
is easily found to be 



W(z d ,z d ;q,ctv F q) 







4irz d 



z d <0 
z d >0 



(8.39) 



Hence, in the presence of a 3D substrate that is spa- 
tially separated from the 2D sheet (z d > 0), introduction 
of I(zd) = 47rz d into Eq. (|8.38|) yields at z d >> 1: 

a = (8.40) 



which is the result first obtained by Chaplik in his study 
of cha rge-c arrier crystallization in low-density inversion 
layers [306l |. 

If the 2D sheet is located inside the 3D substrate 
( z d < 0), I{zd) — 0, which means that the effective di- 
electric function of Eq. (|8.34p has no zero at low energies 
(u> < lo s ). Hence, in a local picture of the 3D response 
the characteristic collective oscillations of the 2D electron 
gas would be completely screened by the surrounding 3D 
substrate, and no low-energy acoustic mode would exist. 
This result had suggested over the years that acoustic 
plasmons should only exist in the case of spatially sepa- 
rated plasm as, as pointed out by Das Sarma and Mad- 
hukar [3071 ]. 

Nevertheless, Silkin et al. |l28| have demonstrated 
that metal surfaces where a partially occupied quasi-2D 
surface-state band coexists in the same region of space 
with the underlying 3D continuum can indeed support a 
well-defined acoustic surface plasmon. This acoustic col- 
lective oscillation has been found to appear as the result 
of a combination of the nonlocality of the 3D dynamical 
screening an d the spill out of the 3D electron density into 
the vacuum [131 1 . 

b. Self- consistent 3D response. Figures [25k and 1 2 5b 
exhibit self-consistent RPA calculations of the effective 
dielectric function of Eq. (|8.34|) corresponding to the 
(0001) surface of Be, with q = 0.01 1 (Fig. d5k) and 
q = 0.1 Gq 1 (Fig. I2"5b) . The real and imaginary parts 
of the effective dielectric function (thick and thin solid 
lines, respectively) have been displayed for Zd = 0, as ap- 
proximately occurs with the quasi-2D surface-state band 
in Be(0001). Also shown in these figures (dotted lines) 
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lium edge (zd — 0), and outside the metal at Zd — Ay/2 
and Zd = Xf- 

Collective excitations 



zero of 
is small, 

which yields a maximum in the energy-loss function 



are related to a 

2D 



Reeeff(q,u)) in a region where Im e^fJq, 



Im 



C 2D 



0?X 



For the 2D electron density un- 
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der study (r 2D 
stratc a 2D plasmon would occur at u>2D 



1.22 cV for 



q = O.Ola^ and oj 2 d = 3.99 cV for q = O.Ioq . How- 
ever, in the presence of the 3D substrate (see Fig. [25]) 
a well-defined low-energy acoustic plasmon occurs at en- 
ergies (above the upper edge uj u = vpq + q 2 /2 of the 
2D e-h pair continuum) dictated by Eq. (|8.37[) with the 
sound velocity cxvf being just over the 2D Fermi velocity 
vf (at > 1). When the 2D sheet is located far inside the 
metal surface, the sound velocity is found to approach the 
Fermi velocity (a — ► 1). When the 2D sheet is located 
far outside the metal surface, the coefficient a approaches 
the classical limit (a — ► y/2zd) of Eq. (|8.40|) . 

Finally, we note that apart from the limiting case Zd = 
Xp and q = 0.01 Oq 1 (which yields a plasmon linewidth 
negligibly small) the small width of the plasmon peak 
exhibited in Fig.[2S]is entirely due to plasmon decay into 
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FIG. 25: Effective dielectric function of a 2D sheet that is lo- 
cated at the jellhim edge (za = 0), as obtained from Eq. (|8.34[) 
and self-consistent RPA calculations of the 3D screened inter- 
action W (zd, Zd', q,oj) and the 2D density-res pons e function 
X 2D (q,u) with (a) q = 0.01 and (b) q = 0.1 [nj]. The real 
and imaginary parts of e e ff(q,u>) are represented by thick 
and thin solid lines, respectively. The dotted lines repre- 
sent the effective 2D energy-loss function Im [— e~}f (q, w)l for 
Zd = — Xf, Zd — Xf/2, and Zd — Xf- The vertical dashed line 
represents the upper edge u> u = VFq + q 2 /2m of the 2D e- 
h pair continuum, where 2D e-h pairs can be excited. The 
calculations presented here for Zd = Xp and q — 0.01 a^ 1 
have been carried out by replacing the energy wbya complex 
quantity u) + in with n = 0.05 eV. All remaining calculations 
have been carried out for real frequencies, i.e., with r\ = 0. 
The 2D and 3D electron-density parameters have been taken 
to be r\ D = 3.14 and r s = 1.87, corresponding to Be(0001). 



c-h pairs of the 3D substrate. 
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FIG. 26: Energy loss function lm[g(q, w)]/u) of Be(0001) ver- 
sus the excitation energy u), obtained fro m Eq . (|7,73p with 
the use of t he ID model potential of Ref. [26 ll ] , as reported 
in Ref. [l28t |. The magnitude of the wave vector q has been 
taken to be 0.05 (top panel), 0.1 (middle panel), and 0.15 
(bottom panel), in units of the inverse Bohr radius a^ 1 . In 
the long- wavelength limit (q — > 0), g(q,u>) is simply the total 
electron density induced by the potential of Eq. (|7.75[) . 



2. ID model calculations 

For a more realistic (but still simplified) description of 
electronic excitations at the (0001) surface of Be and the 
(111) surface of the noble metals C u, A g, and Au, the 
ID mod el po t enti al vmp(z) of Ref. [26l| was employed 
in Refs. [l28l Il29j |. The use of this model potential al- 
lows to assume translational invariance in the plane of 
the surface and to trace, therefore, the presence of col- 
lective excitations to the peaks of the imaginary part of 
the jellium-likc surface response function of Eq. (|7.73[) . 
The important difference between the screened interac- 
tion W(z, z';q,Lu) used here to evaluate g(q,oj) and the 
screened interaction used in Eq. (|8.34[) in the framework 
of the simple model described above lies in the fact that 
the single-particle orbitals ipi( z ) an d energies are now 
obtained by solving Eq. (|7.65p with the jellium Kohn- 
Sham potential fi<"s[no](z) replaced by the model poten- 
tial vmp(z). Since this model potential has been shown 
to reproduce the key features of the surface band struc- 
ture and, in particular, the presence of a Shockley surface 
state within an energy gap around the Fermi level of the 
materials under study, it provides a realistic description 
of surface-state electrons moving in the presence of the 
3D substrate. 

Figure [55] shows the imaginary part of the surface- 
response function of Eq. (|7.73|) . as obtained for Be(0001) 
and for increasing values of q by using the ID model 
potential vmp{z) of Ref. |26l| . As follows from the fig- 
ure, the excitation spectra is clearly dominated by two 
distinct features: (i) the conventional surface plasmon 
at htu s ~ 13 eV, which can be traced to the character- 
istic pole that the surface-response function g(q,w) of a 
bounded 3D free electron gas with r^ D = 1.87 exhibits at 
this energy |308l | , and (ii) a well-defined low-energy peak 




Energy (eV) 

FIG. 27: Energy loss function Im[g{q,u>)}/u) of Be(0001) ver- 
sus the excitation energy uj obtained fro m Eq . (|7.73p with the 
use of the ID m odel potential of Ref. 261] and for various 
values of q |310l |. 
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FIG. 28: Energy-loss function lm\g(q, uj) ]/ui of the (111) sur- 
faces of the noble metals Cu, Ag, and Au 129], shown by solid, 
dashed, and dashed-dotted lines, respectively, versus the ex- 
citation energy uj, as obtained f rom Eq. (|7.73|l with the use of 
the ID model potential of Ref. [26 1| and for q — 0.01 a^ 1 and 
?7 = 1 meV. The vertical solid lines are located at the energies 
uj — Vp D q, which would correspond to Eq. (|8.37|l with a = 1. 



with linear dispersion. 

That the low-frequency mode that is visible in Fig. [26] 
has linear dispersion is clearly shown in Fig. 1271 where the 
imaginary part of the surface-response function g{q,uj) 
of Be(0001) is displayed at low energies for increasing 
values of the magnitude of the wave vector in the range 
q = 0.01 — 0.12 Og . The excitation spectra is indeed 
dominated at low energies by a well-defined acoustic peak 
at energies of the form of Eq. (|8.37[) with an a coefficient 
that is close to unity, i.e., the sound velocity being at 
long wavelengths very close to the 2D Fermi velocity vjP 
(see Table E| . 

The energy-loss function lmg(q,oj) of the (111) sur- 
faces of the noble metals Cu, Ag, and Au is displayed 
in Fig. [55] for q = 0.01 . This figure shows again the 
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FIG. 29: The solid line sho ws t he energy of the acoustic sur- 
face plasmon of Be(0001) [128| |. as obtained from the max- 
ima of the calculated surface-loss function Img(ij,ii;) shown 
in Fig. I27I The thick dotted line and the open circles repre- 
sent the maxima of the energy-loss function Im [— 1 /e e / / (q, u>)] 
obtained from Eq. (|8.34[) with Zd far inside the solid (thick 
dotted line) and with Zd — (open circles). The dashed line 
is the plasmon dispersion of a 2D electron gas in the absence 
of the 3D system. The gray area indicates the region of the 
(q,u>) plane (with the upper limit at u^, = vjfq + q 2 /2rri2D) 
where e-h pairs can be created within the 2D Shockley band of 
Be(0001). The area below the thick solid line corresponds to 
the region of momentum space where transitions between 3D 
and 2D states cannot occur. The quantities w™™ r and q"" n 
are determined from the surface band structure of Be(0001). 
2D and 3D electron densities have been taken to be those cor- 
responding to the Wigner radii rl D = 3.14 and r j. D = 1.87, 
respectively. 



- L - L . L . 

FIG. 30: A periodic grating of constant L. The grating peri- 
odic structure can provide an impigning free electromagnetic 
radiation with additional momentum 2n/L. 



presence of a well-defined low-energy collective excitation 
whose energy is of the form of Eq. (|8.37p with a ~ 1. 

Figurel251shows the energy of the acoustic surface plas- 
mon of Be(0001) versus q (solid line), as derived from 
the maxima of the calculated lmg(q, u>) of Fig. [27] (solid 
line) and from the maxima of the effective energy-loss 
function Im [— l/e e ff(q, ui)] (see Fig. obtained from 
Eq. (|8.34[) for Zd = (dashed line). Little discrepancies 
between these two calculations should be originated in (i) 
the absence in the simplified model leading to Eq. (|8.34|) 
of transitions between 2D and 3D states, and (ii) the 
nature of the decay and penetration of the surface-state 
orbitals, which in the framework of the model leading to 
Eq. (|8.34[) are assumed to be fully localized in a 2D sheet 
at z = Zd- 



3. First-principles calculations 

First-principles calculations of the imaginary part of 
the surface-response function g g =n s =p(q_, ui) of Be(0001) 
have been carried out recently [l32|, and it has been 
found that this metal surface is indeed expected to sup- 
port an acoustic surface plasmon whose energy dispersion 
agrees with the solid line represented in Fig. [55] (if the 
dispersion of Fig. [55] calculated for surface state effective 
mass m = 1 i s scaled according to the ab initio value 
m — 1.2 309]). Furthermore, these calculations have 
been found to agree closely with recent high-resol ution 
EELS measurements on the (0001) surface of Be jl32| 
(under grazing incidence), which represent the first evi- 
dence of the existence of acoustic surface plasmons. 

4- Excitation of acoustic surface plasmons 

As in the case of the conventional surface plasmon 
at the Ritchie's frequency u> s , acoustic surface plasmons 
should be expected to be excited not only by moving 
electrons (as occurs in the EELS experiments reported 
recently |132| ) but also by light. Now we focus on a 
possible mechanism that would lead to the excitation of 
acoustic surface plasm ons b y light in, e.g., vicinal sur- 
faces with high indices |3 1 Of ] . 

At long wavelengths (q — > 0), the acoustic surface- 
plasmon dispersion curve is of the form of Eq. 1|8.37[) with 
a ~ 1. As the 2D Fermi velocity vjP is typically about 
three orders of magnitude smaller than the velocity of 
light, there is, in principle, no way that incident light 
can provide an ideal surface with the correct amount of 
momentum and energy for the excitation of an acous- 
tic surface plasmon to occur. As in the case of conven- 
tional surface plasmons, however, a periodic corrugation 
or grating in the metal surface should be able to provide 
the missing momentum. 

Let us consider a periodic grating of constant L (see 
Fig. [3D]). If light hits such a surface, the grating periodic 
structure can provide the impigning free electromagnetic 
waves with additional momentum arising from the grat- 
ing periodic structure. If free electromagnetic radiation 
hits the grating at an angle 9, its wave vector along the 
grating surface has magnitude 

q=-sm8±—n, (8.41) 
c L 

where L represents the grating constant, and n = 1, 2, . . .. 
Hence, the linear (nearly vertical) dispersion relation of 
free light changes into a set of parallel straight lines, 
which can match the acoustic-plasmon dispersion rela- 
tion as shown in Fig. 1311 

For a well-defined acoustic surface plasmon in Be(0001) 
to be observed, the wave numb er q needs to be smaller 
tha n q ~ 0.06 Op 1 (see Fig. [55]) [31l| . For q = 0.05 a^ 1 , 
Eq. (|8.41|) with n = 1 yields a grating constant L = 66 A. 
Acoustic surface plasmons of energy uj ~ 0.6 eV could be 
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FIG. 31: A schematic representation of the dispersion relation 
of acoustic surface plasmons (solid line) and free light impign- 
ing on a periodic grating of constant L (essentially vertical 
dotted lines). 



excited in this way. Although a grating period of a few 
nanometers sounds unrealistic with present technology, 
the possible control of vicinal surfaces with high indices 
could provide appropriate grating periods in the near fu- 
ture. 



microsco 
nology [3 




IX. APPLICATIONS 



Surface plasmons have been employed over the years 
in a wide spectrum of studies ranging from condensed 
matter and surface physics |p, fiol ITU H3 . [H, [3, 
EE Eip, H SIS M, [llMJplto electrochem- 
istrv [25l] wetting |26l | b iosensing |27l l28l, [2^ | , scanning 
tunne ling microscopy [3fJ, the ejection of ions from sur- 
faces |3l| nanoparticle growth [32|, [H[ , surface-plasmon 
3al , and surface-plasmon resonance tech- 
H HI, HO, EH, \M- Renewed interest in 
surface plasmons has come from recent advances in the 
investigation of the optical properties of nanostructured 
materials [HI, Hi| , one of the most attractive aspects of 
these collective excitations now being their use to concen- 
trate light in subwavclcngth structures and to enhance 
transmission through periodic arrays of subwavelcngth 
holes in optically thick metallic films [45|, |46[, as well 
as the possible fabrication of nanoscale photonic circuits 
operating at optical frequencies [48| and their use as me- 
diators in the transfer of energy from donor to acceptor 
molecules on opposite sides of metal films (67j . 

Here we focus on two distinct applications of collective 
electronic excitations at metal surfaces: The role that 
surface plasmons play in particle-surface interactions and 
the new emerging field called plasmonics. 



A. Particle-surface interactions: energy loss 

Let us consider a recoilless fast point particle of charge 
Zi moving in an arbitrary inhomogeneous many-electron 
system at a given impact vector b with nonrclativistic ve- 
locity v, for whic h ret ardation effects and radiation losses 
can be neglected |312j . Using Fermi's golden rule of time- 
dependent perturbation theory, the lowest-order proba- 
bility for the probe particle to transfer moment um q to 
the medium is given by the following expression [313J: 



47T 9 

P « = ~LA 



dq' 



^ib-(q+q' 



J (27T) 3 

x JmW(q, q';w) 5(w-q- v) 8(w + q' ■ v), (9.1) 

where L and A represent the normalization length 
and area, respectively, and W(q,q';uj) is the double 
Fourier transform of the screened interaction W(r, r';u>) 
of Eq. (HI: 

W(q,q';u) = j dr J dr' e- 4(q r+q ' r ' ) lF(r, r'; w). 

(9.2) 

Alternatively, the total decay rate r 1 of the probe 
particle can be obtained from the knowledge of the imagi- 
nary part of the sel f-ene rgy. In the GW approximation of 
many-body theory [314| , and replacing the probe-particle 
Green function by th at of a non-interacting recoilless par- 
ticle, one finds [263| : 

r- 1 = -2 Z\ I dr j dr' #(r) cf>* f (r') 
f J 

xlmW(r, r', e z -e f ) &(r') ^(r), (9.3) 

where 4n{r) represents the probe-particle initial state of 
energy £,;, and the sum is extended over a complete set 
of final states / (r) of energy e / . Describing the probe- 
particle initial and final states by plane waves in the di- 
rection of motion and a Dirac <5 function in the transverse 
direction, i.e., 



1 



e lv ' r v^(r±-b), 



(9.4) 



where rj^ represents the position vector perpendicular to 
the projectile velocity, one finds that the decay rate of 
Eq. (|9.3[) reduces indeed to a sum over the probability 
P q of Eq. (OH]), i-e.: 



1 = -Tp 



(9.5) 



T being a normalization time. 

For a description of the total energy AE that the 
moving probe particle loses due to electronic excitations 
in the medium, one can first define the time-dependent 
probe-particle charge density 



p e * t (r,t) = Z 1 6(r-h-vt), 



(9.6) 
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and one then obtains the energy that this classical parti- 
cle loses per unit time as follows [88| 

where V md (r,t) is the potential induced by the probe 
particle at position r and time t, which to first order in 
the external perturbation yields [see Eq. (|7.f [) ]: 



p f + OQ p 

V md (r,t) = J dr' J dt' J 



du iul ( t -t') 
2ir 



x W(r,r';u)p ext (r',t' 



with 



W(r, r'; u) = W(r, r'; u) - v(r, r'), 
Finally, one writes: 



AE 



+ 00 



dt 



(9.8) 



(9.9) 



(9.10) 



dE 
~dt 



Introducing Eqs. QHJ) and (|978[) into Eq. (j9J|) . and 
Eq. (|9~7) into Eq. (PTTUj) . one finds that the total en- 
ergy loss AE can indeed be written as a sum over the 
probability P q of Eq. (|9~Tj) . i.e.: 



A£ = ]T(q.v)P q , 



(9.11) 



where q • v is simply the energy transferred by our recoil- 
less probe particle to the medium. 



1. Planar surface 

In the case of a plane-bounded electron gas that is 
translationally invariant in two directions, which we take 
to be normal to the z axis, Eqs. (|9.6[) - (|9.8|) yield the fol- 
lowing expression for the energy that the probe particle 
loses per unit time: 



dE . Z\ 

dt 7T 



d 2 q 

(2tt)= 



dt' I duo (J 
la 



-i("-q-vii)(t-t') W[ z (t), z(t');g,o>],(9.12) 



where q is a 2D wave vector in the plane of the surface, 
vy represents the component of the velocity that is par- 
allel to the surface, z(t) represents the position of the 
projectile relative to the surface, and W(z, z'\ q, ui) is the 
2D Fourier transform of W(r, r'; u>). 

In the simplest possible model of a bounded semi- 
infinite electron gas in vacuum, in which the screened 
interaction W(z, z'; q,w) is given by the classical expres- 
sion Eq. (|7.8p with ei being the Drude dielectric function 
of Eq. (|2.15p and e% = 1, explicit expressions can be found 
for the energy lost per unit path length by probe parti- 
cles that move along a trajectory that is either parallel 



or normal to the surface. 




FIG. 32: Particle of charge Z\ moving with constant velocity 
at a fixed distance z from the surface of a plane-bounded 
electron gas. Inside the solid, the presence of the surface 
causes (i) a decrease of loss at the bulk-plasmon frequency 
lj p varying with z as ~ Ko(2u p z/v) and (ii) an additional 
loss at the surface-plasmon frequency u s varying with z as 
~ Ko(2ui s z/v). Outside the solid, energy losses are dominated 
by a surface-plasmon excitation at u) s . 



a. Parallel trajectory. In the case of a probe particle 
moving with constant velocity at a fixed distance z from 
the surface (see Fig. [32l . introduction of Eq. (|7.8| into 
Eq. (PT2"j) yields [31of~ 



dE 
dx 



■1 1 — 



Up [ln(k c v/ujp) - K (2uj p z/v)} + u 2 s K (2u> s z/v), z < 



u 2 s K (2u s \z\/v), 



z>0, 



(9.13) 



wher e K o(a) is the zero-order modified Bessel func- able. 

tiongll, and k c denotes the magnitude of a wave vector For particlc trajectories outside the solid (z > 
above which long-lived bulk plasmons arc not sustain- ^ Eq ^ reproduces the classical expression of 
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Echcnique and Pendry |22j , w hich was found to describe 
correctly EELS experiments [l35j and whic h was ex- 
tended to include relativistic corrections [315| . For par- 
ticle trajectories inside the solid [z > 0), Eq. (I9.13 p re- 
produces the result first obtained by Nunez et al. ]316t 
Outside the solid, the energy loss is dominated by the 
excitation of surface plasmons at lu s . When the parti- 
cle moves inside the solid, the effect of the boundary is 
to cause (i) a decrease in loss at the bulk plasma fre- 
quency Lo p , which in an infinite electron gas would be 
—dE/dx = Z 2 u; 2 \n(k c v / ! uo p ) / 'v 2 and (ii) an additional loss 
at the surface-plasma frequency to s . 

Nonlocal effects that are absent in the classical 
Eq. (|9.13|) were incorporated approximately by several 
authors in the framework of the hydrodynamic ap- 
proach and the specular-re flection model described in 
Sections lWBl1 and lVlTB"2l (3171 [HI, [HI [32l|, HH . 
More recently, extensive RPA and ALDA calculations of 
the energy-loss spectra of charged parti cles moving near 
a jellium surface were carried out [257( within the self- 
consistent scheme described in Section fVII B 31 At high 
velocities (of a few Bohr units) and for charged particles 
moving far from the surface into the vacuum, the actual 
energy loss was found to converge with the classical limit 
dictated by the first line of Eq. (|9. 13|) . However, at low 
and intermediate velocities substantial changes in the en- 
ergy loss were observed as a realistic description of the 
surface response was considered. 

Corrections to the energy loss of charged particles 
(moving far from the surface into the vacuum) due to 
the finite width of the surface-plasmon resonance that is 
not present, in principle, in jellium s elf-co nsistent calcu- 
lations, have been discussed recently 323]. These correc- 
tions have been included to investigate the energy loss 
of highly charged ions undergoing distant collisions at 
grazing incidence angles with the internal surface of mi- 
crocapillary materials, and it has been suggested that the 
correlation between the angular distribution and the en- 
ergy loss of transmitted ions can be used to probe the 
dielectric properties of the capillary material. 

For a more realistic description of the energy loss of 
charged particles moving near a Cu(lll) surface, the 
Kohn-Sham potential vks Iz) u sed in the self-consi stent 
jellium calculations of Ref. [257j | was replace d in Ref . 324 
by the ID model potential vmp of Ref. [26l| . It was 



shown, however, that although the Cu(lll) surface ex- 
hibits a wide band gap around the Fermi level and a 
well-defined Shockley surface state the energy loss ex- 
pected from this model does not differ significantly from 
its jellium counterpart. This is due to the fact that the 
presence of the surface state compensates the reduction 
of the energy loss due to the band gap. 



Existing first-principles calculations of the interaction 
of charged particles with solids invoke periodicity of the 
solid in all directions and neglect, therefore, surface ef- 
fects and, in p articular, the excitation of surfaces plas- 
mons |325l 13261 ] . An exception is a recent first-principles 
calculation of the ene rgy l oss of ions moving parallel with 
a Mg(0001) surface [3271 ] . which accounts naturally for 
the finite width of the surface-plasmon resonance that 
is present nei ther in the self-consistent jellium calcula- 
tions of Ref. |257| nor in the ID model calculations of 
Ref. 



324] 



A typical situation in which charged particles can be 
approximately assumed to move along a trajectory that is 
parallel to a solid surface occurs in the glancing-incidcncc 
geometry, where ions penetrate into the solid, they skim 
the outermost layer of the solid, and are then specularly 
repelled by a repulsi ve, sc reened Coulomb potential, as 
discussed by Gemell [328| . By first calculating the ion 
trajectory under the combined influence of the repulsive 
planar potential and the attractive image potential, the 
total energy loss can be obtained approximately as fol- 
lows 



AE = 2v 



dE 
dx 



(z) [v z (z)] 1 dz, 



(9.14) 



Zt p and v z (z) denoting the turning point and the value of 
the component of the velocity normal to the surface, re- 
spectively, which both depend on the angle of incidence. 

Accurate measurements of the energy loss of ions being 
reflected from a variety of solid surfaces at gr azing inci- 
den c e have be en reported by several autho rs 13291 133C 



33l | . In particular, Winter et al. [331| carried 
out measurements of the energy loss of protons being re- 
flected from Al(lll). From the analysis of their data at 
120 keV, these authors deduced the energy loss dE/dx(z) 
and found that at large distances from the surface the en- 
ergy loss follows closely the energy loss expected from the 
excitation of surface plasmons. Later on, RPA jellium 
calculations of the energy loss from the excitation of va- 
lence electrons were combined with a first-Born calcula- 
tion of the energy loss due to the excitation of the inner 
shells and reasonably good agreement with the e xperi - 
mental data was obtained for all angles of incidence [3341 ] . 

b. Normal trajectory. Let us now consider a situa- 
tion in which the probe particle moves along a normal 
trajectory from the vacuum side of the surface (z > 0) 
and enters the solid at z = t = 0. The position of the 
projectile relative to the surface is then z(t) = —vt. As- 
suming that the electron gas at z < can be described 
by the Drude dielectric function of E q. (I2.15[) and intro- 
ducing Eq. (O into Eq. O yields [310j: 



dE 

dz 



r?2 I LO 2 [ln(k c v/ujp) — h(uj p z/v)] + uj 2 h(uj s z/v), z<0 
u; 2 f(2u; s \z\/v), z > 0, 



(9.15) 
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which would exist in its absence, as noted by Ritchie 
A more accurate jellium self-consistent description of the 
energy loss of charged particles passing through thi n foil s 
has been performed recently in the RPA and ALDA [335| | . 



B. STEM: Valence EELS 



FIG. 33: Particle of charge Z\ passing perpendicularly 
through a finite foil of thickness a. The presence of the bound- 
aries leads to a decrease in the energy loss at the bulk-plasmon 
frequency u v and an additional loss at the surface-plasmon 
frequency ui 3 , the net boundary effect being an increase in 
the total energy loss in comparison to the case of a particle 
moving in an infinite medium with no boundaries. 



where 



h(a) = 2cos(a)f(a)-f(2a), 



(9.16) 



with f(a) being given by the following expression: 



/(«) = 



xe 



1 



-dx. 



(9.17) 



Eq. (|9.15|) shows that (i) when the probe particle is mov- 
ing outside the solid the effect of the boundary is to cause 
energy loss at the surface-plasmon energy lo s , and (ii) 
when the probe particle is moving inside the solid the ef- 
fect of the boundary is to cause both a decrease in loss at 
the bulk-plasmon energy lo p and an additional loss at the 
surface-plasmon energy u s , as predicted by Ritchie 

Now wc consider the real situation in which a fast 
charged particle passes through a finite foil of thickness 
a (see Fig. 133]) . Assuming that the foil is thick enough 
for the effect of each boundary to be the same as in the 
case of a semi-infinite medium, and integrating along the 
whole trajectory from minus to plus infinity, one finds 
the total energy that the probe particle loses to collec- 
tive excitations: 



7 2 

ae = % 

v 



a to p hi- 



ke 



7T 

2 Uf 



(9.18) 



This is the result first derived by Ritchie in a different 
way [l|, which brought him to the realization that sur- 
face collective excitations exist at the lowered frequency 
co s . The first term of Eq. (|9.18j) . which is proportional 
to the thickness of the film represents the bulk contribu- 
tion, which would also be present in the absence of the 
boundaries. The second and third terms, which are both 
due to the presence of the boundaries and become more 
important as the foil thickness decreases, represent the 
decrease in the energy loss at the plasma frequency u) p 
and the energy loss at the lowered frequency u) s , respec- 
tively. Eq. (|9.18|) also shows that the net boundary effect 
is an increase in the total energy loss above the value 



The excitation of both surface plasmons on solid sur- 
faces and localized Mie plasmons on small particles has 
attracted great interest over the years i n the fields of 
scanning tran smission electron microscopy fl3l. lT35l.fl36l 
Il37l Il38l Il39j and near-field optical spectroscopy [l4fl jT^ 

EELS of fast electrons in STEM shows two types of 
losses, depending on the nature of the excitations that 
are produced in the sample: atomically defined core- 
electron excitations at energies uj > 100 eV and valence- 
electron (mainly collective) excitations at energies up to 
~ 50 eV. Core-electron excitations occur when the probe 
moves across the target, and provide chemic al in forma- 
tion about atomic-size regions of the target [336| . Con- 
versely, valence-electron excitations provide information 
about the surface structure with a resolution of the order 
of several nanometers. One advantage of valence EELS is 
that it provides a strong signal, even for non-penetrating 
trajectori es ( t he s o-called aloof beam energy loss spec- 
troscopy |337l . l338| h and generates less specimen dam- 
age |l35| . 

The central quantity in the interpretation of valence 
EELS experiments is the total probability P(to) for the 
STEM beam to exchange energy u with the sample. In 
terms of the screened interaction W(r, r';cu) and for a 
probe electron in the state </>o(r) with energy £o, first- 
order perturbation theory yields: 

PM = ~2j2j dr J ^>}(r')0o(r')^/(r)^(r) 



ImW^r, r'; w) 5{u — £q — £f] 



(9.19) 



where the sum is extended over a complete set of final 
states <^/(r) of energy Ef. For probe electrons moving 
on a definite trajectory (and having, therefore, a charge 
density of the form of Eq. (|9.6p with no beam recoil), 
the total energy loss AE of Eqs. (|9.10[) and (|9.11[) can be 
expressed in the expected form 



AE 



dw lo P(w). 



(9.20) 



Ritchie and Howie [339j | showed that in EELS experi- 
ments where all of the inelastic scattering is collected 
treating the fast electrons as a classical charge of the 
form of Eq. (|9.6p is indeed adequate. Nonetheless, quan- 
tal effects due to the spatial extens ion of th e bea m have 
been addressed by several authors [l37l l34dl . l34l[ | . 
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Planar surface 



2. Spheres 



In the case of a classical beam of electrons moving with 
constant velocity at a fixed distance z from a planar sur- 
face (see Fig. |3"2")) . the initial and final states can be de- 
scribed by taking a S function in the transverse direction 
and plane waves in the direction of motion. Neglecting 
recoil, the probability P(u) of Eq. (|9.19p then takes the 
following form: 



p(u) 



dq x ImW{z, z;q,u), 



(9.21) 



with q = \/q^ + {u/v) 2 . 

In a classical model in which the screened interac- 
tion W(z, z';q,u) is given by the classical expression 
Eq. (|7.8[) . the probability P(u) is easy to calculate. In 
particular, if the beam of electrons is moving outside the 
sample, one finds: 



W(z, z; q, u>) 



2tt 



[1 



-2qz 



(9.22) 



with the surface-response function g given by Eq. (|7.9p . 
Introduction of Eq. (|9~2"2"|) into Eq. (f9~2Tj) yields the clas- 
sical probability 



PM 



2L 



Kq(2ujz/v) Img(u), 



(9.23) 



which in the case of a Drude metal in vacuum yields 



P(u) = L -4-S(u - u 8 ), 



(9.24) 



and, therefore [see Eq. (|9.20p ]. the energy loss per unit 
path length given by Eq. (|9.13[) with Z\ = — 1 and z > 0. 

The classical Eq. (|9.22[) can be easily extended to the 
case in which the sample is formed by a semi-infinite 
medium characterized by e\ and covered by a layer of 
dielectric function £3 and thickness a. One finds: 



2ir 

W(z,z;q,u) = — 

q 



-2q Z £32 +gi3e 
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where 



and 



,32 



£1 



£3 ~ £2 
£3 + £2 



£1 - e 3 
ei + £3 ' 



(9.25) 



(9.26) 



(9.27) 



Eqs. (|9.25j) - ()9.27j) show that while for a clean surface 
(a = 0) the energy-loss function P(u) is dominated by the 
excitation of surface plasmons at u = u s [sec Eq. (|9.24|) ]. 
EELS should be sensitive to the presence of sub -surf ace 
structures. This effect was observed by Batson jl34l | in 
the energy-loss spectra corresponding to an Al surface 
coated with an AI2O3 layer of increasing thickness. 



a. Definite trajectories. In the case of a classical 
beam of electrons moving with constant velocity at a 
fixed distance b from the center of a single sphere of ra- 
dius a and dielectric function £i(w) that is immersed in 
a host medium of dielectric function £2(0;), the initial 
and final states can be described (as in the case of the 
planar surface) by taking a S function in the transverse 
direction and plane waves in the direction of motion. If 
recoil is neglected and the classical screened interaction 
of Eqs. f7TT|) - ([7T3|) is used, then Eq. (j9T9]) yields the 
following expression for the energy-loss probability [l36| : 



OO I 



ttv* * — ' * — ' (l — m)(l + m) 

1=1 m=0 V ;V ' 
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x K^(u b /v)lmgi(uj), 
for outside trajectories (a < b), and [34 



(9.28) 



P H = E E » 



(I — m) 



(I + m) 
1 1 

im] ' 

for inside trajectories (a > 6), with 



Al 



1=1 m=0 

x lm 9l (u) + A\ m [A° lm + A\ m ] Im [e^H] } , (9.29) 
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dx 
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Here, c = y/aF—W, P[ a {x) denote Legendre func- 
tions, Ci m {x) = cos(x), if (I + m) is even, os sin(x), if 
(I + m) is odd, and gi(u) denotes the classical function 
of Eq. (|7.13[) . which has poles at the classical surfacc- 
plasmon condition of Eq. (|3.33p . 

Equations (|9.28|) and (|9.29|) show that an infinite num- 
ber of multipolar modes can be excited, in general, which 
contribute to the energy loss of moving electrons. It is 
known, however, that for small spheres with a << v/u s 
the dipolar mode (1 = 1) dominates, which in the case of 
a Drude sphere in vacuum occurs at Ui = u p For 
spheres with a ~ v /u s , many multipoles contribute with 
similar weight. For very large spheres (a >> v/u s ) and 
a < b, the main contribution arises from high multipo- 
lar modes occurring approximately at the planar surface 
plasmon energy uj Si since in the limit of large I gi(u) of 
Eq. (|7.13p reduces to the energy- loss function of Eq. (|7.9p . 
Indeed, this is an expected result, due to the fact that for 
very large spheres the probe electron effectively interacts 
with an almost planar surface. 
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b. Broad beam. We now consider a broad beam ge- 
ometry, and we therefore describe the probe electron 
states by plane waves of the form: 



and 



ikn-r 



(9.32) 



electron gas but with the inverse dielectric function 
e~ 1 (Q 1 ui) replaced by the effective inverse dielectric func- 
tion e-j f {Q,w) of Eq. liTTHjl . 



3. Cylinders 



(9.33) 



where f2 represents the normalization volume. Then, 
intr oducin g Eqs. ([9732]) . (|9~33"1) and l[7JT|) - ([7l3|l into 
Eq. (|9.19|1 . and neglecting recoil, one finds 

P(w) = i J ^ica[-e-} f {Q,w)]8{u-Cl-v). (9.34) 

This is precisely the energy-loss probability correspond- 
ing to a probe electron moving in a 3D homogeneous 



a. Definite trajectory. In the case of a classical beam 
of electrons moving with constant velocity at a fixed dis- 
tance b from the axis of a single cylinder of radius a 
and dielectric function £i(uj) that is immersed in a host 
medium of dielectric function e 2 (u)), the initial and fi- 
nal states can be described (as in the case of the planar 
surface) by taking a S function in the transverse direc- 
tion and plane waves in the direction of motion. If re- 
coil is neglected and the classical screened interaction of 
Eqs. dZHD-dZ^lD is used, then Eq. (PH?)) yields the fol- 
lowing expression for the energy-loss probability [l39j : 
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(9.35) 



(9.36) 



for outside trajectories (a > b). In particular, for axial 
trajectories (b = 0), the energy-loss probability P(cu) is 
due exclusively to the m = mode. 

b. Broad beam. For a broad beam geometry the 
probe electron states can be described by plane waves of 
the form of Eqs. (|9.32|) and (|9.33p . which after introduc- 
tion into Eq. (|9.19[) . neglecting recoil, and using the clas- 
sical screened interaction of Eqs. (|7.19[) - (|7.22[) yield an 
energy-loss probability of the form of Eq. (|9.34[) but now 
with the inverse dielectric function e^ 1 (Q,u}) replaced 
by the effective inverse dielectric function e~^(Q,uj) of 
Eq. (TT251 . 



C. Plasmonics 

Renewed interest in surface plasmons has come from 
recent advances in the investigation of the electromag- 
netic properties of nanostructured materials (43l . [44| , one 
of the most attractive aspects of these collective excita- 
tions now being their use to concentrate light in subwavc- 



lcngth structures and to enhance transmission through 
periodic arrays of subwavelength holes in optically thick 
metallic films (4f| H^]. Surface-plasmon polaritons are 
tightly bound to metal-dielectric interfaces penetrating 
around 10 nm into the metal (the so-called skin-depth) 
and typically more than 100 nm into the dielectric (de- 
pending on the wavelength), as shown in Fig. 1121 In- 
deed, surface plasmons of an optical wavelength concen- 
trate light in a region that is considerably smaller than 
their wavelength, a feature that suggests the possibility 
of using surface-plasmon polaritons for the fabrication of 
nanoscale photonic circuits operating at optical frequen- 
cies [47l.l4"8j. 

Here we discuss only a few aspects of the most recent 
research that has been carried out in the so-called field 
of plasmonics. This constitutes an important area of re- 
search, since surface-plasmon based circuits are known 
to merge the fields of photonics and electronics at the 
nanoscale, thereby enabling to overcome the existing dif- 
ficulties related to the large size mismatch between the 
micrometer-scale bulky components of photonics and the 
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nanometer-scale electronic chips. Indeed, the surface- 
plasmon polariton described in Sec. |TT] can serve as a 
base for constructing nano-circuits that will be able to 
carry optical signals and electric currents. These opto- 
electronic circuits would consist of various components 
such as couples, waveguides, switches, and modulators. 

In order to transmit optical signals to nanophotonic 
devices and to efficiently increase the optical far-field to 
near-field conversion, a nanodot coupler (fabricated from 
a linear array of closely spaced metallic nanoparticles) 
has been combined recently with a surfacc-plasmon po- 
lariton condenser (working as a phase array) fabricated 
from hemispherical metallic nanoparticles [49|. By fo- 
cusing surface-plasmon polaritons with a spot size as 
small as 400 nm at A = 785 nm, their transmission length 
through the nanodot coupler was confirmed to be three 
times longer than that of a metallic-core waveguide, ow- 
ing to the efficient near-field coupling between the local- 
ized surface plasmon of neighboring nanoparticles. 

Achieving control of the light-solid interactions in- 
volved in nanophotonic devices requires structures that 
guide electromagnetic energy with a lateral mode confine- 
ment below the diffraction limit of light. It was suggested 
that this so-called subwavclcngth-sized wave guiding can 
occur along chains of closely spaced metal nanoparticles 
that convert the optical mode into surface (Mie) plas- 
mons [5C|. The existence of guided long-range surface- 
plasmon waves was observed experimentally by using 

U 



closely 



thin metal films 1511 52 . 153j , nanowires 
spaced silver rods [56 1, and metal nanoparticles [571 ]. 

The propagation of guided surface plasmons is subject 
to significant ohmic losses that limit the maximum prop- 
agation length. In order to avoid these losses, various 
geometries have been devised using arrays of features of 
nanosize dimensions [58l [59j . The longest propagation 
length (13.6 mm) has been achieved with a structure 
consisting of a thin lossy metal film lying on a dielec- 
tric substrate and covered by a different dielectric super- 
strate [60] . The main issue in this context is to strongly 
confine the surface-plasmon field in the cross section per- 
pendicular to the surface-plasmon propagation direction, 
while keeping relatively low propagation losses. Recently, 
it has been pointed out that strongly localized channel 
plasmon polaritons (radiation waves guided by a chan- 
nel cut into a planar surface of a solid characterized by 
a negative dielectric function [6l| ) exhibit relatively low 
propagation losses [62]. The first realization and char- 
acterization of the propagation of channel plasmon po- 
laritons along straight subwavelength metal grooves was 
reported by Bozhevilnyi et al. [631 ]. More recently, the 
design, fabrication, and characterization of channel plas- 
mon polariton based subwavelength waveguide compo- 
nents have been reported; these are Y-splittcrs, Mach- 
Zehnder interferometers, and waveguide-ring resonators 
operating at telecom wavelengths [64| . 

In the framework of plasmonics, modulators and 
switches have also been investigated. Switches should 
serve as an active element to control surface-plasmon po- 



lariton waves [65j, [66j. This approach takes advantage 
of the strong dependence of the propagation of surface- 
plasmon polaritons on the dielectric properties of the 
metal in a thin surface layer that may be manipulated 
using light. This idea was realized introducing a few- 
micron longgallium switching section to a gold-on-silica 
waveguide [661 ] . An example of an active plasmonic device 
has been demonstrated by using a thin silver film covered 
from both sides by thin polymer films with molecular 
chromophores [63]. In this case, coupled surface-plasmon 
polaritons provide an effective transfer of excitation en- 
ergy from donor molecules to acceptor molecules on the 
opposite sides of a metal film up to 120 nm thick. 

Another emerging area of active research in the field of 
plasmonics is based on the generation and manipulation 
of electromagnetic radiation of various wavelengths from 
microwave to optical frequencies. For instance, the coat- 
ing of semiconductor quantum wells by nanometer metal 
films results in an increased spontaneous emission rate in 
the semiconductor that leads to the enhancement of light 
emission. This enhancement is due to an efficient en- 
ergy transfer from electron-hole pair recombination in the 
quantum well to surface-plasmon polaritons at the sur- 
face of semiconductor heterostructures coated by metal. 
Recently, a 32-fold increase in the spontaneous emission 
rate at 440 nm in an InGaN/GaN quantum well has 
been pro bed by time-resolved photolumincsccncc spec- 
troscopy [3431 ] . Also probed has been the enhancement of 
photoluminescence up to an order of magnitude through 
a thin metal film from organic light emitting diodes, by 
removing the surfacc-plasmon polar iton quenching with 
the use of a periodic nanostructure 344]]. 

Recent theoretical and experimental work also sug- 
gest that surface-plasmon polaritons play a key role both 
in the transmission of electromagnetic waves through a 
single aperture and the enhanced trans miss i on o f ligh t 
through subwavelength hole arrays [45|, 345, 346, 13471 ], 
This enhanced transmission has also been observed at 
millim eter- waves and micro- waves [348l |349| and at THz- 
waves [3501 . 



Finally, we note that surface-plasmon polaritons have 
been used in the field of nanolithography. This surfacc- 
plasmon based nanolithography can produce subwave- 
length structures at surfaces, such as sub-100 nm lines 
with visible light [HI SH, HI SH . On the other hand, 



the enhancement of evanescent waves thro ugh the excita- 
tion of surface plasmons l ed P endry [355| to the concept 
of the so-called superlens [3561 ] . 
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